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SUMMARY 


This report is a supplemental NASTRAN document for a new 
capability to determine the vibratory response of turbosystems 
subjected to aerodynamic excitation. Supplements to NASTRAN 
Theoretical, User's, Programmer's, and Demonstration Manuals are 
included . 

Turbosystems such as advanced turbopropellers with highly 
swept blades, and axial-flow compressors and turbines can be 
analyzed using this capability, which has been developed and 
implemented in the April 1984 release of the general purpose 
finite element program NASTRAN. 

The dynamic response problem is addressed in terms of the 
normal modal coordinates of these tuned rotating cyclic 
structures. Both rigid and flexible hubs/disks are considered. 
Coriolis and centripetal accelerations, as well as differential 
stiffness effects are included. 

Generally non-uniform steady inflow fields and uniform flow 
fields arbitrarily inclined at small angles with respect to the 
axis of rotation of the turbosystem are considered as the sources 
of aerodynamic excitation. The spatial non-uniformities are 
considered to be small deviations from a principally uniform 
inflow. Subsonic and supersonic relative inflows are addressed, 
with provision for linearly interpolating transonic airloads. 

A stand-alone program, AIRLOADS, has been additionally 



developed to generate the vibratory airloads on the blades of the 
turbosystems • This program can be used as a pre— processor to 
NASTRAN. 

Both NASTRAN and pre-processor capabilities are operational on 
the CRAY 1-s computer system at NASA's Lewis Research Center. 

The work was conducted under Contract NAS3-24387 from NASA 
LeRC, Cleveland, Ohio, with Mrs. Marsha Nall as the Program 
Monitor. 
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SECTION 1 


SUPPLEMENT TO 
NASTRAN THEORETICAL MANUAL 
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MODAL FORCED VIBRATION ANALYSIS OP 
AERODYNAMICALLY EXCITED TURBOSYSTEMS 


1.1 INTRODUCTION 

A new capability has been developed and implemented in NASTRAN 
to perform modal forced vibration analysis of turbosystems 
subjected to aerodynamic excitation. 

Single- and counter-rotating advanced turboprops with 
significantly swept blades, and axial-flow compressors and 

turbines are examples of turbosystems that can be analyzed using 
this capability. 


Generally non-uniform steady inflow fields and uniform flow 
fields arbitrarily inclined at small angles with respect to the 
axis of rotation of the turbosystem are considered as the 
aerodynamic sources of excitation. Subsonic and supersonic 
relative inflows are recognized, with a provision for linearly 
interpolating transonic aerodynamic loads. 

The capability has been implemented in the April 1984 release 

of NASTRAN on the CRAY 1-s computer system at NASA's Lewis 
Research Center. 

Highlights of the theoretical basis of the new capability 
(Ref. 1) are presented in this section. Details of the User's, 
Programmer's, and Demonstration Manuals are contained in Sections 
2, 3, and 4, respectively. 
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1.2 PROBLEM DESCRIPTION 


Figure 1.1 shows a single-rotation advanced turboprop^ as an 
example of turbosystems ^ operating in a generally non— uniform 
steady inflow field. 

Although the absolute inflow field does not change with time, 
the rotation of the turboprop results in velocities with 
oscillatory components relative to the blades. Relative 
velocities with harmonic components at the rotational frequency 
also exist in uniform flow fields when the turboprop axis of 
rotation is misaligned with the absolute flow direction. 

Given such operating conditions, it is desired to, 

1. determine the oscillatory loading distributions over the 

blades of the turboprop at various excitation frequencies, 
and 

2, determine the resulting vibratory response (displacements, 
stresses, etc.) of the turboprop. 

As per NASA's needs, the problem of determining the applied 
oscillatory airloads on the turbosystem blades has been addressed 
in a stand-alone development outside, and independent of, NASTRAN 
(Ref. 2). The stand-alone program, AIRLOADS, however, can also 
function as a pre-processor to NASTRAN analyses. 

Determination of the modal vibratory response is discussed 
further in the following subsections. 
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1.3 THEORY 


1.3.1 Structural Aspects 

1. Structures of turbosystems are treated as tuned cyclic 
structures with identical mass^ stiffness/ dampin 9 / nnd 
constraint properties for all cyclic sectors. The 
structural modelling capabilities of NASTRAN for 
rotationally cyclic structures are fully admitted. 

2. Turbosystems with both rigid and flexible hubs/disks are 
considered. 

3. Differential stiffness effects due to centrifugal loads 

and any (externally specified) steady state airloads are 
included. 

4. Coriolis and centripetal acceleration (stiffness 
softening) effects are taken into account. 

5. Circumferential harmonic-dependent normal modes of tuned 
cyclic structures are used in formulating and solving the 
dynamic response problem in the frequency domain. 

1.3.2 Ae rodynamic Aspects 

1. Aerodynamic modelling is essentially dictated by the 
unsteady aerodynamic theories used to determine the 
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unsteady blade loading distribution. Due to the use of 
two-dimensional cascade aerodynamic theories, the blade 
aerodynamic model comprises a series of chordwise strips 

stacked spanwise to cover the entire blade surface as shown 
in Figure 1.2. 

2. Two-dimensional subsonic and supersonic cascade aerodynamic 
routines are utilized for generating the reactionary 
airloads on turbosystem blades due to oscillatory blade 
motions. Blade sweep effects are Included in both 
routines. Transonic airloads are linearly interpolated. 

3. Externally specified aerodynamic loads can be applied to 

any degree of freedom of the structural model. These 

degrees of freedom are not restricted to those used in 

generating reactionary airloads mentioned in point (2) 
above . 


Equations of Motion and their Solution 

The total translational and rotational displacements at any 
fixed point of the rotating and vibrating turbosystem, expressed 
in body-fixed coordinate systems, consist of 

steady state components due to the steady airloads and 
centrifugal loads, and 

vibratory components due to the vibratory excitation, 
superposed on the steady displacements. 
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The aerodynamic vibratory response problem of the turbosystem is 
treated herein in terms of the vibratory components of total 
displacements • 


For an K-bladed tuned turbosystem, with structural coupling 
between blades via a relatively flexible hub, the equations of 
forced motion can be written as (Ref, 1) 


- [:<■’] 






and 
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The forcing term on the right hand side of equation (1) is 
entirely due to aerodynamic excitation. Cyclic sector numbers 

and their sides referred to in equation (2) are illustrated in 
Figure 1.3 . 


In seeking solutions for the vibratory displacements uT , 
for all n , based on the qualitative and quantitative nature of 
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the right hand side forcing functions, the following steps, 
applicable to tuned cyclic structures with rotational cyclic 
symmetry, are considered. 


1. The displacements u" (and the loads P” ) can be written 
as (Ref. 1) 






M-/ N/i 

+ (-1) ]. _ 


(3) 


This IS a standing wave representation wherein the 

coefficients and 

^ ^ ^ are, in general, 

functions of time or frequency. For a given circumferential 
harmonic index k, by defining appropriate relations between 
U and u. ® the above equation can be transformed to a 

travelling wave form ( ' Ref . 1 ) . The constants k^ and 

are given by. 
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2. For a given circumferential harmonic index 
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3. With the use of inter-segment compatibility constraint 
conditions (equation 2) 




(«) 

(7) 


-K 


where u. is an independent displacement vector consisting 
of LL and U ^ degrees of freedom from the interior and 
side 1 of a cyclic sector. The transformations ^ . and 

Cr 

functions of the circumferential harmonic k, and 
express the side 2 degrees of freedom in terms of those on 
side 1. 


4. Equation (5) can then be written as 


{u.' I = ko.; 4- ka.) 


(8) 


5. The real eigenvalue problem, for a given k, can then be 
stated from equation (1) as 


C9; 


where 
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The eigenvectors j^ j* J obtained from the solution of 

equation (9), can be used to introduce the modal coordinates 
^ as 




00 


7. The modal equations of forced motion of the turbosystem, for 

a given circumferential harmonic index k, can then be written 
from equation (1) 

- _ (n; 

where 




{13) coy^td. 
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[ 1 ""] = 




t^cfe3 CQ 3C^ck.3 CQ”3L^_5t3 




8. For a given circumferential harmonic index k, P and 

— ks 

P are the circumferential harmonic components of the 
total external excitation. Such excitation due to 
aerodynamic sources is discussed in R^fs. 1 and "27^^ 

9. Tl^ generalized oscillatory aerodynamic reaction matrix 

"ZZ ^ 

Q can be written as 

[Q j = [gi 3Cq3C«* J+C?* 3 Cgi 11^ 3 0+; 


where 


"1 = 3 


(jB’) Covy-td, 
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are the 'cosine' and 'sine' component mode shapes of 
eigenvectors ^ Equation (14) is rewritten as 


— K 

] = 




0 &) 


For turbosystem structures with flexible hub/disk, 

^*1 «*ist when the circumferential harmonic 
index k ^ 0 and :jb n/ 2 when N, the total number of cyclic 
segments in the structure/ is even/ and 

b) only exists when k = 0 or N/2, n even. 

For turbosystem structures with rigid hub/disk, each cyclic 
segment of the structure behaves structurally independent of 
its adjacent segments. Degrees of freedom at segment 
boundaries are completely constrained to zero. The only 
possible structural modes are those akin to k = 0 modes with 
fixed inter-segment boundaries. Only exists. 

Derivation of is discussed in detail in Ref: 1. ' 
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Equation U2) can now be solved for Substitution in 

equation (11), and equation (8) yields 

11 . Repeating steps 2 through 10 for all applicable 

circumferential harmonic indices, and substitution in 
equation (3), result in a” for all n 

12 . Other dynamic responses such as stresses, etc., can be 

obtained for all sectors of the turbosystem by current 
NASTJRAN procedures* 

The procedure described above for the solution of 

circumferential harmonic components of dynamic response, for a 

given circumferential harmonic index, has been implemented in the 

April 1984 release of NASTRAN on the CRAY 1-s computer system at 
NASA LeRC. 

An overall flowchart of the solution procedure is shown in 
Figure 1.4 , 

Supplements to the NASTRAN User's, Programmer's, and 
Demonstration manuals are presented in the following sections. 
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Figure 1 .2 NASTRAN Aerodynamic Model 
of Turboprop Blade for 2-D 
Cascade Theories {Ref. 1) 
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SECTION 2 


SUPPLEMENT TO 
NASTRAN USER'S MANUAL 
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MODAL FORCED VIBRATION ANALYSIS OF 
AERODYNAMICALLY EXCITED TURBOSYSTEMS 

2.1 INTRODUCTION 

All aspects of conducting modal forced vibration analysis of 
aerodynamically excited turbosystems rotating about their axis of 
symmetry are presented from a NASTRAN user *s viewpoint in this 
section. 

Generally stated, the complete problem of modal forced 
vibration analysis of aerodynamically excited turbosystems can be 
considered in four distinct phases (Figure 2.1); 

Phase 1 to conduct differential stiffness analysis. 

Phase 2 to investigate aeroelastic stability at the operating 
point. 

Phase 3 to generate the applied oscillatory airloads at the 
operating condition, and 

Phase 4 to finally conduct the modal forced vibration 
.analysis. 

All of the four phases are illustrated by the example 
presented in the Demonstration manual supplement. Section 4. 

Phase 1 uses DISP APP Rigid Format 4 with ALTERS to save the 
total stiffness matrix. 

Phase 2 uses AERO APP Rigid Format 9 (Refs. 3 and 4). 
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Phase 3 uses the pre-processor program AIRLOADS (Ref. 2). This 
phase can be conducted using any alternative procedure for 
airloads generation. 

Phase 4 uses DISP APP Rigid Format 8 with the MFVAAET Alter 
package. This phase is discussed in detail in the following 
subsections. 

2.2 NASTRAN MODEL 

The user models one rotationally cyclic sector (segment) of 
the entire structure as shown in Figure 2,2. This modelled sector 
is considered as the n = 1 st, sector. Each cyclic sector is 
defined with two sides which identify its boundaries with the two 
adjacent cyclic sectors (Figure 1,3) . 

The side 2 degrees of freedom are related to those on side 1 
via the circumferential harmonic index. The modelling of rigid 
hub/disk conditions is accomplished by completely constraining 
all degrees of freedom on both sides of the cyclic sector to 
zero. Although the circumferential harmonic index is irrelevant 
in such situations, it should be selected as zero for 
computational efficiency. 

The structural model is prepared using the general 
capabilities of NASTRAN for modelling rotationally cyclic 
structures . 

The aerodynamic model for the generation of reactionary 
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airloads comprises a grid defined by the intersection of a series 
of chords and "computing stations" (Figure 1.2). The chords are 
selected normal to any spanwlse reference curve such as the blade 
leading edge. The choice of the number and location of the 
chords and the computing stations is dictated by the expected 
variation of the relative flow properties across the blade span, 
and the complexity of the mode shapes exhibited by the propeller 
blade. Due to its resemblance to the structural model of the 
blade, and the adequacy of a relatively coarse grid to describe 
the spanwise flow variations, the aerodynamic model is chosen as 
a subset of the structural model as shown in Figure 1.2. 

The aerodynamic grid is specified on STREAMLl bulk data 
cards . 


2.2«1 Coordinate Systems 

In order to conveniently pose and solve the aetodynamically 
forced vibration problem of turbosystems, a number of coordinate 
systems have been defined. Figure 2.2 illustrates these 
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coordinate systems for an advanced turbopropeller with its axis 


of rotation mounted at an angle with respect to the tunnel mean 
flow. 

Each of these coordinate systems is described as follows: 

~ Tunnel coordinate system 

* This is defined to conveniently specify the velocity 
components of the spatially non-uniform tunnel free 
stream. It can be suitably oriented based on the 
available tunnel data. In the special case of 
^^^o^ynamic excitation in uniform inflow^ the tunnel 
coordinate system is oriented such that the X^Z_^ plane 
is parallel to the X^Z^ plane of the inertial 
coordinate system as shown in Figure 2.3. The origin of 
the X^y^Z^ system is arbitrarily located. The 
inclination angle of the turbosystem axis of rotation 
with respect to the tunnel flow also lies in a plane 

parallel to X^Z^ plane. The uniform flow is directed 
along +X^ axis . 

- X^Y^Z^ Inertial coordinate system 

In the present problem, this coordinate system is used 
to relate the quantities in the tunnel and the basic 
coordinate systems. The orientation of this coordinate 
system is completely arbitrary except for the X^ axis 
to be parallel to, and in the direction of, X^, axis of 
the basic coordinate system described next. The zero 
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reference for time/phase measurements is defined when 
the inertial and the basic coordinate systems are 
parallel • 

All of the following NASTRAN coordinate systems are fixed to 
the rotating turbosystem. 

- Basic coordinate system 

* This coordinate system has its axis coincident with 
the turbosystem axis of rotation / and directed 
aftward. Location of the origin is arbitrary. The 

plane contains (approximately) the maximum 
planform of the modelled blade. The definition of this 
coordinate system is consistent with the theoretical 
developments of the 2— d cascade unsteady aerodynamics 
presently incorporated in the Bladed Disks Computer 
Program (Ref. 1 ) . 

“ (Blade) shank-fixed coordinate system 

* The principal advantage of this shank-fixed coordinate 
system is in modelling changes in the blade setting 
angles by a simple 3x3 transformation matrix relating 
to the basic coordinate system. coincides with the 
blade shank axis. The definition of the coordinate 
system otherwise is arbitrary. 

Grid point location and displacement coordinate 

systems 
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* Any number of such rectangular, cylindrical, or 
spherical coordinate systems can be completely 
arbitrarily defined to locate grid points of the 
NASTRAN model, as well as request output at these grid 
points. All of the coordinate systems used for 

output requests collectively form the NASTRAN global 
coordinates system. 


Internally generated coordinate system on swept 

chord s 

* This coordinate system is generated within the present 
Bladed Disks Computer Program, and is used to define 
flow and motion properties for the unsteady aerodynamic 
theories on a given swept chord s. It is located at 
the blade leading edge with the directed aftward 

along the chord s. y^ is defined normal to the blade 
local mean surface. 
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2.3 EXECUTIVE CONTROL DECK 


The salient points are noted as follows: 

1. The NASTRAN card is required immediately preceding the ID 
card in the Executive Control Deck, and must contain, at 
least, the following operational parameter: 

NASTRAN SYSTEM (93) « 1 

This invokes the sweep effects in subsonic and supersonic 
reactionary aerodynamic routines, and is suggested for use 
even when sweep effects are negligible. In all cases where 
S TREAML2 bulk data cards are obtained from AIRLOADS program, 
this card is required. 

2. SOL 8 and APR DISP must be selected. 

3. The alter package, MFVAAET , (Modal Forced Vibration ^alysis 
of Aerodynamically Excited Turbosystems) must be included. 
The READFILE capability of NASTRAN can be utilized as 
follows: 

READFILE MFVAAET 


2.4 CASE CONTROL DECK 
2.4.1 Subcase Definitions 

The PARAMeter KMAX ( 0, NSEGS/2 for even NSEGS, 

^ (NSEGS-D/2 for odd NSEGS) determines the number. 
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order and meaning of subcases as follows: 

The number of subcases is equal to FKMAX, where 


FKMAX 


1, if KMAX 
1 + 2*KMAX, 
1 + 2*KMAX, 
and 


if 

if 


0 1 

0<KMAX<(NSEGS-I)/2, NSEGS odd 
0<KMAX ^(NSEGS-2)/2 NSEGS even 


9 

9 


NSEGS, if KMAX - NSEGS/2, NSEGS even 


SUBCASE 1 Ck' = 0) 
SUBCASE 2 Ck' = ic) 
SUBCASE 3 Ck' = Is) 
SUBCASE 4 Ck' = 2c) 
SUBCASE 5 Ck' = 2s) 


SUBCASE FKMAX ('k' - KMAXs) . 

in the event that NSEGS is even and KMM - NSEGS/2, Subcase 
FKMAX will represent 'k' . KMAXc as KMAXs does not exist. 

Circumferential harmonic components of directly applied loads 
are specified under the appropriate subcases, with RLOADi bulk 
data cards, null loads need not be specified by the user. 


Other Data Selection Items 

1. The SPC and Wc request must appear above the subcase level 
and may not be changed, 

2. me thod must be used to select an EIGR bulk data card. 

3. fr equency must be selected and must be above the subcase 
level . 
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^^Q^ENCY nuisb be used to select one and only one FREQ^ FREQl 


or FREQ2 card from the Bulk Data deck. 

5. Direct input matrices are not allowed. 

6. OFREQ must not be used. 

7. pLOAD must be used to define a frequency-dependent loading 
condition for each subcase. For frequency-dependent loads , 
subcases without loads need not refer to a PLOAD card. 

The following printed output, sorted by frequency (SORTl) or 
by point number or element number (S0RT2) , is available, either 
as real and imaginary parts or magnitude and phase angle (0* - 
360 lead) , for the list of frequencies specified: 

1. Displacements, velocities, and accelerations for a list of 
PHYSICAL points (grid points and extra scalar points 
introduced for dynamic analysis) or SOLUTION points (points 
used in formulation of the general K system) . 

2. Nonzero components of the applied load vector and 

single-point forces of constraint for a list of PHYSICAL 
points. 

3. Stresses and forces in selected elements (ALL available only 
for SORTl) . 

The following plotter output is available for Frequency 
Response calculations: 

1. Undeformed plot of the structural model. 
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2. x-Y plot Of any component of displacement, velocity, or 
acceleration of a PHYSICAL point or SOLUTION point. 

3. X-Y plot of any component of the applied load vector or 
single-point force of constraint, 

4. X-Y plot of any stress or force component for an element. 

The data used for preparing X-Y plots may be punched or 
printed in tabular form . Also, a printed summary is prepared for 

each X Y plot which includes the maximum and minimum values of 
the plotted function. 

The following items relate to Bulk Data restrictions; 

!• SUPORT cards are not allowed. 

2. jPOINT cards are not allowed. 

3. SPQINT cards are not allowed. 

4. CYJOIN cards are required. 


The following parameters are used in Modal Forced Vibration 

Analysis of Aerodynamically Excited Turbosystems: 

1. Gg pPNT - optional - A positive integer value of this 
parameter will cause the Grid Point Weight Generator to be 
executed and the resulting weight and balance information to 
be printed. All fluid related masses are ignored, 

2. j gMASS - optional - The terms of the structural mass matrix 
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are multiplied by the real value of this parameter when they 
are generated in EMA. Not recommended for use in hydroelastic 
problems. 

3. COUPMASS - fixed - Only lumped mass matrices must be used. 

4. GRAB optional - The BCD value of this parameter is used to 
tell the GKAD module the desired form of matrices KDD, BDD 
and HDD. The BCD value can be FREQRESP or TRANRESP. The 
default is TRANRESP. 

NO^: Remember to define parameters G, W3 and W4. See Section 
9.3.3 (DIRECT DYNAMIC MATRIX ASSEMBLY) Pages 9.3-7 and 9.3-8 
of the NASTRAN Theoretical manual for further details. 

5. LGKAD - optional - The integer value of this parameter is 
used in conjunction with parameter GKAD. If GKAD « FREQRESP 
then set LGKAD = 1, if GKAD = TRANRESP set LGKAD = -1. The 
default value is -1. 

optional - The real value of this parameter is used as a 
uniform structural damping coefficient in the formulation of 
dynamics problems. Not recommended for use in hydroelastic 
problems (use GE on MATl) . 

7. - optional - The real value of this parameter is used as a 

pivotal frequency for uniform structural damping if parameter 
GKAD « TRANRESP. In this case W3 is required if uniform 
structural damping is desired. The default value is 0.0. 

” optional - The real value of this parameter is used as a 
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pivotal frequency for element structural damping if parameter 
GKAD = TRANRESP. In this case W3 is required if structural 
damping is desired for any of the structural elements. The 
default value is 0.0. 

9. NSEGS - required - The integer value of this parameter is the 
number of identical segments in the structural model. 

10. CYCIO - required - The integer value of this parameter 
specifies the form of the input and output data. A value of 
+1 is used to specify physical segment representation ^ and a 
value of -1 for cyclic transform representation. The value 
of CYCIO must be set to -1. 

11. CYCSEQ - fixed - The integer value of this parameter 
specifies the procedure for sequencing the equations in the 
solution set. A value of +1 specifies that all cosine terms 
should be sequenced before all sine terms ^ and a value of “I 
for alternating the cosine and sine terms. The value of 
CYCSEQ has been set to -1. 

12. £TYPE - fixed - The BCD value of this parameter defines the 
type of cyclic symmetry as follows: 

(1) ROT — rotational symmetry 

13. jCMAX - required - The integer value of this parameter 
specifies the maximum value of the harmonic index, and is 
used in subcase definition. There is no default for this 

parameter. The maximum value that can be specified is 
NSEGS/2 . 
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14. ^IN optional - The integer value of this parameter 
specifies the minimum value of the harmonic index. 

If KMIN (^0/ default = 0} equals KNAX^ then Parameter 
KINDEX is internally defined equal to KMIN and KMAX. User 
supplied KINDEX is ignored. 

If KMIN differs from KMAX, then KINDEX (KMIN ^KINDEX <KMAX) 
must be specified. 

15. j iNDEX optional, but see Parameter KMIN - Circumferential 
harmonic index. No default. 

16. jj l*QAD - fixed - The integer value of this parameter is the 
number of static loading conditions. The value of NLOAD is 
internally computed. 

17. ^KPRT - optional - An integer value of +1 for this 

parameter will cause the current harmonic index, KINDEX, to 

be printed at the top of the harmonic loop. The default is 

+ 1 . 

” optional - The real value of this parameter defines 
the rotational speed of the structure in revolutions per 
unit time. The default is 0.0. 

19. LFREQ and HFREQ - required unless LMODES is used - The real 
values of these parameters give the frequency range (LPREQ 
IS lower limit, and HFREQ is upper limit) of the modes to be 
used in the modal formulation. 


31 



20. LMODES - used unless set to 0 - The integer value of this 
parameter is the number of lowest modes to be used in the 
modal formulation. The default is to use all modes. 

^NMACH optional - This is the minimum Mach number at and 
above which the supersonic unsteady cascade theory is 
valid. The default is 1.01. 

22. M ^MACH - optional - This is the inaximuro Mach number at and 
below which the subsonic unsteady cascade theory is valid. 
The default is 0.80. 

23. - optional - This defines the reference streamline 
number, iref must be equal to an SLN on a STREAML2 bulk 
data card. The default value, -1, represents the blade tip 

streamline, if iref does not correspond to a valid SLN, the 
default is taken. 

24. KGGIN - optional - A positive integer value of this 
parameter indicates that the user-supplied stiffness matrix 
is to be read from tape (GINO file INPT) via the INPUTTl 
module in the rigid format. The default is -1. 

25. required - The real value of this parameter specifies 
the inflow dynamic pressure, based on the density and 
velocity on STREAML2 card for reference (PARAM IREF) 
Streamline. 

26. _BOV- required - The real value of this parameter equals the 
ratio of the semichord to the velocity on STREAML2 card for 
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reference (PARAM IREF) streamline. 


2.5 BULK DATA DECK 

No new bulk data cards have been introduced to conduct modal 
aerodynamically forced vibration analysis of turbosystems. 

Some remarks on the use of some of the bulk data cards are 
Offered as follows: 

OfJOlN. This card is used to list the corresponding GRID 
points on sides 1 and 2 of the modelled cyclic sector. 

In case of rigid hub/disk conditions, the GRID points listed 
on this card must be totally fixed. The Parameters WIM, KMIN, 
and KINDEX must be identically zero. 

In case of flexible hub/disk, the data on this card must 
reflect such boundary connections. Parameters KMM, KMIN , and 
truly active and meaningful. The displacement 
coordinate systems for any pair of corresponding GRID points must 
be axi-symmetrically compatible, i.e., the coordinate system for 
side 1 GRID point must completely coincide with that for the 
corresponding GRID point on side 2, when the side 1 coordinate 

system is rotated as a rigid body about the axis of rotation, and 
moved to side 2. 

The variables on this card represent the conditions for 
the entire blade/turbosystem as a whole. The values of these 
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variables on the reference streamline are assigned to also 
represent those for the entire blade/turbosystem. 

The reference streamline is picked by the user (PARAM IREF ) , 
and defaults to tip streamline otherwise. 

S TREAML2 . This card defines the unsteady aerodynamic data for 
a given streamline. 

Figure 2.4 illustrates some of the definitions pertinent to 
swept blade aerodynamics. 

HKAEROi. The reduced frequency on these cards is based on the 
semichord and velocity on STREAML2 card for reference 
Streamline. 

Positive inter-blade phase angle is taken when, in the 
following sketch, blade 1 LEADS the reference blade. 
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BULK DATA DECK 


Input Data Card CYJ0IN 

^ _ scription . Defines the boundary points of a seament for cyclic symmetry structural models 
Format and Example : 


5 6 7 8 


CYJ0IN 

SIDE 

C 

G1 

G2 

G3 

G4 I 

G5 ' 

66 

1 u 

abc 

CYJ0IN 

1 


7 

9 

16 

25 _J 

33 

64 

ABC 


+bc 

G7 

G8 

G9 

-etc.- 






— 

+BC 

72 










Alternate Form 


CYJ0IN 

SIDE 

C 

GIDl 

"THRU" 

GID2 





CYJ0IN 

2 

S 

6 

THRU 

32 






Field 

SIDE 

C 

Gi.GIDi 


Contents 

Side identification (Integer 1 or 2) 

Coordinate System (BCD value R,C or S or blank) 

Grid or scalar point identification numbers (Integer > 0) 


Remarks : 1 . 


symetry problems. 

parameter (CTYPE) must specify rotational or dihedral syiime try. 

one "'“st be one logical card for side 1 and 

lists must hav^ the sa« lenj?h!‘’“^^ bence both 

aSS s,-Sfrie?^'?o^he°m]dS?; *1"® ’ between segments 

referenced^ irilllfar^hr ^ “iShe^^-ca, . 
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BULK DATA DECK 


Input Data Card A^ Aerodynamic Physical Data 

De scription : Gives basic aerodynamic parameters. 


Format and Example*. 


AER0 

AER0 

ACS ID 
3 

VELOCITY j 
1 . 3+4 

REFC 

100. 

0 

RH0REF 

1.-5 

6 

SYMXZ 

7 

SYMXY 

1 

8 

9 

10 


Field 

ACS ID 

VEL0CITY 

REFC 

RH0REF 

SVMX2 

SYMXY 

Remarks : 1 

2 . 

2. 


Contents 

Aerodynamic coordinate system Identification {Integer > 0). See Remark 2 
Velocity (Real). 

Reference length (for reduced frequency) (Real). 

Reference density (Real). 

-no%‘'^n«^:sy'S).“®™ (-1 for sym. =0 for no sym 

effects^( Integer) !*sMe'code"as*SYH« Sround 

Aorodynamlc problems. Only one AER0 card 1s 

' rectangular coordinate system. Flow Is In the positive x 

I, _ gGFc/7. 


k-_ 


V 
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BULK DATA DECK 


Input Data Card HKAER^I Mach Number - Frequency Table 

Description : Provides a table of Kach numbers or interblade phase angles (m) 

and reduced frequencies (k) for aerodynamic matrix calculation. 


Format and Example : 


^ , 5 6 7 8 9 10 


HKAER01 

mi 


— 

m 3 

?4 

mg 

me 

m7 

me 

ABC 

HKAER01 

.1 

.7 







+ABC 

+ec 

^1 


^3 

■ 

k4 

^5 

k6 

^7 


1 

+BC 

.3 

.6 

1.0 





t 

1 


Field 


Contents 


m. 

1 


List of Mach numbers (Real; 1 ^ i 8). 

List of reduced frequencies (Real > 0.0, 1 j i 8), 


Remarks: 1. Blank fields end the list, and thus cannot be used for 0.0. 

2. All combinations of (m,k) will be used.' 

3. The continuation card is required. 

sich%?'0.000K* with a very small nuaber 

5. 

^ ^re input for wing flutter and interblade phase angles for 
oiaae ti utter and response. 
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BULK DATA DECK 


Input Data Card hKAER^3 Mach Number - Frequency Table 


Description: 


Provides a list of Mach numbers or interbladc phase angles (ni) and 
reduced frequencies (k) for aerodynamic matrix calculation. 


Format and Example : 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

MKAER02 

"1 


m, 

^2 






MKAER02 

.10 

.30 

.10 

.60 

.70 

.30 

.70 

1.0 



Field Contents 

m^. List of Mach numbers (Real > 0.0). 

List of reduced frequencies (Real > 0.0). 


Remarks : 1. This card will cause the aerodynamic matrices to be computed for a set of parameter 

pairs. 

2. Several MKAER02 cards may be in the deck. 

3. Imbedded blank pairs are skipped. 

4. Mach numbers are input for wing flutter and interblade phase angle for blade 
flutter and response. 


38 















BULK DATA DECK 


Input Data Card smMkL . B1 ade Straaml (na Data 


Oescri pti on : 


Defines grid points on 
edge to blede trailing 


the blade streamline 
edge . 


from blade leading 


Format and Example : 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

STREAMLl 

SLN 

Gl 

62 

G3 

G4 

G5 

G6 

G7 

+ABC 

STREAHLl 

3 

2 

4 

6 

6 

10 








‘-I 1—1 rJ 


+ABC 

G8 

G9 

-etc- 







+ABC 

.. 








i 


Alternate Form: 


STREAMLl 

SLN 

GlOl 

"THRU" 

GI02 






STREAMLl 

5 

6 

THRU 

12 







Field 

SLN 

Gi, GIDi 

Rema rk s : 

1 . 


Contents 

Streamline number (integer > 0). 

Grid point identification numbers (integer > 0). 


This card is requ 
and response pfo 


blade steady aeroelastic, 


blade flutter, 

O^rou^on ' sections closest to and farthest from the axis 

4, All grid points should be unique. 

5. All grid points referenced by GIOl through GI02 must exist. 

nJSJs'IusHe V.lttTrL'X lu.lT 9-*^ PO<"ts. The 

trailing edg: I^^JJe^TrecI ^otuio':?’ 


39 






BULK DATA DECK 


Input Data Card STREAML2 Blade Streamline Data 

Description: Defines aerodynamic data for a blade streamline. 

Format and Example : 


STREAML2 

SLN 

NSTNS 

STAGGER 

CHORD 

“radius/ 

0CBD2B 

BSPACE 

MACH 

DEN 

+abc 

STREAML2 

2 

3 

23.5 

1.85 

6.07 

.806 

.934 

.066 



+abc 

VEL 

FLOWA/ 

SWEEP 








+ABC 

1014.2 

55.12 









Contents 

Streamline number (Integer >0) 

Humber of computing stations on the blade streamline. 

(3 < NSTNS <10, Integer) 

Blade stagger angle (-90.0 <stagger <90.0, degrees) 

Blade chord (real >0.0) 

Radius of streamli^ie for dyiuiaic .analysis without sweep effects 
(real >0.0) or 


Field 

SLH 

NSTNS 

STAGGER 

CHORD 

RADIUS/DCBDZB 


BSPACE 

MACH 

DEH 

VEL 

FLOWA/ SWEEP 


3C/3I for dynamic analysis with sweep effects. C is the swept 
chord and 7 is the (local) spanwise reference direction (real) 
Blade spacing (real >0.0) 

Relative flow mach number at blade leading edge (real >0.0) 

Gas density at blade leading edge (real >0.0) 

Relative flow velocity at blade leading edge (real >0.0) 
Relative flow angle at blade leading edge for dynamic^ analysis 
without sweep effects (-90.0 <FL0WA <90.0 degrees) w 
Blade sweep angle fordynamic analysis with sweep effects 
(“90,0 <SWEEP <90.0 degrees) 
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At least three (3) and no more than fifty (50) STREAf-1L2 cards are required 
for a blade dynamic analysis. 


The streamline number. SLN. must be the same as Its corresponding SLN on 
a STREAMLl card. There must be a STREAMLl card for each STREAH2 card. 

It IS not required that all streamlines be used to define the aerodynamic 
matrices used in blade dynamic analysis. 

For dyn«.ic analysis with sweep effects, the use of the NASTRAR card is 
required as follows: 

NASTRAN SYSTEM (93) = 1 

Refer to Section 2.1 of the User's Manual and Section 6.3.1 of the Progra-e 
Manual for description and placement in the Executive Control Deck. 

Dynamic analysis refers to both flutter and response 
analyses. 



NASTRAN Model 
and 

Operating Point 


RUN 1 I 

Differential 
Stiffness 
Analysis 
DISP RF 4 


Total Stiffness 
KTOTAL 


-H 


RUN 2/3 I 

Aeroelastic 
Stability 
Analysis 
AERO RF 9 


-/f\ 


± 


I 


AERO and 
STREAMLJ 

^ I 

(If desired) 


Natural Frequencies 
Mode Shape Plots 
Flutter Summary 
V~g V-f Plots 


RUN 3/2 4 

Applied 

Oscillatory 

Airloads 

Generation 

AIRLOADS 


Cose Control 
Bulk Data 


RIes 


RUN 4 

Modal 

Forced Vibration 
Analysis 

DISP RF 8 with 

MFVAAET After Pockoje 




Natural Frequencies 
Mode Shape Plots 
Dynamic Response 
Stresses to compare with 
Wind Tunnel Dota 


JFigure ;'2.1 Overall Flowchart of Steps to Conduct 
Modal Forced vibration Analysis of 
Aerodynamically Excited Turbosystems 












NOTES 


Planes 
to XjZ^. 




only be parallel 


2. axis is parallel to 

3. axis is parallel to 

4. Uniform Inflow is along +X^ 

coord, system 





Tunnel 

coord. 

system 



Uniform 

Inflow 


Figure 2. 3 Turboprop Axis Inclination Angle and 
Tunnel Coordinate System Orientation 
in Uniform Inflow Case 
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MODAL FORCED VIBRATION ANALYSIS OF 
AERODYNAMICALLY EXCITED TURBOSYSTEMS 

3.1 INTRODUCTION 

The NASTRAN programming effort for conducting modal forced 
vibration analysis of aerodynamically excited turbosystems 
involved extensive modifications of, and additions to, the source 
code in the CRAY version of the April 1984 release of NASTRAN. in 
addition, the effort also included the adaptation of an existing 
DMAP ALTER package to this version of NASTRAN, as well as the 
development of a new DMAP ALTER package. 

The above changes fall into two groups. The first group 
involves changes made in order to incorporate the UNIVAC version 
of Bladed Shrouded Disks Computer Program enhancements into CRAY 
NASTRAN. (Complete details of these enhancements are contained in 
References 3 through 7,), The second group involves changes made 
in order to incorporate the Updated Bladed Disks Program 

enhancements developed under the present contract into CRAY 
NASTRAN. 

Details of all of the above changes are discussed in the 
following subsections. 


3-2 ONIVAC TO CRAY CONVERSION OF BLADED DISKS PROGRAM 

The incorporation of the ONIVAC version of the Bladed Disks 
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Program into CRAY involved extensive changes to the CRAY version 

of April 1984 NASTRAN. These changes are due to the following 
reasons : 


!• Addition of two new bulk data cards: 


STREAMLl — Defines grid points on a blade 

streamline from the leading edge 
to the trailing edge 

STREAML2 — Defines aerodynamic data for a 
blade streamline 


2 . 


Addition of the 15 new bulk data parameters (PARAMs) : 


APRESS 

KGGIN 

SIGN 


r ATEMP , FXCOOR , PYCOOR , FZCOOR , IREP 
, KTOUT , MAXMACH, MINMACH, MTYPE ^ PGEOM 
r STREAML and ZORIGN 


f 

9 


3. Addition of four new functional modules: 


4. 


ALG 

APDB 

PVRSTRl 

PVRSTR2 


Aerodynamic load generator 
Aerodynamic pool distributor for blades 
Forced vibration response analysis of 
rotating cyclic structures - Phase 1 
Forced vibration response analysis of 
rotating cyclic structures - Phase 2 


Addition of two new rigid formats: 


DISP RF 16 
AERO RF 9 


Static aerothermoelastic analysis 
with differential stiffness 
Compressor blade cyclic modal 
flutter analysis 


5. 


aerodynamic effects via the 
®3rd word of COMMON 

/SYSTEM/ on the NASTRAN card) 


6 . 


7. 


Adaptation of a DMAP ALTER package for DISJ 
developed for the forced vibration analysii 
cyclic structures using the direct approacl 


RF 8 

of rotating 


Various other miscellaneous changes 


It is noted here that the addition of the two new rigid formats 
mentioned above involved not only modifications to the source 
code, but also expansion of NASTRAN rigid format data base by the 
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addition of the DMAP sequences and the restart tables for these 
two new rigid formats. 


Modifications to the Source Cod e 

A total of 26 existing subprograms were modified in the UNIVAC 
to CRAY conversion process. The modified subprograms, the reasons 
for their modification, and the extent of the changes are 
summarized in Table 3.1. 

The actual changes to the code have been made such that they 
are easily identified. Thus, lines of code that have been 
logically deleted have not been physically deleted, but have been 
commented out by using the string CBELD in columns 1 through 5. 
Lines of code that have been logically replaced have been 
commented out by using the string CBELR in columns 1 through 5, 
with the replacement line (or lines) immediately following the 
replaced code. New lines added to the code have been bracketed by 
the comment line CBELNB at the beginning, and the comment line 
CBELNE at the end of the group of lines added to the code. 


Additions to the Source Code 

A total of 82 new UNIVAC subprograms were added to CRAY NASTRAN 

in the conversion process. The added subprograms are identified 
in Table 3.2. 
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^•2-3 Adaptation of the DMAP ALTER Package 


As part of the UNIVAC version of the Bladed Disks Program, a 
DMAP ALTER package for DISP RF 8 had been developed for the forced 
vibration analysis of rotating cyclic structures using the direct 
approach. This ALTER package was adapted to the CRAY version of 
April 1984 NASTRAN by changing all the ALTER numbers to correspond 
to the latest version of DISP RP 8. In addition, the ALTER package 
was streamlined by making the following changes: 


1 . 


All C^PNT instructions were removed as the already 

format causes 

the CHKPNTing of all required data blocks. 


2. All SAVE instructions of the form SAVE FLAG were 

features of the form 

/i>,w,FLAG/ in the immediately preceding DMAP 
instructions. 


constant parameters of the form 
/C,N,M/ were replaced by /M/. 

4. All BCD constant parameters of the form /C#N#B/ 

were replaced by /*bV. / / / / 

5. All variable parameters of the form /V,N,VAR/ 

immediately SAVEd were replaced by 


The DMAP ALTER package resulting from the above changes has 
been named DPVARCS (Direct Forced Vibration ^alysis of R^otating 
Cyclic Structures) , and is listed in Table 3.3. 


^•3 NEW UPDATES IN BLADED DISKS PRCXSRAM ON CRAY 

These new enhancements also involved extensive changes to the 
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CRAY version of April 1984 NASTRAN. These chan9es are primarily 
due to the following reasons: 


1. Development of a new DMAP ALTER package for DISP 

RF 8 for the forced vibration analysis of aerodynamically 
excited turbosystems using the modal approach, including 
the dynamic effects on these rotating cyclic structures. 

2. Ability to rearrange, rather than interpolate, the 
reactionary aerodynamic matrix data in the FRRD2 module 
in analyses using the newly-developed DMAP ALTER 
package. 

3. Expansion of the CURV and OFP modules to compute 
and process complex element stresses in material 
coordinate systems for the CQUADl/2 and CTRIAl/2 
elements . 

4. Various other miscellaneous changes. 


3.3.1 Modifications to the Source Code 


A total of 26 subprograms were modified in order to incorporate 
the new Updates to the Bladed Disks Program into the CRAY version 
of April 1984 NASTRAN. The modified subprograms, the reasons for 

their modification, and the extent of the changes are summarized 
in Table 3.4. 


The actual changes to the code have been made such that they 
are easily identified. Thus, lines of code that have been 
logically deleted have not been physically deleted, but have been 
commented out by using the string CAERD in columns 1 through 5. 
Lines of code that have been logically replaced have been 
commented out by using the string CAERR in columns 1 through 5, 
with the replacement line (or lines) immediately following the 
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replaced code. New lines added to the code have been bracketed by 
the comment line CAERNB at the beginning, and the comment line 
CAERNE at the end of the group of lines added to the code. 

3.3.2 Additions to the Source Code 

Two new subroutines (OPPCCl and OFPCC2) were added to the OPP 
module. These subroutines, in conjunction with the other changes 
made to the OPP module, are designed to process the complex 
element stresses in material coordinate systems for the CQUADl/2 
and CTRlAl/2 elements, in both SORTl and S0RT2 output formats. 

3-3.3 Development of a New DMAP ALTER Package 

As part of the new updates to the Bladed Disks Program, a new 
DMAP ALTER package for DISP RP 8 was developed for the forced 
vibration analysis of aerodynamically excited turbosystems using 
the modal approach. The DMAP ALTER package for DISP RP 8 
developed for the direct approach forced vibration analysis of 
rotating cyclic structures (Table 3.3) was used as the starting 
point for this effort. This ALTER package was then modified and 
considerably expanded by making the following major changes: 


1. The READ module was added in order to compute the 
eigenvalues and eigenvectors of rotating cyclic 
structures. Provision to obtain plots of modal 
information was made by the inclusion of the CYCT2, 
SDRl, SDR2 and PLOT modules just after the READ module. 
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2. The GKAM module was added to generate the modal 
matrices . 

3. The APDB, AMG and AMP modules were inserted in 

order to generate the reactionary aerodynamic matrix data 
tor subsequent use in the solution phase in the FRRD2 
module. 

4. Various other changes to the DMAP were made in 
order to define various control flags and to 
permit the execution of the solution procedure in 
an elegant manner. 


The DMAP ALTER package resulting from the above changes has 
been named MFVAAET (Modal Forced Vibration Analysis of 
Aerodynamically ^xcited ^rbosystems) / and is listed in Table 3.5. 


A complete listing of the DISP RF 8 Rigid Format as modified by 
the MFVAAET Alter package is presented in Table 3.6. 
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TABLE 3.1 


CRAY NASTRAN Subprograms Modified to Incorporate 
UNIVAC Version of Bladed Disks Program 


Subroutine Subprograms (Total: 18) 


Subroutine 


Reason for modification Extent of changes 


AMG 


AMPC 

CYCT2 

PA2 

PRD2B 

PRD2C 

IPP 


IPPPAR 


IPS3P 

NASTOl 

TTLPGE 


XCSA 

XLNKDD 


XMPLDD 

XRGDPM 

XSEN06 

XSEM07 


XSEM09 


Extend analysis to 

compressor blade and 

swept turboprop blade 

methods 

Same as AMG 

Miscellaneous 

Miscellaneous 

Miscellaneous 

Miscellaneous 

To incorporate two new bulk 
data cards (STREAMLl and 
STREAML2) 

To incorporate the two new 
rigid formats (DISP RP 16 
and AERO RP 9) 

Same as IPP 

Miscellaneous 

To change the title page 

in order to reflect the 

current changes 

Same as IPPPAR 

To incorporate four new 

modules (ALG, APDB, PVRSTRl 

and PVRSTR2) 

Same as XLNKDD 
Same as IPPPAR 
To incorporate the new 
ALG module into Link 6 
To incorporate the new 
PVRSTRl and PVRSTR2 
modules into Link 7 
To incorporate the new 
APDB module into Link 9 


Moderate 


Moderate 

Moderate 

Moderate 

Minor 

Minor 

Moderate 


Moderate 


Moderate 

Minor 

Moderate 


Moderate 

Moderate 


Moderate 

Moderate 

Moderate 

Moderate 


Moderate 


(continued) 
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TABLE 3.1 (continued) 


Block Data Subprograms (Total: 8) 


Block data 


GPTABD 

IFXIBD 


IFX2BD 


IFX3BD 

IFX4BD 

IFX5BD 

IFX6BD 

SEMDBD 


Reason for modification Extent of changes 


Miscellaneous Minor 

To incorporate the new Moderate 

STREAMLl and STREAML2 

bulk data cards and 

15 new bulk data 

parameters (PARAMs) 

To Incorporate the new Minor 

STREAMLl and STREAML2 
bulk data cards 

Same as IFX2BD Minor 

Same as IFX2BD Minor 

Same as IFX2BD Minor 

Same as IFX2BD Minor 

To permit the addition Minor 

of the new ALG, APDB, 


FVRSTRl and FVRSTR2 modules 


Total Subprograms: 26 
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TABLE 3.2 


New Subprograms Added to CRAY NASTRAN to Incorporate 
ONIVAC Version of Bladed Disks Program 


Subroutine Subprograms (Total: 74) 


ALAMDA 

AKAPM 

AKAPPA 

AKP2 

AMGBl 

AMGBIA 

AMGBIB 

AMGBIC 

AMGBID 

AMGBIS 

AMGB2 

AHGB2A 

AMGTl 

AMGTIA 

AMGTIB 

AMGTIC 

AMGTID 

AMGTIS 

AMGTIT 

AMGT2 

AMGT2A 

APDB — Driver for APDB module 

APDBl 

APDB2 

APDB2A 

ASYCON 

DLKAPM 

DRKAPM 

FVRSTl — Driver for FVRSTRl module 

FVRST2 — Driver for FVRSTR2 module 

FVRSIA 

FVRSIB 

FVRSIC 

FVRSID 

FVRSIE 

FVRS2A 

GAUSS 

INTERT 

RETNL06 (See Note below) 

SOBA 

SUBBB 

SUBC 


(continued) 
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TABLE 3.2 (continued) 


SUED 

UDGl 

UD03AN 

UD03AP 

UD03AR 

UD03PB 

UD03PO 

UD03PR 

UD0300 — Driver for ALG module 

UD0301 thru UD0319 

UD0325 

UD0326 

UD0329 

UD0330 


Function Subprograms (Total: 8) 


UDG2 thru DDG9 


Total Subprograms: 82 


Note: Subroutine RETNL06 is required in order to provide 

dunuay returns for certain subroutines referenced in 
Link 6. 
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ALTER 3 % 

FILE UXVF=APPEND/PDT-APPEND/PD=APPEND * 

* PERFORM INITIAL ERROR CHECKS ON NSEGS AND KMAX. 

COND ERRORCl, NSEGS % IF USER HAS NOT SPECIFIED NSEGS. 

COND ERRORCltKMAX S IF USER HAS NOT SPECIFIED KMAX, 

PARAM //^^EQ*/CYCIOERR /V»Y,CYCI0=0 /O « 

COND ERRORCl ,CYCIOERR $ IF USER HAS NOT SPECIFIED CYCIO. 

PARAM //«DIV*/NSEG2 /V,Y, NSEGS /2 * NSEG2 = NSEGS/2 

PARAM //■JtSUB^/KMAXERR /NSEG2 /V,Y,KMAX * 

COND ERRORCl, KMAXERR % IF KMAX .GT. NSEGS/2 

« SET DEFAULTS FOR PARAMETERS. 

PARAM //*NOP*/V, Y,NOKPRT=+i /V, Y, LGKAD=-1 « 

% CALCULATE OMEGA, 2*0MEGA AND OMEGA**2 FROM RPS. SET DEFAULT RPS. 
PARAMR //*MPY*/OMEGA /V,Y,RPS«0.0 /6. 283185 $ 

PARAMR //*MPY*/0MEGA2 /2.0 /OMEGA » 

PARAMR //*MPY*/OMEGASQR /OMEGA /OMEGA $ 

^ GENERATE NORPS FLAG IF RPS IS ZERO. 

PARAMR //*EQ*//V, Y,RPS /O.O ////NORPS % 

% MAKE SURE COUPLED MASSES HAVE NOT BEEN REQUESTED. 

PARAM //*NOT*/NOLUMP /V,Y,COUPMASS=-l * 

COND ERR0RC2, NOLUMP $ 

ALTER 21,21 * ADD. SLT TO OUTPUT FOR TRLG. 

GP3 6E0M3,EQEXIN,6E0M2 / SLT,GPTT / NOGRAV % 

ALTER 23 ^ 

i SINCE MULTIPLE CONSTRAINTS ARE NOT ALLOWED EXECUTE GP4 NOW SO THAT 
% MORE ERROR CHECKS CAN BE MADE BEFORE ELEMENT GENERATION. 
t ADD YS NEEDED FOR PSF RECOVERY IN SSG2. 

PARAM //*MPY*/NSKIP /O/O $ 

GP4 CASECC , 6EOM4 , EQEX IN, GPDT , BGPDT , CSTM , /RG , YS , USET , ASET/LUSET / 

S, N, MPCFl/S, N, MPCF2/S, N, SINGLE/S, N, OMIT/S, N, REACT/S, N, NSKIP/ 

S, N, REPEAT/S, N, NOSET/S, N, NOL/S, N, NOA/C, Y, ASETOUT/S, Y, AUTOSPC t 
PURGE GM,GMD/MPCFl/GO,GOD/OMIT/KFS,PSF,QPC/SINGLE * 

i SUPORT BULK DATA IS NOT ALLOWED, 

PARAM //*NOT*/REACDATA /REACT $ 

COND ERR0RC3,REACDATA % 

* EXECUTE DPD NOW SO CHECKS CAN BE MADE. ADD TRL TO OUTPUT DATA BLOCKS. 

DPD DYNAMICS, GPL, SIL, USET / GPLD, SILD, USETD, TFPOOL, DLT, PSDL, FRL, , 

TRL, ,EQDYN / LUSET/S, N, LUSETD/NOTFL/S, N, NODLT/ 

S, N, NOPSDL/S, N, NOFRL/NONLFT/S, N, NOTRL/NOEED/ / 

S,N,NOUE $ 

MUST HAVE EITHER FREQ OR TSTEP BULK DATA. 

PARAM //*AND*/FTERR /NOFRL /NOTRL $ 

COND ERR0RC5,FTERR * NO FREQ OR TSTEP BULK DATA. 

* ONLY FREQUENCY OR TSTEP IS ALLOWED IN THE CASE CONTROL 

PARAML CASECC //*DTI*/1/14//FREQSET % 

PARAML CASECC / /*DTI */ 1 /38//T IMESET $ 

PARAM //*MPY*/ FREQTIME /FREQSET /TIMESET i 

PARAM //*NOT*/FTERRl /FREQTIME $ 

PARAM //*LE*/NOFREQ /FREQSET /O % 

PARAM //*LE*/NOTIME /TIMESET /O $ 

COND ERR0RC6, FTERRl i BOTH FREQ AND TSTEP IN CASE CONTROL DECK. 

$ EPOINT BULK DATA NOT ALLOWED 

PARAM //«NOT*/EXTRAPTS /NOUE t> 

COND ERR0RC4, EXTRAPTS $ 

GENERATE DATA FOR CYCT2 MODULE. 

GPCYC GE0M4,E0DVN, USETD /CYCDD /CTYPE^^ROT /3,N,NOGO 
COND ERRORS l.NGGO 
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M2GG AND BASEGX. 

/ FRLX,B1GG,M1G6, 


PARAML 

COND 

EQUIV 

LABEL 


ALTER 33 $ 

$ PRE-PURGE DATA BLOCKS THAT WILL NOT BE GENERATED 
PARAM //*OR*/NOBMl /NOMGG /NORPS $ 

PURGE B1GG»M1GG /NOBMi $ 

PURGE M2GG,M2BASEXG /NOMGG % 

ALTER 36 $ 

% GENERATE' DATA BLOCKS FRLX, BIGG, MIGG, 

% GENERATE PARAMETERS FKMAX AND NOBASE X. 

FVRSTRl CASECC, BGPDT, CSTM, DIT, FRL, MGG. 

M^GG, BASEXG, PDZERO, , /NOMGG/V, Y, CYCIO/V, Y, NSEGS/ 

V, Y,KMAX/S,N,FKMAX/V, Y,BXTID=-1/V,Y,BXPTID=-1/ 

V , Y , BYT I D=- 1 / V , Y , B YPT I D=- 1 / V , Y , BZ T I D=- 1 / 

V, Y, BZPTID--1/S, N, NOBASEX/NOFREQ/OMEGA $ 

FRLX //*PRESENCE*////NOFRLX t 
LBLFRLX,NOFRLX $ 

FRLX , FRL $ 

LBLFRLX % 

ALTER 43 % 

PARAM //*ADD*/NOBGG /NOBMI /O $ 

ALTER 53 % 

% REDEFINE BGG AND KGG. 

COND LBLllA.NOBMl $ 

PARAMR //*COMPLEX*// OMEGA2 /O.O/ CMPLXl * 

PARAMR //*SUB*/ MOMEGASQ / 0.0 / OMEGASQR * 

PARAMR //^COMPLEX*// MOMEGASQ / 0.0 / CMPLX2 $ 

BGG, BIGG / BGGl / (1.0, 0.0) / CMPLXl $ 

BGG 1, BGG $ 

KGG, MIGG / KGGl 
KGG1,KGG $ 

LBLllA 

ALTER 54,56 % GP4 HAS BEEN MOYED-UP. 

ALTER 88,88 $ DPD HAS BEEN MOVED-UP* 

LGKAD1,LGKAD $ BRANCH IN NOT FREQRESP. 


RESET NOBGG. 


ADD 

EQUIV 

ADD 

EQUIV 

LABEL 


PARAM 

COND 


/ <1.0,0. 0) / CMPLX2 % 


ALTER 114 SEE ALTER 113 COMMENT. 


JUMP 

LABEL 

EQUIV 

LABEL 


LGKAD2 $ 

LGKADl $ 

KAA^KDD/^%EKA''$^^^ B2DD/N0A/K2PP, K2DD/N0A/MAA, MDD/MDEMA/ 
LGKAD2 $ 


BkIT AND N4 TO GKAD. 


ALTER 

COND 


I lctr-rT\ ’ n-T i OISMU. 

GOn^K-nn’ K2PP, M2PP, B2PP/KDD, BDD, MDD, GMD, 

GOD,K^DD,M.iDD,B^DD/C,Y,GKAD=TRANRESP/*DISP*/*DIRECT*/ 

C, Y, G-0. 0/C, Y , W3=0. 0/C, Y, W4=0. 0/N0K2PP/N0M2PP / 

N0B2PP/MPCF1/SINGLE/0MIT/N0UE/N0K4Gg7 

N0BGG/KDEK2/-1 % 

117 * SEE ALTER 113 COMMENT, 

LGKAD3,LGKAD % BRANCH IF NOT FREQRESP. 


11® * SEE ALTER 113 COMMENT. 


JUMP 

LABEL 

EQUIV 

LABEL 


LGKAD4 * 

LGKAD3 

B2DD,BDD/N06PDT/M2DD,MDD/N0SIMP/K2DD,KDD/KDEK2 $ 

LGKAD4 $ ■ 

ALTER 119,123 % 

* NEW SOLUTION LOGIC 


LBLTRLl , NOTIME * 


27 S^^^'^SES FOR TIME-DEPENDENT LOADS. 


PARAM 

PARAM 

LABEL 

CASE 

PARAM 

TRLG 


//*MPY-«/REPEATT /I /-i * 

//#ADD*/APPFLG /I /O % INITIALIZE FOR SDRl. 

TRLGLOOP t 

NOLoopi 

//*MPY*/NC0L /O /I $ 

CASEYY, LISETD, DLT, SLT, BGPDT, S I L , CSTM, TRL, DIT, GMD, GOD, , EST, MGG/ 


59 



SDRl 

SDRl 

PARAM 

COND 

REPT 

JUHP 

LABEL 

fvrstr: 


EQUIV 

JUMP 

LABEL 


/APPFLG/*DYNAMICS* $ 
* APPFLG=APPFLG+1 . 


5 , PDTl , PDl , , TOL/ NOSET/NCOL » 

” ’ ” ” ^ »PDT, /APPFLG/*DYNAMICS* 

TRL,PD1 1 , / ,PD , - 

//#ADD*/APPFLG /APPFLG /I 
TRLGDONE, REPEAT! $ 

TRLGLOOP, 100 $ 

ERR0R3 % 

TRLGDONE $ 

^ FQLZ, REORDER!, RE0RDER2, , , 

V, Y, NSEGS/V, Y, CYCIO/S, Y, LMAX<=-1/FKMAX/ 

S, N, FLMAX/S, N, NTSTEPS/S, N, NOROl /S, N, N0R02 $ 
FRLZfFRL // FOLZ,FOL « 

LBLFRL2 $ 

LBLTRLI t 


.^GENeR«TE^PREQUENCV-OEPE^^^^^^ SEUECTED IN CC. 


COND 
MPYAD 
EQUIV 
LABEL 
it FORM NEW LOADS. 
COND 
MPYAD 
ADD 
EQUIV 
COND 


ractrvv I.oi— r, E , r Y WHb selected IN Cr 

*DIRECT*/FREQY/*FREQi°*°’^°^’^^^’ ^ PSF, PDF, FOL, PHFDUM / 

PPFrPD2ERS?''^^PFX /f 2 ™>-X GENERATED. 

PPFX,PPF % 

LBLFRLXl % 


LBLFRL1,N0BASEX $ 

M2GG,BASEXG, / M2BASEXG /O t 

ppf;?ppPT° ' Xd. 0,0.0) /(-1. 0,0.0) t 

LBLBASEl, NOSET % 

/ »P0DUM1,PSF1,PDF1 $ 
PSF1,PSF // PDF1,PDF $ ^ j. ^ 

LBLBASEl % 

LBLFRLl ♦ 

PPF, PDF/NOSET % 

* LOADS ARE FREQUENCY-DEPENDENT 

IF CYCIO=.l. 

f'HKHnL PDF //*TRAILER*/1 /PDFCOLS * 

number of loads FOR CYCI0=-1. 

mC?v PDpfpxF^S = -^^XEKNAX 

COND LBLPD0NE,CYCI0 $ 

* number of loads for CYC 1 0=1. 


SSG2 

EQUIV 

LABEL 

LABEL 

EQUIV 


PARAM 

CYCTl 

COND 

JUMP 

LABEL 


PDF‘'/''pxF'-GrvrF?“^F?bL^'^’^’^®^®® * = NF/NSEGS 

V Y I'MaJ-’? T ^,L ^ /»fore»/v,y,nsegs=-i / 

V,Y,kMAX — 1 / NLOAD /S,N,NOGO $ 

ERR0RC1,N0G0 $ 

LBLPDONE % 

LBLFRL2 it 


* LOADS ARE TIME-DEPENDENT 

PARAM //*N0T*/N0TCYCI0 /V,Y,CYCI0 * 

% BRANCH DEPENDING ON VALUE OF CYCIO 

OOKir\ I m.t 


COND 
$ CYCIO 
EQUIV 
COND 
MPYAD 
LABEL 
CYCTl 

COND 
EQUIV 
COND 
MPYAD 
LABEL 
EQUIV 
JUMP 
LABEL 
$ CYCIO = +1 


LBLTRL2,N0TCYCI0 $ 

-1 

PD, PDTRZl /NOROl * 

LBLROIA, NOROl $ 

PD,RE0RDER1, / PDTRZl / 0 % 

LBLROIA % 

^ F>^TRZ1 ,GCYCF2 /CTYPE/*FORE*/NTSTEPS/ 
V, Y, LMAX/FKMAX/S, N, NOGO % 

ERRORC1,NOGO $ 

PXTRZl ,PXFZ1/N0R02 $ 

LBLRQ2A,N0R02 % 

PXTRZ 1 , RE0RDER2, / PXFZl 
LBLR02A $ 

PXFZl, PXFl $ 

LBLTRL3 $ 

LBLTRL2 % 


/O $> 
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MPYAD ■ POiREQRDtRl,' / PDVk‘Z2 7 0 % 

CYCTl PDTRZ2 /PXTRZ2, GCYCF3 /CTYPE/*FORE*/NTSTEPS/V, Y» LMAX/ 
V»Y,NSEGS/S,N,NOGO t 
CQND ERRORCliNOGO % 

EQUIV PXTRZ2,PXTR2/N0R02 « 

COND LBLR02B,N0R02 % 

MPYAD PXTRZ2,RE0RDER2, / PXTR2 /O * 

LABEL LBLR02B $ 

CYCTl PXTR2 / PXFZ2,GCYCF4 / CTYPE/*FORE*/V, Y, NSEGS/V, Y, KMAX/ 

FLMAX/S,N,NOGO $ 

COND ERRORCUNOGO * 

EQUIV PXFZ2,PXF1 $ 

LABEL LBLTRL3 ♦ 

% TIME-DEPENDENT LOADS ARE REAL. MAKE LOADS COMPLEX TO CORRESPOND 
$ TO FREQUENCY DEPENDENT LOADS. ALSO SDR2 EXPECTS LOADS TO BE COMPLEX 

* IN FREQRESP TYPE PROBLEMS. 

COPY PXFl / PXF2 $ CONVERT REAL PXFl TO COMPLEX PXF. 

ADD PXF1,PXF2 / PXF / (0.5, 1.0) / (0.5, -1.0) t 

% DEFINE NLOAD FOR CYCT2. 

PARAM //*ADD*/NLOAD /FLMAX /O * NLOAD = FLMAX 
LABEL LBLPDONE * 

PARAM //*ADD*/KINDEX /V,Y,KMIN=0 /O $ INTITIALIZE KINDEX, 

$ 

* INITIALIZE UXVF IF KMIN IS NOT ZERO. 
i 

PARAM //*ADD*/KMINL /V,Y,KMIN /-I % 

COND NOKMINL,KMINL % 

PARAM • //*ADD*/KMINV /O /O % 

LABEL KM INLOOP % 

CYCT2 CYCDD, , ,PXF, , /,,PKFZ,, /*FORE*/V, Y, NSEGS/ 
KMINV/CYCSEQ/NLOAD/S,N,NOGO * 

COND ERR0RC1,N0G0 * 

ADD PKFZ, / UKVFZ / (0.0, 0.0) $ 

PRTPARM //0/*KINDEX* % 

CYCT2 CYCDD, ,, UKVFZ, , /,,UXVF,, /*BACK*/V, Y, NSEGS/ 
KMINV/CYCSEQ/NLOAD/S,N,NOGO * 

PRTPARM //0/*KINDEX* * 

COND ERR0RC1,N0G0 * 

PARAM //*ADD#/KMINV /KMINV /i % 

REPT KM INLOOP, KM I NL « 

LABEL NOKMINL % 

LABEL TOPCYC % LOOP ON KINDEX 
COND NOKPRT,NOKPRT $ 

PRTPARM //O /*K INDEX* * 

LABEL NOKPRT $ 

CYCT2 CYCDD, KDD, MDD, , , /KKKF, MKKF, , , /*FORE*/V, Y, NSEGS / 

KINDEX/CYCSEQ=-1/NL0AD/S,N, NOGO % 

COND ERR0RC1,N0G0 % 

CYCT2 CYCDD, BDD, , PXF, , /BKKF,,PKF,, /*FORE*/V, Y, NSEGS/ 
KINDEX/CYCSEQ/NLOAD/S,N,NOGO « 

COND ERR0RC1,N0G0 « 

$ SOLUTION 

FRRD2 KKKF, BKKF, MKKF, ,PKF,FOL / UKVF /O. 0/0. 0/-1 . O % 

CYCT2 CYCDD, ,, UKVF, , /,,UXVF,, /*BACK*/V, Y, NSEGS/KINDEX/ 
CYCSEQ/NLOAD/S,N,NOGO % 

COND ERR0RC1,N0G0 it 

PARAM //*ADD*/KINDEX/KINDEX/1 % KINDEX = KINDEX + 1 
PARAM //*SUB*/D0NE / V,Y,KMAX / KINDEX $ 

COND LCYC2,D0NE IF KINDEX .GT. KMAX THEN EXIT 

REPT TOPCYC, 100 % 

JUMP ERR0R3 $ 

LABEL LCYC2 « 

EQUIV UXVF,UDVF / CYCIO % 

COND LCYC3, CYCIO $ IF CYCIO . GE. 0 THEN TRANSFORM TO PHYSICAL. 

CYCTl UXVF / UDVF,GCYCB1 / CTYPE/*DACK* /V , Y , NSEGS/ V , Y , KMAX/ 
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LABEL 

COND 

EQUIV 

COND 

CYCTl 

LABEL 


NLOAD^ 

LCYC3 % 

LBLTRL4, NOTIME $ 

PXF,PDF2 / CYCIO $ 

CYCIO .GE. 0 THEN TRANSFORM TO PHYSICAL 

NL^D ^ CTYPE/*BACK*/V,Y,NSEGS/V,Y,KMAX/ 

LCYC4 $ 


sDRi PSP ppoM pxp. 


isG2 .'.TV pp?^r; /r/«?Nfi«;csn 

lauiV PP™r;r;sPZ;PSp’r'^^' ' .POSUM.PSFZ.PLDUM * 

LABEL LBLTRL4 ♦ 

^QUDVCI, /.FREQRESP./ 

At TPP ,*?i^^J^*/S»N,N0S0RT2/S,N, NOD/S, N,N0P/0 $ 

qnp^ instead of PPF TO GET OUTPUT FREQUENCY LIRT 

^ , ,BGPDP,FOL,QPC,UPVC EST Yvrnp 

srN:SoSOR;°T* ’ ’ '^UPVc“/«FREQRESP;f 

^CTER 161 $ add LABEL FOR ERR0R3. 

LABEL ERR0R3 % 

aI~tpc * REMOVE ERROR 1 AND ERR0R2. 

ALTER 169 $ FORCED VIBRATION ERRORS 

PR?PAR« 

PR?PRRM //HcfYfcoSpMRsf 

JUMP FINIS $ 

PR?PAR« //T/«?cs?at?cLT" 
pRTPkRM /%°z^;:o2E/r'^" 

JUMP FINIS % 

PR?PAR« /%°/“oFRLn™^'^ deck. 

PRTPARM //O /*NOTRL* % 

JUMP FINIS $ 

PR?PAR« //f /“OFREQ^^ 'N CASE CONTROL. 

PRTPARM //O /*NOTIME* $ 

ENDALTER $ 
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TABLE 3.4 


CRAY NASTRAN Subprograms Modified to Inco rporate 
Updated Bladed Dieks Program 
Subroutine Subprograms (Total: 23) 


Subroutine 

Reason for modification 

Extent of changes 

APDB 

To permit the FLIST data 
block to be purged if the 
approach is not AERO 

Moderate 

CURV 

To give better diagnostics 

Minor 

CURVl 

To permit the computation of 
complex element stresses in 
material coordinate systems 
for the CQUADl/2 and 
CTRIAl/2 elements 

Moderate 

CURV2 

Same as CURVl as well as 
to correct a problem with 
the handling of incoming 
stresses/strains when the 
number of frequencies, 
times or subcases involved 
is more than one 

Extensive 

CYCT2A 

To correctly write a matrix 
trailer 

Minor 

FRD2I 

To avoid interpolation and 
and to permit just rear- 
rangement of the QHHL data 
in the FRRD2 module when the 
third parameter is specified 
as a negative number 

Extensive 

FVRSTl 

To set the Coriolis acceler- 
ation coefficient matrix and 
the base acceleration coeffi- 
matrix as square matrices 
rather than as square 
symmetric matrices 

Minor 

OFP 

To process complex element 
stresses in material coor- 
dinate systems for the 
CQUADl/2 and CTRIAl/2 
elements in both SORTl 
and S0RT2 output formats 

Moderate 

OFPIA 

Same as OFP 

Moderate 


(continued) 
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TABLE 3,4 (continued) 


Subroutine 

SDR3A 


AMGBIB 

GAUSS 

AMGBIC 


AMGTIB 


AMGTl 

AMGTIA 

AMGTIC 

AMGTIT 

SUBA 

SUBBB 

SUBC 

SUED 


Reason for modification Extent of changes 


To permit transformation 
of data from SORTl format 
to SORT2 format when the 
number of frequencies, 
times or subcases involved 
is one 

To improve unsteady pressure 
calculations in subsonic aero- 
dynamics without sweep effects 
To rectify and rearrange AJJ 
computations in supersonic 
aerodynamics without sweep 
effects 

To improve unsteady pressure 
calculations in subsonic 
aerodynamics with sweep effects 


Minor 


Minor 


Extensive 


Minor 


To incorporate sweep effects Extensive 

in supersonic aerodynamics 


AMGTID 


To permit interpolation of Moderate 

transonic aerodynamic matrices 
from subsonic and supersonic 
matrices with sweep effects 


(continued) 
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TABLE 3.4 (continued) 


Block Data Subprograms (Total: 3) 


Block data 

Reason for modification 

Extent of changes 

OFPIBD 

To process complex element 
stresses in material coor- 
dinate systems for the 
CQUADl/2 and CTRIAl/2 
elements in both SORTl 
and S0RT2 output formats 

Moderate 

OF2PBD 

To correct the output 
heading of complex 
element stresses in 
element coordinate 
systems for CTRIA2 
elements in SORTl 

Minor 


format 


SEMDBD 

To permit the processing 
of larger DMAPs than was 

Minor 


possible earlier 



Total Subprograms: 26 


Note: Subroutines APDB, FVRSTl, AMGBlBr GAUSS, AMGBIC, AMGTlB, 
AMGTl, AMGTIA, ^GTlC, AMGTlT, SUBA, SUBBB, SUBC, SUBD 
and AMGTID mentioned above were added to the code as part of 
the UNIVAC version of the Bladed Disks Program 
(see Table 3.2} • The rest of the subprograms existed in 
CRAY April 1984 NASTRAN, 
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ALTER 3 $ 

FILE UXVF=APPEND/PDT=APPEND/PD=APPEND % 

* PERFORM INITIAL ERROR CHECKS ON NSEGS, KMAX, KM IN AND K INDEX. 

COND ERRORCl , NSEGS $ IF USER HAS NOT SPECIFIED NSEGS. 

COND ERRORCl, KMAX $ IF USER HAS NOT SPECIFIED KMAX. 

COND ERRORCl, KM IN % IF USER HAS SPECIFIED NEGATIVE KM IN. 

PARAM //*NE*/KTEST/V, Y,KMAX/V, Y,KMIN=0 * 

COND LBL1KIND,KTEST $ 

% KMIN = KMAX 

PARAM //*ADD*/KINDEX/V,Y,KMAX/0 % SET KINDEX = KMAX (=KMIN) 
JUMP LBL2KIND $ 

LABEL LBLIKIND 

$ KMIN .NE. KMAX 

COND ERRORCl , KINDEX % IF USER HAS NOT SPECIFIED KINDEX. 

PARAM //*LT*/KTEST/V, Y,KINDEX/V,Y,KMIN « 

COND ERRORCl, KTEST $ 

PARAM //*GT*/KTEST/V, Y,KINDEX/V, Y,KMAX * 

COND ERRORCl, KTEST $ 

LABEL LBL2KIND * 

PARAM //*EQ*/CYCIOERR /V,Y,CYCI0=0 /O « 

COND ERRORCl , CYC I OERR * IF USER HAS NOT SPECIFIED CYC 10. 

PARAM //*DIV*/NSEG2 /V,Y, NSEGS /2 $ NSEG2 = NSEGS/2 

PARAM //*SUB*/KMAXERR /NSEG2 /V,Y,KMAX $ 

COND ERRORCl , KMAXERR * IF KMAX .6T. NSEGS/2 

$ CHECK FOR KINDEX = O 

PARAM //*EQ*/KTEST/V, Y,KINDEX/0 $ 

COND LBL3KIND, KTEST $ 

$ CHECK FOR KINDEX = NSEGS/2 (NSEGS EVEN. ONLY) 

PARAM //*ADD*/NSEGS1/V, Y,NSEGS/1 $ 

PARAM //*DIV*/NSEG21/NSEGSl/2 * 

PARAM //*EQ*/KEVEN/NSEG21/NSEG2 * 

PARAM //*EQ*/KNSEG2/NSEG2/V,Y, KINDEX $ 

PARAM //*EQ*/KTEST/KNSEG2/KEVEN $ 

COND LBL3KIND,KTEST % 

$ KINDEX IS -NE.O AND .NE. NSEGS/2 (NSEGS EVEN ONLY) 

PARAM //*ADD*/KTEST/1 /O * 

LABEL LBL3KIND $ 

PARAM //*GT*/KFLAG/KTEST/0 % 

% SET DEFAULTS FOR PARAMETERS. 

PARAM //*NOP*/V, Y,N0K:pRT=+1 /V, Y, LGKAD=-1 $ 

t CALCULATE OMEGA, 2*0MEGA AND 0MEGA**2 FROM RPS. SET DEFAULT RPS 
PARAMR //*MPY*/OMEGA /V,Y,RPS=0.0 /6. 283185 $ 

PARAMR //*MPY*/0MEGA2 /2.0 /OMEGA * 

PARAMR //*MPY*/OMEGASQR /OMEGA /OMEGA % 

$ GENERATE NORPS FLAG IF RPS IS ZERO. 

PARAMR //*EQ*//V, Y, RPS /O.O ////NORPS $ 

S MAKE SURE COUPLED MASSES HAVE NOT BEEN REQUESTED. 

PARAM //»NOT*/NOLUMP /V, Y, COUPMASS=-l % 

COND ERR0RC2, NOLUMP $ 

ALTER 21,21 $ ADD SLT TO OUTPUT FOR TRLG, 

6P3 GE0M3, EQEXIN, GE0M2 / SLT,GPTT / NOGRAV $ 

ALTER 23 % 

% SINCE MULTIPLE CONSTRAINTS ARE NOT ALLOWED EXECUTE GP4 NOW SO THAT 
% MfDRE ERROR CHECKS CAN BE MADE BEFORE ELEMENT GENERATION. 

$ ADD Y3 NEEDED FOR PSF RECOVERY IN SS62. 

PARAM //■»<M1--Y*/NSI: IP /O/O $ 
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GP4 CASECC , GEOM'4 , EQEX IN, 6PDT , BGPDT , CSTM , /RG , YS, USET , ASET/LUSET/ 

S, N, MPCFl/S, N,MPCF2/S, N, SINGLE/S, N, OMIT/S, N, REACT/S, N, NSKIP/ 

S, N, REPEAT /S, N, NOSET/S, N, NOL/S, N, NOA/C, Y, ASETOUT/S, Y, AUTOSPC $ 
PURGE GM, GMD/MPCFl/GO,GOD/OMIT/KFS,PSF,QPC/SINGLE $ 

* SUPORT BULK DATA IS NOT ALLOWED. 

PARAM //*NOT*/REACDATA /REACT * 

COND ERR0RC3, REACDATA % 

* EXECUTE DPD NOW SO CHECKS CAN BE MADE. ADD TRL TO OUTPUT DATA BLOCKS. 

DPD DYNAMICS, GPL, SIL, USET / GPLD, SILD, USETD, TFPOOL, DLT, PSDL, FRL, , 

TRL,EED,EQDYN / LUSET/S, N, LUSETD/NOTFL/S, N, NODLT/ 

S, N, NOPSDL/S, N, NOFRL/NONLFT/S, N, NOTRL/S, N, NOEED// 

S,N,NOUE $ 

% CHECK FOR EIGENVALUE EXTRACTION DATA 
COND ERR0RC7, NOEED * 

* MUST HAVE EITHER FREQ OR TSTEP BULK DATA. 

PARAM //*AND*/FTERR /NOFRL /NOTRL $ 

COND ERR0RC5,FTERR $ NO FREQ OR TSTEP BULK DATA. 

S ONLY FREQUENCY OR TSTEP IS ALLOWED IN THE CASE CONTROL 
PARAML CASECC //*DTI*/1/14//FREQSET * 

PARAML CASECC //*DTI*/1 /38//TIMESET « 

PARAM //*MPY*/FREQTIME /FREQSET /TIMESET * 

PARAM //*NOT*/FTERRl /FREQTIME ♦ 

PARAM //*LE*/NOFREQ /FREQSET /O * 

PARAM //*LE*/NOTIME /TIMESET /O * 

COND ERR0RC6, FTERRl % BOTH FREQ AND TSTEP IN CASE CONTROL DECK. 

* EPOINT BULK DATA NOT ALLOWED 

PARAM //*NOT*/EXTRAPTS /NOUE ♦ 

COND ERR0RC4,EXTRAPTS * 

* GENERATE DATA FOR CYCT2 MODULE. 

GPCYC GE0M4,EQDYN, USETD /CYCDD /CTYPE=ROT /S,N,NOGO % 

COND ERRORCl,NOGD t 

ALTER 28 $ 

PARAM //*NOP*/V, Y,KGGIN=-1 * 

COND JMPKGGIN,KGGIN * 

PARAM //*ADD*/NOKGGX/-l/0 ^ 

INPUTTl /KTOTAL, , , , /C, Y, LOCATION=-l/C, Y, INPTUNIT=0 % 

EQUIV KTOTAL, KGGX $ 

LABEL JMPKG6IN S 

ALTER 33 $ • . 

% PRE-PURGE DATA BLOCKS THAT WILL NOT BE GENERATED 
PARAM //*OR*/NOBMl /NOMGG /NORPS $ 

PURGE B1GG,M1GG /NOBMi % 

PURGE M2G6,M2BASEX6 /NOMGG ♦ 

ALTER 36 $ 

% GENERATE DATA BLOCKS FRLX, BIGG, MIGG, M2GG AND BASEGX. 

* GENERATE PARAMETERS FKMAX AND NOBASEX. 

FVRSTRl CASECC, BGPDT, CSTM, DIT, FRL, MGG, , / FRLX , BIGG, MIGG, 

M2GG, BASEXG, PDZERO, , /NOMGG/V, Y, CYCIO/V, Y, NSEGS/ 

V, Y,KMAX/S,N,FKMAX/V, Y,BXTID=-1/V, Y, BXPTID=-1/ 

V, Y, BYTID=-1/V, Y,BYPTID=-1/V, Y,BZTID=-1/ 

V, Y, BZPTID=-1/S, N, NOBASEX/NOFREQ/OMEGA $ 

PARAML FRLX //*PRESENCE*////NOFRLX $ 

COND LBLFRLX,NOFRLX $ 

EQUIV FRLX, FRL $ 

LABEL LBLFRLX % 

ALTER 43 t 

PARAM //*ADD*/NOBGG /NOBMI /O $ RESET NOBGG. 

ALTER 53 % 

% REDEFINE BGG AND KGG, 

COND LBLllA, NOBMI $ 

PARAMR //^COMPLEX*// 0MEGA2 /O.O/ CMPLXl t 

PARAMR //*SUB*/ M0ME6ASQ / 0.0 / OMEGASQR % 

PARAMR //*COMPLEX*// MOMEGASQ / 0.0 / CMPLX2 % 

BGG, BIGG / BGGl / (1.0, 0.0) / CMPLXl % 

EQUIV BGGl, BGG $ 

KGG, MIGG / KGGl / (1.0, 0.0) / CMPLX2 $ 


67 



£Q01V kGu1,KGG $ 

LABEL LBLllA 

ALTER 54,56 % GP4 HAS BEEN MOVED-UP. 
alter 88,88 % DPD HAS BEEN MOVED-UP. 

ALTER 113 $ PARAM AND EQUIV LOGIC DEPENDING ON LGKAD FOR FREQ OR TRAN 
//*AND»/KDEKA/N0UE/N0K2PP $ ^ TRm, 

LGKAD 1, LGKAD % BRANCH IN NOT FREQRESP. 

114 $ SEE ALTER 113 COMMENT. 

LGKAD2 * 

LGKAD 1 $ 

TAr^nn B2DD/N0A/K2PP, K2DD/NOA/MAA, MDD/MDEMA/ 

, kuu/KDc.KA % 

LGKAD2 $ 

ALTER 116,116 $ 

PARAMETERS GKAD, W3 AND W4 TO GKAD. 

GkAD * 


PARAM 

COND 

ALTER 

JUMP 

LABEL 

EQUIV 

LABEL 


ALTER 

COND 


GOD^-nn’ K4AA, K2PP, M2PP, B2PP/KDD, BDD, MDD, GMD, 
’ r B2DD/C , Y , GKAD=TRANRESP/*D I SP* /*D I RECT*/ 

C, Y , G 0. 0/C, Y , W3— 0. 0/C, Y, W4»=0. 0/N0K2PP/N0M2PP / 

N0B2PP/MPCF1/SINGLE/0MIT/N0UE/N0K4GG/ 

N0BGG/KDEK2/-1 * 

117 $ SEE ALTER 113 COMMENT. 

LGKAD3, LGKAD $ BRANCH IF NOT FREQRESP. 


118 $ SEE ALTER 113 COMMENT. 


JUMP 

LABEL 

EQUIV 

LABEL 


LGKAD4 $ 

LGKAD3 $ 

LGKAD4 % 

ALTER 119,123 % 

$ NEW SOLUTION LOGIC 

t UsfPSL^INS?EAD^orppp^k°rr? REQUESTED IN CASE CONTROL. 

CONf cS^fN^TIMr* frequency LIST. 

PAPAS'’ '"°R time-dependent loads. 


PARAM 
PARAM 
LABEL 
CASE 
PARAM 
TRLG 

SDRl 

SDRl 

PARAM 

COND 

REPT 

JUMP 

LABEL 

FVRSTR2 


EQUIV 

JUMP 

LABEL 


//*MPY*/REPEATT /I /-I « 

/I /O » INITIALIZE FOR SDRl. 

I HLGL.OOP ♦ 

/%MPy;/NCOL"%*lrr"®’ '^°E00P1 , 

??PD^^: PDfy?;SI:J’NOSE?%cJ:^^’“"”’ - ^^T. mss/ 

TP^^onT^” ” ” ” ’ ^ »PDT, /APPFLG/*DYNAMICS* ifc 
L,PD1 ,,,,,,,,, / ,PD , /APPFLG/*DYNAMICS* ^ 

//.ADD./APPFL6 /APPFLG /I * APPFLG=APPFLGfI 
TRLGDONE,REPEATT * nt-ri-LG+l. 

TRLGLOOP, 100 % 

ERR0R3 % 

TRLGDONE * 

^ 2 »F0LZ, REORDER 1, RE0RDER2,, ,, / 

V, Y, NSEGS/V, Y, CYCIO/S, Y, LMAX=-1/FKMAX/ 

S, N, FLMAX/S, N, NTSTEPS/S, N, NOROl /S, N, N0R02 t 
FRLZ,FRL // FOLZ.FOL % 

LBLFRL2 $ 

LBLTRLl * 


SELECTED IN CC. 


COND 
MPYAD 
EQUIV 
LABEL 
t FORM NEW LOADS. 


CASEXX HGPTn nix iroi i wHa atLtu I tu IN CC. 

“?rec;“/freq?/Ifr^^^^ ^ fpf,psf,pdf,fol,phfdum / 

PPFrPoZERS?'"/''pPFX IF FRLX WAS GENERATED. 

PPFX,PPF * 

LBLFRLXl * 


COND 

MPYAD 

ADD 

EQUIV 

COND 


LBLFRLl , NOBASEX $ 
M2GG,BASEXG, / M2BASEXG /O $ 
PPF,M2BASEXG / PPFl /(1 
PPF 1 , PPF $ 

LBLBASEl , NOSET $ 


0,0.0) /(-I. 0,0.0) $ 
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SSG^ USETD, GMD, YS, KFS, GOD, , PPF / , PODUMl , PSF 1 , PDh 1 % 
EQUIV PSF1,PSF // PDF 1, PDF 

LABEL LBLBASEl $ 

LABEL LBLFRLl * 

EQUIV PPF, PDF/NOSET $ 

* LOADS ARE FREQUENCY-DEPENDENT 

♦ PERFORM CYCLIC TRANSFORMATION ON LOADS IF CYCIO=+l 
PARAML PDF //^TRAILER*/ 1 /PDFCOLS t 

$> CALCULATE THE NUMBER OF LOADS FOR CYCIO«=-l. 

PARAM //*DIV*/NLOAD /PDFCOLS /FKMAX $ NLOAD = NF/FKMAX 
EQUIV PDF, PXF/CYCIQ $ 


COND LBLPDONE,CYCIO $ 

* CALCULATE THE NUMBER OF LOADS FOR CYC I 0=1. 


PARAM 

CYCTl 


//*DIV*/NLOAD /PDFCOLS /V,Y,NSEGS NLOAD >= NF/NSEGS 
PDF / PXF,GCYCF1 /CTYPE /*FORE*/V, Y, NSEGS=-1 / 
V,Y,KMAX=-1 / NLOAD /S,N,N060 « 

COND ERR0RC1,N0G0 % 

JUMP LBLPDONE % 

LABEL LBLFRL2 $ 

% LOADS ARE TIME-DEPENDENT 

PARAM //*NOT*/NOTCYCIO /V,Y,CYCIO $ 

* BRANCH DEPENDING ON VALUE OF CYC 10 

COND LBLTRL2,N0TCYCI0 $ 

$ CYCI0=-1 


EQUIV 
COND 
MPYAD 
LABEL 
CYCTl 

COND 
EQUIV 
COND 
MPYAD 
LABEL 
EQUIV 
JUMP 
LABEL 
* CYCIO = +1 
MPYAD 


PD,PDTRZ1/N0R01 $ 

LBLR01A,N0R01 % 

PD, REORDER 1, / PDTRZ 1 / 0 * 

LBLROIA $ 

PDTRZ 1 / PXTRZ1,GCYCF2 /CTYPE/*FORE*/NTSTEPS/ 
V, Y, LMAX /FKMAX /S, N, NOGO « 

ERR0RC1,N0G0 $ 

PXTRZ1,PXFZ1/N0R02 $ 

LBLR02A, N0R02 * 

PXTRZ1,RE0RDER2, 

LBLR02A * 

PXFZ1,PXF1 $ 

LBLTRL3 « 

LBLTRL2 $ 


/ PXFZl /O $ 


PD,RE0RDER1, / PDTRZ2 / O * 

PDTRZ2 /PXTRZ2,GCYCF3 /CTYPE/*FORE*/NTSTEPS/V, Y , LMAX/ 

V, Y,NSEGS/S,N,NOGO $ 

ERR0RC1,N0G0 * 

PXTRZ2,PXTR2/N0R02 $ 

LBLR02B,N0R02 $ 

PXTRZ2,RE0RDER2, / PXTR2 /O t 
LBLR02B $ 

PXTR2 / PXFZ2,GCYCF4 / CTYPE/*FDRE*/V, Y, NSEGS/V, Y, KMAX/ 
FLMAX/S, N, NOGO $ 

ERR0RC1,N0G0 $ 

PXFZ2,PXFI $ 

LBLTRL3 % 

% TIME-DEPENDENT LOADS ARE REAL. MAKE LOADS COMPLEX TO CORREPPOMn 
i DEFINE NLOAD FOR CYCT2. 

PARAM //*ADD*/NLOAD /FLMAX /O % NLOAD = FL MAX 
LABEL LBLPDONE % f-L-n«X 

$ 

$ INITIALIZE UXVF IF KM IN IS NOT ZERO. 

$ 


CYCTl 

COND 

EQUIV 

COND 

MPYAD 

LABEL 

CYCTl 

COND 

EQUIV 

LABEL 


PARAM //*ADD-« /KMINL /V, Y, KINDEX=-1 /-I $ 

COND NOKMINL, KMINL $ 
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PARAM 

LABEL 

CYCT2 

COND 

ADD 

PRTPARM 

CYCT2 

PRTPARM 

COND 

PARAM 

REPT 

LABEL 

COND 

PRTPARM 

LABEL 

CYCT2 

COND 

CYCT2 

COND 

CYCT2 

COND 

READ 

OFP 

COND 

OFP 

COND 

CYCT2 

COND 

SDRl 

SDR2 

PLOT 

PRTMSG 

LABEL 

GKAM 


PARAML 

SMPYAD 

SMPYAD 

SMPYAD 

MPYAD 

EQUIV 

COND 

% KINDEX 
APDB 


AMG 

AMP 

JUMP 
LABEL 
$ KINDEX 
CYCT2 

COND 

APDB 


//*ADD*/KMINV /O /O ^ 

KM INLOOP % 

CYCDD, , ,PXF, , /i,PKFZ,, / *FORE*/V, Y, NSEGS/ 
KMINV/CYCSEQ/NLOAD/S,N,NOGO $ 

ERRORCl,NOGO * 

PKFZ, / UKVFZ / <0.0, 0.0) $ 

//0/*KINDEX* % 

CYCDD, ,, UKVFZ, , /,,UXVF,, /*BACK*/V, Y, NSEGS/ 
KMINV/CYCSEQ/NLOAD/S,N,NOGO % 

//0/*KINDEX* * 

ERR0RC1,N0G0 % 

//*ADD*/KMINV /KMINV /I % 

KM INLOOP, KM I NL $ 

NOKMINL $ 

NOKPRT,NOKPRT * 

//0/*KINDEX* $ 

NOKPRT ♦ 

CYCDD, KDD, MDD, , , /KKKF, MKKF, , , /*FORE*/V, Y, NSEGS / 

V, Y,KINDEX/CYCSEQ/NLOAD/S,N,NOGO « 

ERRORCi,NOGO $ 

CYCDD, BDD, , PXF, , /BKKF,,PKF,, / *FORE*/V, Y, NSEGS/ 

V, Y, KINDEX/CYCSEQ/NLOAD/S, N, N060 $ 

ERR0RC1,N0G0 $ 

CYCDD, KAA,MAA, , , /KKK,MKK, , , /*FORE*/V, Y, NSE6S/V, Y, KINDEX/ 
C YCSEQ— “ 1 / 1 /S , N , NOGO $ 

ERR0RC1,N0G0 % 

KKK,MKK, , ,EED, » CASECC/LAMK, PHIK, MIK,OEIGS/*MODES*/ 
S,N,NEIGV i 


OEIGS, , , , , //S,N,CARDNO % 
FINIS, NEIGV $ 


LAMK, , , , , //S,N,CARDNO $ 

NOPLOT, JUMPPLOT % 

CYCDD, , , ,PHIK,LAMK/, , , PHIA, LAMA/*BACK*/V, Y, NSEGS/V, Y, KINDEX/ 
CYCSEQ/1/S,N,N0G0 $ 

ERRORCi,NOGO $ 

USET, ,PHIA, , ,GO,GM, ,KFS, , /PHIG, ,QG/1/*REIG* $ 
CASECC,CSTM,MPT,DIT,EQEXIN,SIL, , » BGPDP, LAMA, QG, PHIG, EST, , / 

, OQGl , OPHIG, OESl , OEFl , PPHIG/*REIG* % 

PLTP AR , GPSETS , ELSETS , CASECC , BGPDT , EQE X I N , S I P , , PPH I G , GPECT , / 
PLOTXX/NSIL/LUSEP/JUMPPLOT/PLTFLG/S, N, PFILE * 

PLOTXX// * 

NOPLOT * 

USETD, PHIK, MIK, LAMK, DIT, M2DD, B2DD, K2DD, CASECC/MDUM, BDUM, 
KDUM, PHIKH/NOUE/C, Y, LM0DES=6/C, Y, LFREQ=0. 6/C, Y, HFREQ=-1 . 0/ 
N0M2PP/N0B2PP/N0K2PP/N0NCUP/S,N,FM0DE=0 * 
PHIKH//*TRAILER*/1/S,N,NM0DES $ 

PHIKH,MKKF,PHIKH, , , /MHH/3////1 % 

PHIKH,KKKF,PHIKH, , ,/KHH/3////l * 

PHIKH, BKKF, PHIKH, , , /BHH/3////1 $ 

PHIKH,PKF, /PHF/1 * 


MHH, MKKF/ /BHH, BKKF //KHH, KKKF//PHF, PKF $ 
KLABELi,KFLAG % 

IS EITHER 0 OR NSEGS/2 < NSEGS EVEN ONLY) 

EOT , USET , BGPDT , CSTM, EQEX IN, GM, GO/AERO, ACPT, , GTKA, / 
S,N,NK/S,N,NJ/V, Y,MINMACH/V,Y,MAXMACH/V, Y, IREF// 
NMODES/V, Y, KINDEX * 

AERO, ACPT/AJJL, SKJ, DIJK, D2JK/NK/NJ/1 * 

A J JL , SKJ , D UK , D2 JK , GTKA , PH I KH , , , USETD , AERO/QHHL , , / 
NOUE/1 % 

KLABEL2 % 

KLABELl % 

IS .NE.O AND .NE. NSEGS/2 (NSEGS EVEN ONLY) 

CYCDD, , , , PHIKH, LAMK/, , , PHIAH, LAMAH/*BACK*/V, Y, NSEGS/ 
V, Y, KINDEX /CYCSEQ/1/S,N, NOGO % 

ERRORCljNOGO * 


LDT, USET , BGPDT , CSTM, EQEX IN, GM, GO/AERO, ACPT, , 


GTKAtPVECT/ 
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AMG 

PARTN 

AMP 

APDB 


PARTN 

AMP 

ADD 

LABEL 


OUN, ixfN/ Q, IV, NJ/ V, y 1 niNnHUH/ V, Y, n«XMAUH/V, Y, IREF/*COSINE* / 
NMODES/V,Y,K INDEX % 

AERO, ACPT/AJJL, SKJ, DIJK, D2JK/NK/NJ/ 1 * 
PHIAH,PVECT,/PHIAC, , ,/l * 

AJJL, SKJ , DIJK, D2JK, GTKA, PHIAC, , , USETD, AERO/QHHLC, , / 
NOUE/1 * 

EDT , USET , BGPDT , CSTM, EQEXIN, GM, GO/AERO, ACPT, , GTKA, PVECT/ 
S, N, NK/S, N, NJ/V, Y, MINMACH/V, Y, MAXMACH/V, Y, IREF/*SINE*/ 
NMODES/V, Y,KINDEX « 

PHIAH,PVECT, /PHIAS, , , /I $ 

AJJL, SKJ, DIJK, D2JK, GTKA, PHI AS, ,, USETD, AERO/QHHLS, , / 
NOUE/1 % 

QHHLC,QHHLS/QHHL * 

KLABEL2 ♦ 


$ SOLUTION 


FRRD2 

DDRl 

EQUIV 

CYCT2 

COND 

EQUIV 

COND 

CYCTl 

LABEL 

COND 

EQUIV 

COND 

CYCTl 

LABEL 


/*BACK*/V, Y, NSEGS/V, Y, KINDEX/ 


KKKF,BKKF,MKKF,QHHL,PKF,FOL/UKVF/V, Y, BOV/V,Y,Q/-l 0 $ 
UKVF,PHIKH/UKKVF ♦ 

UKKVF,UKVF % 

CYCDD, , ,UKVF, , /,,UXVF,, 

CYCSEQ/NLOAD/S,N,NOGO % 

ERRORCl,NOGO $ 

UXVF,UDVF / CYCIO * 

LCYC3,CYCI0 % IF CYCIO .GE. 0 THEN TRANSFORM TO PHYSICAL 
UXVF / UDVF,GCYCB1 / CTYPE/*BACK*/V, Y, NSEGS/V, Y, KMAX/ 
NLOAD * 

LCYC3 ^ 

LBLTRL4, NOTIME * 

PXF,PDF2 / CYCIO % 

LCYC4, CYCIO * IF CYCIO .GE. 0 THEN TRANSFORM TO PHYSICAL 
PXF / PDF2,6CYCB2 / CTYPE/*BACK*/V, Y, NSEGS/V, Y, KMAX/ 

NLOAD % 

LCYC4 $ 


$ IF LOADS WERE TIME-DEPENDENT THEN RECOVER PPF AND PSF FROM PXF. 


SDRl 

SSG2 

EQUIV 

LABEL 

VDR 


USETD, , PDF2, ,, GOD, 6MD, ,, , / PPFZ,, /I /^DYNAMICS* $ 

USETD, GMD,YS,KFS, GOD, ,PPFZ / , PODUM, PSFZ, PLDUM * 

PPFZ, PPF // PSFZ, PSF t 
LBLTRL4 t 

CASEXX, EQDYN, USETD, UDVF , FOL, XYCDB, /OUDVCl , /*FREQRESP*/ 
*DIRECT*/S, N, N0S0RT2/S, N, NOD/S, N, NOP/FMODE ♦ 

ALTER 139,139 $ USE FOL INSTEAD OF PPF TO GET OUTPUT FREQUENCY LIST. 

CASEXX , CSTM, MPT , DIT, EQDYN, SILD, , , BGPDP, FOL,QPC, UPVC, EST, XYCDB, 
PPF/OPPCl , OQPCl , OUPVCl , OESCl , OEFCl , PUPVCl /*FREQRESP*/ 
S,N,N0S0RT2 % 

OESCl, MPT, CSTM, EST, SIL,GPL/0ESC1M, /I * 


sdr: 


CURV 


ALTER 141,142 $ 


SDR3 


OFP 


OPPCl , OQPCl , OUPVCl , OESCl , OEFCl , 0ESC1M/0PPC2, 0QPC2,0UPVC2, 
0ESC2,0EFC2,0ESC2M $ 

OPPC2, 0QPC2, OUPVC2,OEFC2,OESC2, OESC2M//S, N, CARDNO $ 


ALTER 153,153 * 

OPPCl, OQPCl, OEFCl, OESCl, 0ESC1M//S,N, CARDNO * 

ALTER 161 t ADD LABEL FOR ERR0R3. 

LABEL ERR0R3 % 

164,167 % REMOVE ERROR 1 AND ERR0R2. 

ALTER 169 $ FORCED VIBRATION ERRORS 

LABEL ERRORCl ♦ CHECK NSEGS, KMAX AND OTHER CYCLIC DATA 
PRTPARM //-7 /*CYCSTATICS* % 

LABEL ERR0RC2 $ COUPLED MASS NOT ALLOWED. 

PRTPARM //O /C, Y, COUPMASS % 

JUMP FINIS % 

LABEL ERR0RC3 $ SUPORT BULK DATA NOT ALLOWED, 

PRTPARM //-6 /^CYCSTATICS* % 

LABEL ERRORC4 % EPOINT BULK DATA NOT ALLOWED 
PRTPARM //O /*NOUE* % 

JUMP FINIS % 

* neither FREQ OR TSTEP WERE IN BULK DATA DECK. 
ERTPARM //O /*NOFRL* $ 


IS 

QUALITY 


QUAury 
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Ph:TJ-'Ati:fi 

JUMP 

LABEL 

PRTPARM 

PRTPARM 

JUMP 

LABEL 

PRTPARM 

ENDALTER 


//O /*rviU’lRL* % 

FI^4IS $ 

/%°/^N0*REqS^2 SELECTED :n CASE CONTROL. 

//O /*NOTIME* ♦ 

FINIS $ 

ERRORC7 $ NO EIGENVALUE EXTRACTION DATA 
//-2/*CYCM0DES* % 

$ 
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TABLE 3 .6 DMA? Sequence of DISP APP RF 8 ALTERed by 
MFVAAET ALTER Pack^ 


1 

2 

~r 


3 

3 

3 



3 



T 


4 


6 

7 

8 
9 

10 
1 1 


14 
1 5 


BESIN 

PRECHK 

FILE 

FILE 

COND 

COND 

COND 

PARAM 

COND 

PARAM 

JUMP 

LABEL 

COND 

PARAM 

COND 

PARAM 

COND 

LABEL 

PARAM 

COND 

PARAM 

PARAM 

COND 

PARAM 

COND 

PARAM 

PARAM 

PARAM 

PARAM 

PARAM 

COND 

PARAM 

LABEL 

PARAM 

PARAM 

PARAMR 

PARAMR 

PARAMR 

PARAMR 

PARAM 

COND 

PARAM 

GPl 

PLTTRAN 

PURGE 

COND 

GP2 

PARAML 

PURGE 

COND 

PLTSET 


PR T MSG 
PARAM 


DISP 08 - DIRECT FREQUENCY/RANDOM RESPONSE ANALYS I S-APR. 1984 $ 
ALL $ 

KGGX=TAPE/KGG=TAPE/GOD=SAVE/GMD=SAVE/MDD=SAVE/BDD=SAVE $ 
UXOF=APPEND/PDT=APPEND/PD~APPEND % 

ERR0RC1,NSEGS * IF USER HAS NOT SPECIFIED NSEGS. 

ERRORCl.KMAX $ IF USER HAS NOT SPECIFIED KMAX. 

ERR0RC1,KMIN $ IF USER HAS SPECIFIED NEGATIVE KMIN. 

//*NE*/KTEST/V, Yi KMAX/V, Y, KMIN=0 $ 

LBL1KIND,KTEST $ 

//*ADD*/KINDEX/V, Y,KMAX/0 * SET KINDEX = KMAX <= KMIN) 
LBL2KIND $ 

LBLIKIND 

ERRORCl , KINDEX $ IF USER HAS NOT SPECIFIED KINDEX. 
//*LT*/KTEST/V,Y,KINDEX/V, Y,KMIN $ 

ERRORCl, KTEST $ 

//*GT*/KTEST/V, Y, KINDEX/V, Y, KMAX $ 

ERRORCl, KTEST $ 

LBL2KIND $ 

//*EQ*/CYCIOERR /V,Y,CYCI0=0 /O % 

ERRORCl, CYC lOERR % IF USER HAS NOT SPECIFIED CYCIO 
//*DIV*/NSEG2 /V,Y, NSEGS /2 $ NSEG2 = NSEGS/2 
//*SUB*/KMAXERR /NSEG2 /V,Y,.KMAX S 
ERRORCl , KMAXERR % IF KMAX .GT. NSEGS/2 
//*EQ*/KTEST/V, Y,KINDEX/0 $ 

LBL3KIND,KTEST $ 

//*ADD*/NSEGS1/V, Y,NSEGS/1 $ 

//*DIV*/NSE621/NSEGSl/2 % 

//*EQ*/KEVEN/NSEG21/NSEG2 $ 

//*EQ*/KNSEG2/NSEG2/V,Y, KINDEX $ 

//*EQ*/KTEST/KNSEG2/KEVEN % 

LBL3KIND,KTEST $ 

//*ADD*/KTEST/l/0 $ 

LBL3KIND $ 

//*GT*/KFLAG/KTEST/0 $ 

//*NOP*/V, Y,N0KPRT=+1 /V, Y, LGKAD=-1 $ 

//*MPY*/OMEGA /V,Y,RPS=0.0 /6. 283185 $ 

//*MPY*/0ME6A2 /2.0 /OMEGA $ 

//*MPY*/OMEGASQR /OMEGA /OMEGA % 

//*EQ*//V, Y,RPS /O.O ////NORPS $ 

//*NOT*/NOLUMP /V, Y,C0UPMASS=-1 $ 

ERR0RC2, NOLUMP % 

//*MPY*/CARDN0/0/0 % 

GEOMl , GE0M2, /GPL, EQEXIN, GPDT, CSTM, BGPDT, SIL/S, N, LUSET/ 
S,N,N0GPDT/ALWAYS=-1 $ 

BGPDT , SIL/BGPDP, SIP/LUSET/S, N, LUSEP $ 

USET, GM, GO, KAA, BAA, MAA, K4AA, KFS, PSF , QPC, EST, ECT, PLTSETX , PLTPAR, 
GPSETS,ELSETS/NOGPDT $ 

LBL5,N0GPDT $ 

GEGM2, EQEXIN/ECT $ 

PCDB//*PRES*////JUMPPLOT * 

F-LTSETX, PLTPAR, GPSETS, ELSETS/ JUMPPLOT $ 

P1,JUMPPL0T $ 

PCDB, EQEXIN, ECT /PLTSET X , PLTPAR, GPSETS , ELSETS/S , N , NS I L/ 

S,N, JUMPPLOT t. 

PLISETX// $ 

//*MPY#/PLTFLG/1 /I $ 
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16 

PARAM 

17 

COND 

18 

PLOT 

19 

PRTMSG 

20 

LABEL 

21 

GP3 

22 

TAl 


PURGE 

23 

PARAM 

23 

GP4 


23 

PURGE 

23 

PARAM 

23 

COND 

23 

DPD 


23 

COND 


PARAM 

23 

COND 

23 

PARAML 

23 

PARAML 

23 

PARAM 

23 

PARAM 

23 

PARAM 

23 

PARAM 

23 

COND 

23 

PARAM 

23 

COND 

23 

GPCYC 


COND 

24 

COND 

25 

PARAM 

26 

PARAM 

27 

PARAM 

28 

PARAM 

28 

PARAM 

28 

COND 

28 

PARAM 

28 

INPUTT 1 

28 

EQUIV 

28 

LABEL 

29 

EMG 


30 

PURGE 

31 

COND 

32 

EMA 

33 

LABEL 

33 

PARAM 

33 

PURGE 

3.3 

PURGE 

34 

COND 

“T cr 

EMA 

36 

LABEL 

36 

FVRSTI 


/ 7*KPY*/Ph iLE/U/O % 

P 1 , JUMPPLOT $ 

PLTPAR, GPSETSi ELSETS, CASECC, BGPDT t EQEX IN, SIL, , ECT, , /PLOT XI / 
NSIL/LUSET/S,N, JUf1PPLOT/S,N,PLTFLG/S,N,PFILE 
PLOTXl//$ ' - 

PI * 

GE0M3,EQEXIN,6E0M2 / SLT,GPTT / NOGRAV % 

ECT, EPT, BGPDT, SIL, GPTT,CSTM/EST,GEI,GPECT, , / 

LUSET/S, N, NOSIMP/l/S, N, N06ENL/S, N, GENEL $ 

K4GG, GPST , OGPST , MGG, BBG, K4NN, K4FF , K4AA, MNN, MFF, MAA, BNN, BFF, BAA, 
KGGX/N0SIMP/06PST/GENEL % 

//*MPY*/NSKIP /O/O $ 

CASECC, GE0M4, EQEXIN, GPDT , BGPDT , CSTM, /RG, YS, USET, ASET/LUSET/ 

S , N , MPCF 1 / S , N , MPCF2 /S , N , S I NGLE/S , N , OM I T /S , N , RE ACT /S , N , NSK I P/ 

S, N, REPEAT /S, N, NOSET /S, N, NOL/S, N, NQA/C, Y, ASETOUT/S, Y, AUTOSPC $ 
GM,GMD/MPCFl/GO,GOD/OMIT/KFS,PSF,QPC/SINGLE % 

//*NOT*/REACDATA /REACT % 

ERR0RC3,REACDATA * 

DYNAMICS, GPL, SIL, USET / GPLD, SILD, USETD, TFPQQL, DLT, PSDL, FRL, , 
TRL, EED, EQDYN / LUSET/S, N, LUSETD/NOTFL/S, N, NODLT/ 

S, N, NOPSDL/S, N, NOFRL/NONLFT /S, N, NOTRL/S, N, NOEED// 

S,N,NOUE % 

ERRORC7, NOEED % 

//*AND*/FTERR /NOFRL /NOTRL $ 

ERR0RC5,FTERR % NO FREQ OR TSTEP BULK DATA. 

CASECC //*DTI*/1/14//FREQSET * 

CASECC //*DTI*/l/38//TIMESET % 

//*MPY*/FREQTIME /FREQSET /TIMESET % 

//*NOT*/FTERRl /FREQTIME % 

//*LE*/NOFREQ /FREQSET /O $ 

//*LE*/NOTIME /TIMESET /O $ 

ERR0RC6,FTERR1 $ BOTH FREQ AND TSTEP IN CASE CONTROL DECK. 
//*NOT*/EXTRAPTS /NOUE $ 

ERR0RC4,EXTRAPTS % 

GE0M4, EQDYN, USETD /CYCDD /CTYPE=ROT /S,N,NOGO $ 

ERRORCl,NOGO $ 

LBLl,NOSIMP $ 

//*ADD*/NOKGGX/l/0 $ 

//*ADD*/NOMGG/l/0 $ 

//*ADD»/NQBGG=-l/l/0 $ 

//*ADD*/N0K4GG/1 /O % 

//*NOP*/V, Y,KGGIN=-1 % 

JMPKGGIN,KGGIN % 

//»ADD*/NOKGGX/-l/0 $ 

/KTOTAL, , , , /C, Y, LOCATION=-l/C, Y, INPTUNIT=0 * 

KTOTAL,KGGX $ 

JMPKGGIN $ 


EST , CSTM, MPT , DIT , GE0M2, /KELM, KDICT, MELM, MDICT, BELM, EDICT/ 
S, N, N0KG6X/S, N, N0MG6/S, N, NOBGG/S, N, N0K4GG//C, Y, COUPMASS/ 

C, Y, CPBAR/C, Y , CPRQD/C, Y , CPQUADl /C, Y, CPQUAD2/C, Y, CPTRI A1 / 

C, Y,CPTRIA2/C, Y, CPTUBE/C, Y, CPQDPLT/C, Y, CPTRPLT/C, Y, CPTRBSC 
GPST/NOKGGX/MGG/NOMGG $ 

LBLKGGX,NOKGGX $ 

GPECT, KDICT, KELM/KGGX, GPST $ 

LBLKGGX % 

//*OR*/NOBMl /NOMGG /NORPS % 

B1GG,M1GG /NOBMl $ 

M2GG,M2BASEXG /NOMGG % 

LBLMGG, NOMGG % 


GPECT , MDICT , MELM/MGG, /-i /C, Y, WTMASS=1 . O % 

LBLMGG $ 

CASECC, BGPDT, CSTM, DIT, FRL, MGG, , / FRLX , BIGG, MIGG, 
M2GG, BASEXG, PDZERO, , /NOMGG/V, Y, CYCIO/V, Y, NSEGS/ 
V, Y, KMAX/S, N, FKMAX/V, Y, BXT I D=-l /V, Y , BXPT I D=- 1 / 

V, Y, BYTID=-1 /V, Y, BYPTID=-1/V, Y, BZTID=-1 / 

V, Y, BZPTID=-1/S, N, NOBASEX/NOFREQ/OMEGA % 


OF POOf? quality 
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36 

PARAML 

36 

COND 

36 

EQUIV 

36 

LABEL 

37 

COND 

38 

EMA 

39 

LABEL 

40 

COND 

41 

EMA 

42 

LABEL 

43 

PURGE 

43 

PARAM 

44 

PURGE 

45 

COND 

46 

COND 

47 

GPWG 

48 

OFP 

49 

LABEL 

50 

EQUIV 

51 

COND 

52 

SMA3 

53 

LABEL 

53 

COND 

53 

PARAMR 

53 

PARAMR 

53 

PARAMR 

53 

ADD 

53 

EQUIV 

53 

ADD 

53 

EQUIV 

53 

LABEL 

57 

COND 

58 

COND 

59 

F’ARAM 

60 

COND 

61 

GPSP 

62 

OFP 

63 

LABEL 

64 

EQUIV 

65 

COND 


66 MCE! 

67 MCE2 

68 LABEL 

69 EQUIV 

70 COND 

7 1 SCE 1 

72 LABEL 

73 EQUIV 

74 EQUIV 

75 EQUIV 

76 EQUIV 

77 COK'D 

78 SMPl 

79 COND 

80 SMP2 

81 LABEL 

82 COND 

83 SMP2 

84 LABEL 

85 COND 

86 SMP2 

87 LABEL 

89 EQUIV 

90 PARAM 

9 1 PARAM 

92 BMG 


/O $ RESET NOBGG. 


f^RllX //*FRESt:NCE*////NOFRLX * 

LBLFRLX,NOFRLX * 

FRLXtFRL * 

LBLFRLX % 

LBLBGG, NOBGG $ 

GPECT,BDICT,BELM/B66, $ 

LBLBGG ♦ 

LBLK46G,N0K4GG * 

GPECT , KDICT , KELM/K4GG, /N0K4GG $ 

LBLK4G6 % 

MNN, MFF, MAA/NOMGG $ 

//*ADD*/NOBGG /NOBMl 
BNN,BFF,BAA/NOBGG $ 

LBL1,GRDPNT % 

ERR0R4,N0MGG % 

LBLl $ 

KGGX,KGG/NOGENL $ 

LBL11,N0GENL $ 

GEI f KGGX/KG6/LUSET /N06ENL/N0S IMP t 
LBLll $ 

LBLl lA, NOBMl $ 

//♦COMPLEX*// 0MEGA2 /O.O/ CMPLXl ♦ 

//*SUB*/ MOMEGASQ / 0.0 / OMEGASQR $ 

//♦COMPLEX*// MOMEGASQ / 0.0 / CMPLX2 $ 

BGGjBlGG / BGGl / (1.0 ,0.0) / CMPLXl $ 

BGG 1 <1 BGG $ 

KGG»M1GG / KGGl / (1. 0,0.0) / CMPLX2 t 
KGG1,KGG * ^ 

LBLllA 

LBL4,GENEL $ 

LBL4,N0SIMP ^ 

//♦EQ*/GPSPFLG/AUT0SPC/0 $ 

LBL4,GPSPFLG $ 

GPL, GPST, USET , SIL/OGPST /S, N, NOGPST $ 

OGPST, , , , , //S,N,CARDNO $ 

LBL4 * 

BGe,BNN/MPCFl/K4eG,K4NN/MPCFl , 

USET,R6/GM $ 

LBL2 ' KNN , MNN , BNN , K4NN $ 

LbS’ K4NN/KFF, KFS, , MFF, BFF, K4FF $ 

KFF,KAA/0MIT % 

MFF,MAA/OMIT $ 

BFF,BAA/0MIT % 

K4FF, K4AA/0MIT $ 

LBL5,0MIT * 

USET, KFF, , , /GO, KAA, KOO, LOO, , , , , $ 

LBLM,NOMGG $ 

USET, 60, MFF/MAA * 

LBLM % 

LBLB, NOBGG * 

USET,GO,BFF/BAA $ 

LBLB $ 

LBL5,N0K4GG $ 

USET,G0,K4FF/K4AA $ 

LBL5 $ 

GO, GOD/NOUE/GM, GMD/NOUE % 

//*ADD*/NEVER/l/0 t> 

//♦MPY*/REPEATF/-1/1 $ 

MATPOOL, BGPDT, EQEX IN , CSTM/BDPOOL/S , N, NOKBFL/S, N, NOABFL/ 


$ 


DRlG^r-^AL 
OF POvjR 


PAGE IS 
QUALITY 
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93 

PARAM 

94 

PURGE 

95 

COND 

96 

MTRXIN 

97 

LABEL 

98 

PURGE 


99 

CASE 

100 

MTRXIN 

101 

PARAM 

102 

PARAM 

103 

EQUIV 

104 

COND 

105 

ADDS 

106 

COND 

107 

TRNSP 

108 

ADD 

109 

LABEL 

110 

PARAM 

111 

PARAM 

112 

PARAM 

113 

PURGE 

1 13 

PARAM 

1 13 

COND 

114 

EQUIV 

114 

JUMP 

114 

LABEL 

114 

EQUIV 

114 

LABEL 

115 

COND 

1 16 

GKAD 


117 LABEL 

117 COND 

118 EQUIV 
118 JUMP 
118 LABEL 
118 EQUIV 
118 LABEL 
123 COND 
123 PAR AM 
123 PARAM 
1 23 LABEL 
123 CASE 
123 PARAM 
123 TRLG 


123 
123 
123 
123 
123 
1 23 
123 
123 


SDRl 

SDRl 

PARAM 

COND 

REPT 

JUMP 

LABEL 

F'v/RSTR2 


S,N,MFACT $ 


//*AND*/NOFL/NOABFL/NOKBFL $ 

KBFL/NOKBFL/ ABFL/NOABFL % 

LBL13,NOFL * 

, BDPOOL, EQDYN, , /ABFL, KBFL, /LUSETD/S, N, NOABFL/S, N, NOKBFL/ 
0 % 

LBL13 ♦ 


OUDVCl , OUDVC2, XYPLTFA, OPPCl , OQPCl , OUPVCl , OESCl , OEFCl , 0PPC2, 
0QPC2, OUPVC2, 0ESC2, 0EFC2, XYPLTF, PSDF, AUTO, XYPLTr7 
K2PP,M2PP, B2PP, K2DD, M2DD, B2DD/NEVER % 


CASECC, PSDL/CASEXX/*FREQ*/S, N, REPEATF/S, N, NOLOOP t 
CASEXXrMATPOOL, EQDYN, , TFP00L/K2DPP, M2DPP, B2PP/LUSETD/S* N. 
N0K2DPP/S,N,N0M2DPP/S,N,N0B2PP $ uubtlD/S,N, 

//*AND*/N0M2PP/N0ABFL/N0M2DPP $ 

//*AND*/N0K2PP/N0FL /NOK2DPP $ 

K2DPP, K2PP/N0FL/M2DPP, M2PP/N0ABFL % 

LBLFL2,N0FL $ 

ABFL, KBFL, K2DPP, , /K2PP/ (-1.0, 0.0) « 

LBLFL2,N0ABFL $ 

ABFL/ABFLT $ 

ABFLT, M2DPP/M2PP/MFACT $ 

LBLFL2 $ 


//*AND*/BDEBA/N0UE/N0B2PP * 

//*AND*/KDEK2/NQGENL/N0SIMP $ 

//#AND*/MDEMA/N0UE/N0M2PP % 
K2DD/N0K2PP/M2DD/N0M2PP/B2DD/N0B2PP $ 
//*AND*/KDEKA/N0UE/N0K2PP t 
LGKAD1,L6KAD $ BRANCH .IN NOT FREQRESP. 

M2PP, M2DD/N0A/B2PP, B2DD/N0A/K2PP, K2DD/N0A/ 

MAA, MDD/MDEMA/BAA, BDD/BDEBA % 

LGKAD2 $ 

L6KAD1 $ 

KAA^KDD/KDEKa'^^^^^’ B2DD/N0A/K2PP, K2DD/N0A/MAA, MDD/MDEMA/ 

LGKAD2 % 

LBL18,N0GPDT % 

MAA, K4AA, K2PP, M2PP, B2PP/KDD, BDD, MDD, GMD, 
GOD, K^DD, M2DD, B2DD/C, Y , GKAD=TRANRESP/*DISP*/*DIRECT* / 

C, Y, 6=0. 0/C, Y, W3=0. 0/C, Y, W4=0. 0/N0K2PP/N0M'^PP/ 
N0B2PP/MPCF1/SINGLE/0MIT/N0UE/N0K4GG/ 

N0BGG/KDEK2/-1 $ 

LBL18 $ 


LGKAD3,LGKAD $ BRANCH IF NOT FREQRESP. 

B2DD, BDD/NOBGG/ M2DD, MDD/NOSIMP/ K2DD, KDD/KDEK2 * 
LGKAD4 ^ 

LGKADo $ 

B2DD, BDD/NOGPDT /M2DD, MDD/N0SIMP/K2DD, KDD/KDEK2 * 
LGKAD4 $ 


LBLTRLl, NOTIME % 

/ /*MPY*/ REPEAT T /I /-1 $ 

//*ADD*/APPFLG /I /O INITIALIZE FOR SDRl 
TRLGLOOP ^ 


CASECC, /CASEYY/*TRAN*/S, N, REPEATT/S, N, NOLOOPl $ 
//*MPY*/NCOL /O /!•$ 


CASEYY, USETD, DLT, SLT , BGPDT , SIL, CSTM, TRL, DIT, GMD, GOD, , 
, ,PDT1,PD1, ,TOL/ NOSET/NCOL $ 

TF;L, F“DT1 ,,,,,,,, , / ,PDT, /APPFLG/*DYNAMICS« $ 

TF:L,PD1 ,,,,,,,,, / ,PD , /APPFLG/#DYNAMICG* $ 
//*ADD*/APPFLG /APPFLG /I $ APPFLG=APPFLG+ I 
TRLGD0NE,REPEATT $ 

TRLGLOOP, 100 $ 


EST,MGG/ 


ERRGR3 S 
TRLGDONE % 

TOL ,,,,,,, / PRL Z , FOL Z , REORDER 1 , RE0RDER2 
O, Y, NSEGS/V ,Y,CYC10/S,Y, LMAX = -- 1 /FKMAX / 


^ 1 1 


/ 
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123 EQUIV 
123 JUMP 
123 LABEL 
123 PRLG 


123 

COND 

123 

MPYAD 

123 

EQUIV 

123 

LABEL 

123 

COND 

123 

MPYAD 

123 

ADD 

123 

EQUIV 

123 

COND 

123 

SSG2 

123 

EQUIV 

123 

LABEL 

123 

LABEL 

123 

EQUIV 

123 

PARAML 

123 

PARAM 

123 

EQUIV 

123 

COND 

123 

PARAM 

123 

CYCTl 

123 

COND 

123 

JUMP 

123 

label 

123 

PARAM 

123 

COND 

123 

EQUIV 

123 

COND 

123 

MPYAD 

123 

label 

123 

CYCTl 

123 

COND 

123 

EQUIV 

123 

COND 

123 

MPYAD 

123 

LABEL 

123 

EQUIV 

123 

JUMP 

123 

LABEL 


123 MPYAD 
123 CYCTl 

1 23 COND 
123 EQUIV 
123 COND 
1 23 MPYAD 
123 LABEL 
123 CYCTl 

1 23 COND 
123 EQUIV 
123 LABEL 
123 COPY 
123 ADD 
1 23 PARAM 
123 LABEL 
1 23 PARAM 
123 COND 
1 23 PARAM 


LBLFRL2 $ 

LBLTRLI $ 

*DIRECT*/FREQY/IfREQ*^$° ^ Pf"f^’PSF,PDF,FOL,PHFDUM / 

PPF,PD2ER0°^/^PPFX /0^° COLUMNS IF FRLX WAS GENERATED. 

ppfx,ppf $ 

LBLFRLXI $ 

LBLFRLI , NOBASEX $ 

M2GG,BASEXG, / M2BASEXG /O $ 

pppI^ppeT*" " /< 1.0, 0.0) /(-1. 0,0.0) * 

LBLBASEl , NOSET * 

^ .POOUMl.PSPl.POP. , 

LBLBASEl $ 

LBLFRLI $ 

PPF , PDF/NOSET % 

PDF //*TRAILER*/1 /PDFCOLS % 

p'i>ZTfZcio 

LBLPDONE,CYCIO $ 

— 1 ^ NLOAD /S,N,NOGO $ 

ERRORCl.NOGO $ 

LBLPDONE t 
LBLFRL2 * 

//*N0T*/N0TCYCI0 /V,Y,CYCI0 * 

LBLTRL2,N0TCYCI0 t 
PD, PDTRZl/NOROl t 
LBLR01A,N0R01 ♦ 

PD, REORDER 1, / PDTRZl / 0 $ 

LBLROIA % 

ERR0RC1,N0G0 ♦ 

PXTRZ1,PXFZ1/N0R02 $ 

LBLR02A,N0R02 $ 

PXTRZ1,RE0RDER2, / PXFZl /O $ 

LBLR02A $ 

PXFZl, PXFl $ 

LBLTRL3 $ 

LBLTRL2 $ 

PD,RE0RDER1, / PDTR22 / 0 $ 

ERRORC 1 , NOGO % 

PXTRZ2,PXTR2/N0R02 % 

LBLR02B,N0R02 $ 

PXTRZ2, RE0RDER2, / PXTR2 /O $ 

EBLR02B $ 

flm«L's;So i 

ERRORC 1 , NOGO t 
PXFZ2,PXF1 * 

LBLTRL3 $ 

P^'pUxPf/ PXF TO COMPLEX FXF. 

//.ADD./NLOAD /FLMAx'/i i NUOAd'^^lAaX * 

LBLPDONE % - FLMAX 

//*ADD*/KMINL /V, Y, KINDEX=-1 /-i * 

N0KMINL,KMINL $ 1/ 1 $ 

//♦ADD# /KM I NV /O /Q $ 
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1 2.5 xTienc 
123 CYCT2 

123 COND 
123 ADD 
123 PRTPARM 
123 CYCT2 

123 PRTPARM 
1 23 COND 
123 PARAM 
123 REPT 
123 LABEL 
123 COND 
1 23 PRTPARM 
123 LABEL 
L23 CYCT2 

2Z COND 
23 CYCT2 

23 COND 
23 CYCT2 

23 COND 
23 READ 

23 OFP 
23 COND 
23 OFP 
23 COND 
23 CYCT2 

23 COND 
23 SDRl 
23 SDR2 

23 PLOT 

23 PRTMSG 
23 LABEL 
23 GKAM 


'RMfNCuOP '% 

CYCDD, , ,PXF, , /,,PKF2,, / *FORE#/V, Y, KSEGS/ 
KMINV/CYCSEQ/NLOAD/S,N,NOGQ $ 

ERRORCl,NOGO i 

PKFZ, / UKVFZ / (0.0, 0.0) % 

M //0/*KINDEX* $ 

CYCDD, ,, UKVFZ, , /,,UXVF,, /*BACK*/V, Y, NSEGS/ 
KMINV/CYCSEQ/NLOAD/S,N,NOGO $ ^^’Y,NSEGS/ 

M //0/itKINDEX* $ 

ERR0RC1,N0G0 % 

//*ADD*/KMINV /KMINV /I $ 

KM INLOOP, KM I NL % 

NOKMINL % 

NOKPRT,NOKPRT $ 

1 //0/*KINDEX* $ 

NOKPRT $ 

CYCDD, KDD, MDD, , , /KKKF, MKKF, , , /*FORE»/V,Y NFFPc; / 

V, Y,KINDEX/CYCSEQ/NLOAD/S,N,NOGO $ ’ ^ 

ERR0RC1,N0G0 $ 

CYCDD, BDD, , PXF, , /BKKF,,PKF,, / *FORE*/V Y MRFrc;/ 

V, Y, KINDEX/CYCSEQ/NLOAD/S, N, NOGO ^ ' 

ERR0RC1,N0G0 t 

^sess/v, y, kindex/ 

ERR0RC1,N0G0 $ 

0EI6S, , , , , //s, N, CARDNO t 
FINIS, NEIGV % 

LAMK, ,,,, //s,N, CARDNO « 
noplot , JUMPPLOT % 

C^CSEi;i;sXNOeri ’ ’ ’ Y, NSEGS/V, Y, KINDEX/ 

ERR0RC1,N0G0 $ 

CASECc’"cs?M’r:,p?’n'?;’^T®’ ’ ''PHIS- -QG/I/.REIG. $ 

,OQ0I,ipHIG:OEs;,OEFi?pPms%RE;Gfr'"’'-''"''’®''’''"'°’^ 

NOPLOT * 

wH;NouE/cr^ ; ms=o;f y°:frm ’ 
NOM2PP/NOB2PP/NOK2PP/NONCUP%?NyPMQDE”o i 
PHIKH//*TRAILER./1/S,N,NM0DES * 

PHIKH,MKKF,PHIKH, , , /MHH/3////1 t 
PHIKH,KKKF,PHIKH, , ,/KHH/3////i % 

PHIKH,BKKF,PHIKH, , , /BHH/3////I % 

PHIKH,PKF, /PHF/1 $ 

KLABEL^rKFLAG’ ^<^KF//PHF, PKF . 

i^M’l:!.?^^’®^^'^''’*^®'’'^’^Q^^^f'*’SM»eO/AERO,ACPT GTKA / 

S,N,Nk/S,N,NJ/V,Y,MINMACH/V, Y,MAXMACH/V Y IRtp//’ 

NMODES/V,Y, KINDEX % v, t , n«xmCH/V, Y, IREF// 

» DUK, D2JK/NK/NJ/ 1 $ 

NOubfi’ , .USETD, AERO/QHHL,, / 

KLABEL2 * 

KLABELI $ 

V, Y, kIndex/cycseq7^^^ ah, LAMAH/*BACK>/V, Y, NSEGS/ 

ERRORC 1 , NOGO % 

S?iL\K%:N?N?JvfY’mNSACH/v?Y;hAx'N^^^^ 

NMODES/V, Y, KINDEX $ ^^^ACH/V, Y, IREF/*COSINE*/ 

aero, ACPT/AJJL, SKJ , D UK , D2JK/NK/NJ/ 1 % 
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123' PARTN 
123 AMP 

123 APDB 


123 PARTN 
123 AMP 

123 ADD 
123 LABEL 
123 FRRD2 
123 DDRl 

123 EQUIV 
123 CYCT2 

123 COND 
123 EQUIV 
123 COND 
123 CYCTl 

1 23 LABEL 
123 COND 
123 EQUIV 
123 COND 
123 CYCTl 

1 23 LABEL 
123 SDRl 
12o> SSG2 
123 EQUIV 
1 23 LABEL 
1 23 VDR 


124 COND 

125 COND 

126 SDR3 

127 OFP 

128 XYTRAN 


129 XYPLOT 

1 30 JUMP 

131 LABEL 

132 OFP 

133 LABEL 

1 34 COND 

135 EQUIV 

136 COND 

137 SDRl 

138 LABEL 

139 SDR2 


139 CURV 

140 COND 
142 SDR3 

1 42 OFP 

143 XYTRAN 

144 XYPLOT 

145 COND 

146 RANDOM 

147 COND 


>Wlfihi;pVECT,VPHIAC, , , / I i • 

N0ue;f ^ PH I AC, , . USETD, AERO/QHHLC, , / 

NMODES/V, Y, KINDEX % ’ ^ ’ ^^^'^ACH/V, Y, IREF/*S INE*/ 

PHIAH,PVECT,/PHIAS,,,/l % 

• - USETD, AERO/QHHLS, , / 

QHHLC, QHHLS/QHHL $ 

KLABEL2 $ 

KKKF,BKKF,MKKF,QHHL,PKF,F0L/UKVF/V Y Pnu/u V 
UKVF,PHIKH/UKKVF « Y, BOV/V, Y, Q/- 1 . o $ 

UKKVFiUKVF $ 

C^CSEi;S;s,N,’NOGr;’ ^«ACK,/V,V,MSEGS/D,V,KINDEX/ 
errorci,nogo $ 

UXVF,UDVF / CYCIO $ 

UXVF /^UDVF,6CYcL^/^CTYPE/^ TRANSFORM TO PHYSICAL. 

NLOAD $ Y, NSEGS/V, Y, KMAX/ 

LCYC3 $ 

LBLTRL4, NOTIME * 

PXF,PDF2 / CYCIO $ 

Pxr^TOP2?eCYCB2 /'"c^^PE^BACkI/V^TeapT ™ 

NLOAD * Y,NSEGS/V, Y,KMAX/ 

LCYC4 $ 

LBLTRL4 t 

LBL15A,NOSORT2 $ 

OUDVCl, , , , ,/0UDVC2, , , , , * 

0UDVC2, , , , , //s,N,CARDNO $ 

SyN°§A°DNO^^ ’ ’ ’ ^>'"^^"^«^«^-^E«*/*DSET,/S,N,PriLE/ 

XYPLTFA// $ 

LBL15 % 

LBL15A $ 

OUDVCl, , , , ,// s,n,CARDNO $ 

LBL15 $ 

LBL20,N0P $ 

UDVF, UPVC/NOA % 

LBL19,N0A $ 

USETD, ,UDVP, - .SOD,eHD,PSF,KPS, ,/UPUC, ,QPC/1/.dyNAHICS. , 

PPF/OPPCl JoQPCU QUPVCl °QESCl''oEFcf'^Pnr’ QPC, UPVC, EST, XYCDB, 
s, N, N0S0RT2 $ ’ OESCl , OEFCl , PUPVC 1 /*FREQRESP*/ 

LBL?7;N0S0R?r;^^^ " 

0ESC2, OEFC2,’oESC2m’$^^^^ ’ O^FCl , 0ESC1M/0PPC2, 0QPC2, 0UPVC2, 

xycdb;opp§2;oqSpvcT§^^^^ * 

s,n,pfile/s,n,cardno $ ^ "’^^‘^^‘^^/^ypltf/^freq^/^pset*/ 

XYPLTF// $ • 

LBL16,N0PSDL $ 

X YCDB, D I T, PSDL 1 0UPVC2. OPPC'^ nppr'^ 

SjNtNORD * "■’°^^^'^’°^-^^^’0EFC2,CASEXX/PSDF,AUT0/ 

LBL16,NQRD 4; 
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I 


148 

XYTRAN 

XYCDB,PSDF; AUTO, , , 
S,N,CARDNO $ 

149 

XYPLOT 

XYPLTR// % 

150 

JUMP 

LBL16 $ 

151 

label 

LBL17 It 

152 

PURGE 

PSDF/N0S0RT2 $ 

153 

OFP 

OUPVC1,OPPC1,OQPC1 

154 

LABEL 

LBL16 $ 

155 

PURGE 

PSDF/NOPSDL * 

1 56 

COND 

LBL20, JUMPPLOT <t 

157 

PLOT 

PLTPAR, GPSETS, ELSE' 
GPECT , OESCl /PL0TX2. 


PRTMSG 

LABEL 

COND 

REPT 

LABEL 

PRTPARM 

JUMP 

LABEL 

PRTPARM 

LABEL 

PRTPARM 

LABEL 

PRTPARM 

JUMP 

LABEL 

PRTPARM 

LABEL 

PRTPARM 

JUMP 

LABEL 

PRTPARM 

PRTPARM 

JUMP 

LABEL 

PRTPARM 

PRTPARM 

JUMP 

LABEL 

PRTPARM 

LABEL 


PLOTX2// $ 

LBL20 % 

FINIS, REPEATF $ 

LBL13, 100 It 
ERR0R3 % 

//-3/*DIRFRRD* $ 

FINIS $ 

ERR0R4 % 

//-4/*DIRFRRD# $ 

//-r%c^cs“ncsff D«Tft. 

ERR0RC2 $ COUPLED MASS NOT ALLOWED 
//O /C, YjCOUPMASS $ 

FINIS $ 

/%°/SnqSe.T"''' allowed. 

FINIS $ 

//0°/“oFRLf IN BULK DATA DECK. 

//O /*NOTRL» $ 

FINIS $ 

ERR0RC6 $ BOTH FREQ AND TSTFP updct cm 

//O /*NOFREQ* $ WERE SELECTED IN CASE CONTROL. 

//O /*NOTIME* Si 
FINIS $ 

//-2/^cycwoDEs“f extraction data 

FINIS $ 


171 PURGE DUMMY/ALWAYS % 

172 END $ 
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SECTION 4 


SUPPLEMENT TO 

NASTRAN DEMONSTRATION MANUAL 
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MODAL FORCED VIBRATION ANALYSIS OF 
AERODYNAHICALLY EXCITED TURBOSYSTEMS 

4.1 INTRODUCTION 

The principal purpose of this section is to demonstrate the 
use of the newly developed capability in NASTRAN to conduct modal 
forced vibration analysis of rotating turbosystems subjected to 
excitation from aerodynamic sources. The demonstration comprises 
a series of four inter-related analysis phases: 

Phase 1 generates a total stiffness matrix consisting of 
elastic plus differential stiffness matrices. 

Phase 2 ascertains the aeroelastlc stability of the 
turbosystem before proceeding with response analysis. 

Phase 3 generates the applied oscillatory airloads on the 
blades of the turbosystein. 

Phase 4 determines the aerodynamically forced response of the 
turbosystein, 

^•2 EXAMPLE PROBLEM DESCRIPTION 

An elght-bladed single-rotation advanced turboprop (Figure 
4.1) is selected as an example of turbosystems. 

The swept blades of the turboprop are set at an angle of 60.8° 
with the plane of rotation, when measured at 3/4 tip radius. The 
prop rotates at a constant 8000 rpm. its axis of rotation is 
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inclined at 2 with the uniform absolute inflow. The freestream 
inflow conditions are given by 0.798 Mach number, 873 fps inflow 
velocity, and 1.9034 x lo'^ Ibf-secVft'^ inflow density. 

This operating condition results in oscillatory airloads 
acting on the blades of the turboprop at an excitation frequency 
of one-per-rev (corresponding to 133.34 Hz.). 

It IS desired to obtain the resultant blade surface vibratory 
stress distribution for comparison with experimental 
Observations • 

4.3 INPUT 

Figure 4.2 illustrates the NASTRAN model of one representative 
blade of the turboprop. For present, the hub is considered 

rigid, and the blade is completely fixed at the bottom of its 
shank. 

The input data decks for all four phases discussed above are 
included in this section. 

Phases 1, 2, and 4 use NASTRAN, while Phase 3 uses AIRLOADS 
program (Ref. 2). 

4.4 RESULTS 

Figure 4.3 compares the analytical and test stresses. 
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ORIGir-JAL PAGE IS 

OF POOR QUALITY 







Figure 4.2 NASTRAN Structural and 

Aerodynamic Models of SR-3 
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MAGNITUDE OF OSCILLATORY STRESS (PSI) 


COMPARISON OF ANALYTICAL AND TEST STRESSES 
SR3; NASA TEST READING NO. 273 


X TEST VALUE 



RANGE COVERED BY FINITE 
ELEMENTS NEIGHBORING THE 
STRAIN GAGE LOCATION 



GAGE 1 gage 2 GAGE 3 GAGE 4 GAGE 5 


(SHE>W) 

STRAIN GAGES 

Figure 4.3 SR-3 One-Per-Rev Stress 
Comparison — Test Readina 
No. 273 
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INPUT DATA DECK FOR 


DIFFERENTIAL STIFFNESS ANALYSTS 


NASTRAN BANDIT = -1, FILES = (INPT, PLT2) 


EXECUTIVE CONTROL DECK 


ID 

NASA,SR3PR0P 

APP 

DISP 

SOL 

4 

DIAG 

8, 14,21,22 

TIME 

10 It CRAY-1 


$ 


i ALTERS TO SAVE ELASTIC PLUS DIFFERENTIAL STIFFNESS 

* (APRIL 1984 VERSION) 

5 

ALTER 149 ^ 

ADD DKDG6,KDGG / KDGGX / (-1.0, 0.0) $ 

ADD KGG, KDGGX / KTOTAL $> 

OUTPUT 1 KTOTAL, ,, ,//-l/0 * 

OUTPUT 1, ,,,,//-3/0 ♦ 

ENDALTER * 

$ 

CEND 



(KTOTAL) 
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CASE CONTROL DECK 


% 

TITLE = SR3 RESPONSE TO 1 PER REV DSC. AIRLOADS 
SUBTITLE = NASA TEST READING NO. 273 
LABEL = DIFFERENTIAL STIFFNESS ANALYSIS 
% 

SPC = 1 

MPC = 1 

LOAD = 1 


SUBCASE 1 

LABEL = DIFF. STIFF. ANAL, — LINEAR SOLUTION 
DISP = ALL 
STRESS = ALL 
SUBCASE 2 

LABEL = DIFF. STIFF. ANAL. — NONLINEAR SOLUTION 
DISP(SORTl, PRINT) - ALL 
STRESS = ALL 

% 

OUTPUT <PLOT) 

SET 1 = all 

PLOTTER NASTPLT, MODEL D, 0 
PAPER SIZE 8.0 BY 8.0 
MAXIMUM DEFORMATION 1.0 
FIND SCALE, ORIGIN 1, SET 1 
PTITLE = SOL 4 
CONTOUR YDISPLAC 

PLOT STATIC DEFORMATION CONTOUR O, 1 , SET i, ORIGIN 1 
CONTOUR YDISPLAC i, 

RPrTM'"DM, DEFORMATION CONTOUR 0, 2, SET 1, ORIGIN 1, 

BEGIN BULK 


OUTLINE 

OUTLINE 
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BULK DATA DECK 


RFORCE 1 

0 

0 

133.33 1 

. 0 

MATl 

1 1 , 

.6 E7 

.35 

0004141 

CTRIA2 

1 

1 10 

9 

8 

PTRIA2 

1 

1 0.01570 


CTRIA2 

2 

2 1 1 

10 

8 

PTRIA2 

2 

1 0,02827 


CTRIA2 

4 

4 12 

1 1 

7 

PTRIA2 

4 

1 0.03380 



CTRIA2 

6 

6 13 

12 

1 

PTRIA2 

6 

1 0.03623 


CTRIA2 

7 

7 1 

2 

13 

PTRIA2 

7 

1 0.02917 


CTRIA2 

8 

8 14 

13 

2 

PTRIA2 

8 

1 0.04440 



CTRIA2 

9 

9 2 

3 

14 

PTRIA2 

9 

1 0.03830 


CTRIA2 

10 

10 15 

14 

3 

PTRIA2 

10 

1 0.04253 


CTRIA2 

11 

1 1 3 

4 

15 

PTRIA2 

11 

1 0.03677 


CTRIA2 

12 

12 16 

15 

4 

PTRIA2 

12 

1 0.03397 



CTRIA2 

13 

13 4 

5 

16 

PTRIA2 

13 

1 0.02740 


CTRIA2 

14 

14 17 

16 

5 

PTRIA2 

14 

1 0.01673 



CTRIA2 

15 

15 5 

6 

17 

PTRIA2 

15 . 

1 0.00823 


CTRIA2 

16 

16 20 

19 

18 

PTRIA2 

16 

1 0.01970 


CTRIA2 

17 

17 21 

20 

18 

PTRIA2 

17 

. 1 0.03550 

CTRIA2 

19 

19 22 

21 

9 

PTRIA2 

19 

1 0.04250 



CTRIA2 

20 

20 9 

10 

22 

PTRIA2 

20 

1 0.03487 


CTRIA2 

21 

21 10 

11 

22 

PTR1A2 

21 

1 0.04743 


CTRIA2 

22 

22 23 

22 

11 

PTRIA2 

22 

1 0.05847 


CTRIA2 

23 

23 11 

12 

23 

PTRIA2 

23 

1 0.05413 


CTRIA2 

24 

24 24 

23 

12 

PTRIA2 

24 

1 0.06033 



CTRIA2 

25 

25 1 2 

13 

24 

PTRIA2 

25 

1 0.05580 



CTRIA2 

26 

26 25 

24 

13 

PTRIA2 

26 

1 0.05663 



CTRIA2 

27 

27 13 

14 


PTRIA2 

27 

1 0.05230 

CTRIA2 

28 

28 14 

15 


PTRIA2 

28 

1 0.04877 




CTRIA2 AND PTRIA2 DATA IDENTICAL TO THAT FOR FORCED RESPONSE ANALYSIS 


ORIGfSIAL PAGE IS 
OF POOR QUALITY 
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CTRIA2 

340 

340 194 

195 

202 

PTRIA2 

340 

1 1 . 06600 



CTRIA2 

341 

341 195 

196 

202 

PTRIA2 

341 

1 1 . 06600 



CTRIA2 

342 

342 203 

202 

196 

PTRIA2 

342 

1 1.13913 



CTRIA2 

343 

343 196 

197 

203 

PTRIA2 

343 

1 1 . 06487 



CTRIA2 

344 

344 204 

203 

197 

PTRIA2 

344 

1 1 . 06620 



CTRIA2 

345 

345 197 

198 

204 

PTRIA2 

345 

1 0.70640 



CTRIA2 

346 

346 205 

204 

198 

PTRIA2 

346 

1 0.55867 



CTRIA2 

3 

3 8 

7 

11 

PTRIA2 

3 

1 0.01897 


CTRIA2 

5 

5 7 

1 

1 2 

PTRIA2 

5 

1 0.02043 


1 4CL 

CTRIA2 

18 

18 18 

9 

21 

PTRIA2 

18 

1 0.02390 



CTRIA2 

143 

143 79 

83 

92 

PTRIA2 

143 

1 0.13113 



CTRIA2 

166 

166 92 

96 

105 

PTRIA2 

166 

1 0. 14623 



CTRIA2 

189 

189 105 

109 

118 

PTRIA2 

189 

1 0.15213 



CTRIA2 

212 

212 118 

122 

131 

PTRIA2 

212 

1 0. 15470 



CTRIA2 

235 

235 131 

135 

144 

PTRIA2 

235 

1 0. 15577 



CTRIA2 

275 

275 166 

165 

153 

PTRIA2 

275 

1 0.51933 



CTRIA2 

277 

277 167 

166 

154 

PTRIA2 

277 

1 0.31070 



CTRIA2 

285 

285 176 

177 

160 

PTRIA2 

285 

1 0.67620 



CTRIA2 

287 

287 168 

177 

160 

PTRIA2 

287 

1 0.72723 



CTRIA2 

288 

288 168 

178 

177 

PTRIA2 

288 

1 0.65833 



C0RD2R 77 

0 

.0 .0 

. 0 


+C2R 10, 

-0. 

618 .0 




o 


GRDSET 77 

GRID 1 

GRID 3 

GRID 4 

GRID 5 

GRID 6 

GRID 7 

GRID 8 

GRID 9 

GRID 10 

GRID 12 


1.808 

1.839 

12.250 

2.376 

2.347 

12.250 

2.625 

2. 558 

12.250 

2.877 

2.765 

12.250 

3. 134 

2,966 

12.250 

1.556 

1.589 

12.033 

1.304 

1 . 339 

11.817 

1.052 

1.088 

1 1 . 600 

1.293 

1.308 

1 1 . 600 

1.791 

1 . 730 

1 1 . 600 


.0 


+C2R 


GRID DATA IDENTICAL TO THAT FOR FORCED RESPONSE ANALYSIS 


OJRIGIMAL PAGE IS 
OF POOR QUALiTY 
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» GRID 
' GRID 

• GRID 

• GRID 
' GRID 
■ GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

grid 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

grid 

GRID 

GRID 

GRID 

grid 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

grid 

GRID 

GRID 

grid 

grid 

GRID 

GRID 

MPC 1 

MPC 1 

HPC 1 

SEQGP 

SEQ6P 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 


.170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 
181 
182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 
2 

1 1 
21 
30 
39 
48 
57 
66 
75 
88 
101 
114 
127 
140 
150 
160 

6 

5 

7 

1 

9 

13 

17 

21 

29 

37 






OF 

p0®e QUALriT 



-0. 183 

-0. 030 

V3U 


0. 000 

0. 000 

2. 930 


0. 183 

0. 030 

2.930 


0.365 

0,059 

2. 930 


0.548 

0.089 

2. 930 


-1.804 

-0. 270 

2.650 


-1. 188 

-0. 182 

2. 650 


-0.750 -0.123 2.740 

-0.550 -0.072 2.600 

-0.367 -0.048 2.600 


4 

4 

4 

9 


20 

16 

3 

28 

ti-X. 

38 

34 

47 


-O. 184 
0 . 000 
O. 184 
0.367 
0.550 
-0.550 
-0.367 
-0. 184 
0.000 
O. 184 
0.367 
0.550 
-0. 550 
-0.367 
-O. 184 
0.000 
0. 184 
0.367 
0.550 
-0.699 
-0.466 
-0.233 
0.000 
0.233 
0.466 
0.699 
-1.496 
2. 129 
1.541 
1.034 
0.550 
0 . 101 
-0.317 
-0.677 
-O. 998 
-1.271 
-1.490 
-1.621 
-1.642 
-1.572 
-1.435 
-1.208 
-0.917 

1 

1 . 

1 . 

6 

10 

14 

18 

26 

30 

34 

38 


-0. 024 
0.000 
0.024 
0.048 
0.072 
-0.072 
-O. 048 
-0.024 
0 . 000 
0.024 
O. 048 
0.072 
-0.072 
-O. 048 
-0.024 
0.000 
0.024 
0.048 
0.072 
-0.091 
-0.061 
-0.030 
0.000 
0. 030 
0. 061 
0.091 
-0.226 
2. 133 
1.520 
1. 014 
O. 564 
O. 182 
-O. 138 
-0.383 
-O. 572 
-O. 702 
-O. 779 
-0. 796 
-0.753 
-O. 662 
-0.538 
-0.374 
-0. 229 
0 4 

O 4 

0 1 

8 
1 

19 

IS 

21 

27 

13 

37 

31 

48 


2. 600 
2.600 
2.600 
2.600 
2.600 
2.350 
2.350 

2 . 350 

2.350 
2. 350 

2.350 
2. 350 
2. 070 
2.070 
2.070 
2.070 
2.070 
2.070 
2, 070 
1.920 
1.920 
1 . 920 
1.920 
1.920 
1.920 
1 . 920 
2.650 

12.250 
1 1.600 
1 1 . 000 

10. 400 
9.800 
9. 187 
8 . 600 
8.000 
7. 400 
6 . 800 
6 . 200 
5. 600 
5. 000 

4 . 400 
3.715 
3. 180 

4 
4 
4 

7 
1 1 
15 
19 

27 
31 

39 


- 1.0 

- 1.0 

- 1.0 

7 

10 

18 

12 

30 

26 

6 

36 

46 


1 


28 

36 

40 


39 


14 

45 


91 


-ti O O M m 



SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SEQGP 

SPCl 1 

SPCl 1 

SPCl 1 

SPCl 1 

SPCl 1 

SPCl 1 

ENDDATA 


41 

44 

45 

24 

49 

54 

53 

41 

57 

65 

61 

58 

65 

77 

69 

72 

73 

86 

77 

83 

81 

61 

85 

70 

89 

97 

93 

92 

97 

94 

101 

1 10 

105 

104 

109 

1 15 

113 

101 

117 

1 19 

121 

132 

125 

102 

129 

133 

133 

150 

137 

169 

141 

146 

145 

154 

149 

136 

153 

153 

157 

125 

161 

160 

165 

165 

169 

172 

173 

176 

177 

159 

181 

183 

185 

181 

189 

192 

193 

195 

197 

200 

201 

202 

205 

206 

cr 

»_j 

17 

5 

23 

vJ 

186 

5 

191 

5 

198 

1234 

56 199 


42 

43 

46 

56 

50 

53 

54 

33 

58 

64 

62 

50 

66 

75 

70 

68 

74 

87 

78 

82 

82 

60 

86 

88 

90 

96 

94 

81 

98 

93 

102 

109 

106 

114 

1 10 

118 

114 

111 

118 

116 

122 

130 

126 

1 12 

130 

129 

134 

157 

138 

113 

142 

140 

146 

167 

150 

147 

154 

155 

158 

137 

162 

161 

166 

164 

170 

173 

174 

177 

178 

171 

182 

184 

186 

187 

190 

193 

194 

197 

198 

201 

202 

203 

206 

138 

27 

36 

54 

37 

192 

187 

193 

194 

THRU 

205 


43 

40 

47 

57 

51 

52 

55 

66 

59 

63 

63 

42 

67 

74 

71 

59 

75 

85 

79 

78 

83 

79 

87 

99 

91 

95 

95 

80 

99 

89 

103 

108 

107 

106 

111 

117 

1 15 

121 

1 19 

127 

123 

139 

127 

122 

131 

131 

135 

144 

139 

123 

143 

141 

147 

179 

151 

151 

155 

166 

159 

148 

163 

162 

167 

168 

171 

174 

175 

126 

179 

180 

183 

185 

187 

190 

191 

194 

195 

198 

199 

189 

203 

204 

14 

10 

112 

146 

188 

189 

195 

196 


44 

32 

48 

55 

52 

49 

56 

67 

60 

62 

64 

76 

68 

73 

72 

51 

76 

84 

80 

69 

84 

71 

88 

98 

92 

91 

96 

103 

100 

100 

104 

107 

108 

105 

112 

90 

1 16 

120 

120 

128 

124 

145 

128 

134 

132 

142 

136 

156 

140 

135 

144 

143 

148 

124 

152 

152 

156 

178 

160 

158 

164 

163 

168 

170 

172 

175 

176 

149 

180 

182 

184 

186 

188 

191 

192 

188 

196 

199 

200 

196 

204 

205 

45 


185 


190 


197 



ORIGINAL PAGE IS 
OF POOR QUALITY 
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INPUT DATA DECK FOR 


aeroelastic stabi lity analysts 

NASTRAN BANDIT *= -1, SYSTEM(93) =1 , FILES = (INPT, PLT2) 


Mecutive control DEOR 


ID 

NASA,SR3PR0P 

APP 

AERO 

SOL 

9 

DIAG 

8, 14,21,22 

TIME 

10 $ CRAY-1 




* ALTERS TO ADD - 0MEGA**2 M TERMS TO TOTAL STIFFNESS 
% 

ALTER 26 

//*MPY*/ OMEGA /V, Y,RPS«=0.0/ 6.283185 * 

//*MPY*/ OMEGASQR / OMEGA /OMEGA ♦ 

//*SUB*/ MOMEGASQ /O.O /OMEGASQR % ' 

//^COMPLEX*// MOMEGASQ /O.O/ CMPLX2 $ 

CASECC, BGPDT, CSTM, DIT , « MGG, y /, DUMDUM, MIGG, , , , , /NOMGG/V, Y, CYCIO/ 

V, Y,NSEGS/0/S,N,FKMAX/-l/-i/-l/-l/-l/-l/s,N,NOBASEX/lVOMEGA * 

KGGX,M1GG / KGGl / (1.0, O.O) / CMPLX2 * 

KGG1,KGGX « 

ENDALTER 

♦ 

CEND 


PARAMR 

PARAMR 

PARAMR 

PARAMR 

% 

FVRSTRl 

« 

ADD 

EQUIV 
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CASE CONTROL DECK 


response TO 1 PER REV OSC. 

reading no. 273 

LABEL = AEROELASTIC STABILITY ANALYSIS 


AIRLOADS 


SPC = 1 

MPC = 1 

METHOD = 1 
FMETHOD = 1 


OUTPUT (XYOUT) 

PLOTTER NASTPLT, MODEL D, O 
X PAPER =8.0 
YPAPER = 11.0 
YAXIS = YES 

XTaJis'^Hes * operating (MEAN RELATIVE) VELOCITY 

XBAXIS = YES 

CURVELINESYMBOL = 5 

^DIVISIONS = 10 

YTDIVISIONS = 10 

YBDIVISIONS = 10 

YTGRID LINES = YES 

YBGRID LINES = YES 

XTGRID LINES = YES 

XBGRID LINES = YES 


XTITLE = 
YTTITLE = 
YBTITLE = 
TCURVE = K 
XYPLOT, XYPRINT 
BEGIN BULK 


VELOCITY VSBAR IN/SEC 
DAMPING G 
frequency F, H2 
.179 SIGMA = “45.0 DEG. 

VG/ 1 <G,F) ,2(G,F> ,3(G,F) ,4(G,F) ,5(G,F) 
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BULK DATA DECK 


C0RD2R 

77 

0 

. 0 

+C2R 

10. 

-0. 

618 .0 

CTRIA2 

1 

1 

10 

CTRIA2 

2 

2 

1 1 

CTRIA2 

3 

3 

8 

CTRIA2 

4 

4 

12 

CTRIA2 

5 

5 

7 

CTRIA2 

6 

6 

13 

CTRIA2 

7 

7 

1 

CTRIA2 

8 

8 

14 

CTRIA2 

9 

9 

2 

CTRIA2 

10 

10 

15 

CTRIA2 

11 

11 

3 

CTRIA2 

12 

12 

16 

CTRIA2 

13 

13 

4 

CTRIA2 

14 

14 

17 

CTRIA2 

15 

15 

CT 

O 

CTRIA2 

16 

16 

20 

CTRIA2 

17 

17 

21 

CTRIA2 

18 

18 

18 

CTRIA2 

19 

19 


CTRIA2 

20 

20 

9 

CTRIA2 

21 

21 

10 

CTRIA2 

22 

22 

23 

CTRIA2 

23 

23 

1 1 

CTRIA2 

24 

24 

24 

CTRIA2 

25 

25 

12 

CTRIA2 

26 

26 

25 

CTRIA2 

27 

27 

13 

CTRIA2 

28 

28 

14 

CTRIA2 

29 

29 

26 

CTRIA2 

30 

30 

15 

CTRIA2 

31 

31 

27 

CTRIA2 

32 

32 

16 

CTRIA2 

33 

33 

29 

CTRIA2 

34 

34 

30 

CTRIA2 

35 

35 

19 

CTRIA2 

36 

36 

31 

CTRIA2 

37 

37 

20 

CTRIA2 

38 

38 

32 

CTRIA2 

39 

39 

21 

CTRIA2 

40 

40 

33 

CTRIA2 

41 

41 

22 

CTRIA2 

42 

42 

23 

CTRIA2 

43 

43 

34 

CTRIA2 

44 

44 

24 

CTRIA2 

45 

45 

35 

CTRIA2 

46 

46 

25 

CTRIA2 

47 

47 

36 

CTRIA2 

48 

48 

26 

CTRIA2 

49 

49 

38 

CTRIA2 

50 

50 

28 

CTRIA2 

51 

51 

39 


. 0 

9 

10 
7 

1 1 
1 

12 

o 

4C. 

13 

3 

14 

4 

15 

vJ 

16 
6 

19 

20 

9 
21 

10 
11 
22 
12 

23 

13 

24 

14 

15 

25 

16 

26 
17 
28 

29 
20 

30 
21 

31 
22 

32 

23 

24 


25 

34 

26 

35 
27 

37 
29 

38 


.0 

8 
8 
1 1 
7 
12 
1 

13 
2 

14 

3 

15 

4 

16 

5 

17 

18 
18 
21 
9 


11 

23 
12 

24 
13 

25 

25 

15 

26 

16 

27 
19 

19 

30 

20 

31 
21 

32 
22 

33 

33 

24 

34 

25 

35 

26 

36 

28 
39 
28 


. 0 


. O 


+C2R 


CTR1A2 DATA IDENTICAL TO THAT FOR FORCED RESPONSE ANALYSIS 
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1 V4 

IV^ 187 

CTRIA2 

327 

327 

187 

188 194 

CTRIA2 

328 

328 

195 

194 188 

CTRIA2 

329 

329 

188 

189 195 

CTRIA2 

330 

330 

196 

195 189 

CTRIA2 

331 

331 

189 

190 196 

CTRIA2 

332 

332 

197 

196 190 

' CTRIA2 

333 

333 

190 

191 197 

CTRIA2 

334 

334 

198 

197 191 

' CTRIA2 

335 

335 

200 

199 192 

CTRIA2 

336 

336 

192 

193 200 

CTRIA2 

337 

337 

201 

200 193 

CTRIA2 

338 

338 

193 

194 201 

CTRIA2 

339 

339 

202 

201 194 

CTRIA2 

340 

340 

194 

195 202 

CTRIA2 

341 

341 

195 

196 202 

CTRIA2 

342 

342 

203 

202 1 96 

CTRIA2 

343 

343 

196 

197 203 

CTRIA2 

344 

344 

204 

203 197 

CTRIA2 

345 

345 

197 

198 204 

CTRIA2 

346 

346 

205 

204 198 

CYJOIN 

1 


199 


CYJOIN 

2 


205 


EIGR 

1 

PEER 



+EIG1 

MAX 




GRDSET 


77 



GRID 


1 

1.808 

1.839 12.250 

GRID 


3 

2.376 

2.347 12.250 

GRID 


4 

2.625 

2.558 12.250 

GRID 


5 

2.877 

2.765 12.250 

GRID 


6 

3- 134 

2.966 12.250 

GRID 


7 

1.556 

1.589 12.033 

GRID 


8 

1.304 

1.339 11.817 

GRID 


9 

1.052 

1 . 088 1 1 . 600 

GRID 


10 

1.293 

1 . 308 1 1 . 600 

GRID 


12 

1.791 

1 . 730 1 1 . 600 

GRID 


13 

2.044 

1 . 937 1 1 . 600 

GRID 


14 

2.298 

2.141 11.600 

GRID 


15 

2.555 

2.343 11.600 

GRID 


16 

2.814 

2,542 11.600 

GR I D 


17 

3.078 

2.733 11.600 

GR I D 


18 

0.737 

0.793 11.300 

GR I D 


19 

0.423 

0.498 11.000 

GRID 


20 

0. 725 

0.761 11.000 

GRID 


22 

1.347 

1 . 264 1 1 . 000 

GRID 


23 

1.663 

1.510 11.000 

GRID 


24 

1.981 

1 . 753 1 1 . 000 

GRID 


25 

2.302 

1 . 993 1 1 . 000 

GRID 


26 

2.626 

2.228 11.000 

GRID 


27 

2.956 

2.454 11.000 

GRID 


28 

-0. 168 

-0.013 10.400 

GRID 


29 

0. 186 

0.281 10.400 

GR I D 


31 

0.917 

0.842 10.400 

GRID 


32 

1.288 

1.115 10.400 


GRID DATA IDENTICAL TO THAT FOR FORCED RESPONSE ANALYSIS 


ORIGINAL PAGE IS 
OF POOR QUALITY 


+EI61 
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GRID 


188 

0. 000 

0.000 

2.350 


GRID 


189 

0. 184 

0. 024 

2.350 


GRID 


190 

0. 367 

0. 048 

2.350 


GRID 


l^L 

0. 550 

0.072 

2.350 


GRID 


192 

-0. 550 

-0.072 

2.070 


GRID 


193 

-0.367 

-0.048 

2.070 


GRID 


194 

-0. 184 

-0.024 

2.070 


GRID 


195 

0.000 

0.000 

2.070 


GRID 


196 

0. 184 

0.024 

2.070 


GRID 


197 

0.367 

0.048 

2.070 


GRID 


198 

0.550 

0.072 

2.070 


GRID 


199 

-0.699 

-0.091 

1.920 


GRID 


200 

-0.466 

-0.061 

1.920 


GRID 


201 

-0.233 

-0.030 

1.920 


GRID 


202 

0.000 

0.000 

1.920 


GRID 


203 

0.233 

0.030 

1.920 


GRID 


204 

0.466 

0.061 

1.920 


GRID 


205 

0.699 

0.091 

1.920 


GRID 


206 

-1.496 

-0.226 

2.650 


GRID 


2 

2. 129 

2. 133 

12.250 


GRID 


11 

1.541 

1.520 

1 1 . 600 


GRID 


21 

1 . 034 

1.014 

1 1 . 000 


GRID 


30 

0.550 

0.564 

10.400 


GRID 


39 

0. 101 

0. 182 

9,800 


GRID 


48 

-0.317 

-0. 138 

9. 187 


GR I D 


57 

-0.677 

-0.383 

8.600 


GRID 


66 

-0.998 

-0.572 

8.000 


GRID 


75 

-1.271 

-0.702 

7. ‘400 


GRID 


88 

-1.490 

-0.779 

6.800 


GRID 


101 

-1.621 

-0.796 

6.200 


GRID 


114 

-1.642 

-0.753 

S.600 


GRID 


127 

-1.572 

-0.662 

5.000 


GRI D 


140 

-1.435 

-0.538 

4.400 


GRID 


150 

-1.208 

-0.374 

3.715 


GRID 


160 

-0.917 

-0.229 

3. 180 


MATl 

1 

1.6 

E7 

35 

0004141 


MPC 

1 

5 

4 1 

.0 4 

4 

-1.0 

MPC 

1 

6 

4 1 

.0 4 

4 

— 1.0 

MPC 

1 

7 

4 1 

. 0 1 

4 

-1.0 

PTRIA2 

1 

1 

.01570 




PTRIA2 

2 

1 

. 02827 




PTRIA2 

3 

1 

.01897 




PTRIA2 

4 

1 

. 03380 




PTRIA2 

5 

1 

. 02043 




PTRIA2 

a 

1 

. 03623 




PTRIA2 

7 

1 

. 029 1 7 




PTRIA2 

Q 

1 

. 04440 




PTRIA2 

9 

1 

. 03830 




PTRIA2 

10 

1 

. 04253 




PTRIA2 

11 

1 

. 03677 




PTRIA2 

12 

1 

. 03397 




PTRIA2 

13 

1 

. 02740 




PTRIA2 

14 

1 

.01673 





PTRIA2 DATA IDENTICAL TO THAT FOR FORCED RESPONSE ANALYSIS 


OF POOR 




QUALITY 
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PTRIA2 

289 

1 

. 76263 

PTRIA2 

290 

1 

.71047 

PTRIA2 

291 

1 

. B0953 

PTRIA2 

292 

1 

. 76303 

PTRIA2 

293 

1 

. 85827 

PTRIA2 

294 

1 

,90487 

PTRIA2 

295 

1 

.83280 

PTRIA2 

296 

1 

. 86773 

PTRIA2 

297 

1 

. 76483 

PTRIA2 

298 

1 

.60217 

PTRIA2 

299 

1 

. 83203 

PTRIA2 

300 

1 

.64773 

PTRIA2 

301 

1 

. 93570 

PTRIA2 

302 

1 

. 90603 

PTRIA2 

303 

1 

. 99047 

PTRIA2 

304 

1 

1.02140 

PTRIA2 

305 

1 

1.01607 

PTRIA2 

306 

1 

1 . 03707 

PTRIA2 

307 

1 

. 89447 

PTRIA2 

308 

1 

. 97073 

PTRIA2 

309 

1 

. 43940 

PTRIA2 

310 

1 

, 62880 

PTRIA2 

311 

1 

. 45607 

PTRIA2 

312 

1 

. 62880 

PTRIA2 

313 

1 

. 89447 

PTRIA2 

314 

1 

. 97073 

PTRIA2 

315 

1 

1 . 06800 

PTRIA2 

316 

1 

1.08900 

PTRIA2 

317 

1 

1 . 08900 

PTRIA2 

318 

1 

1 . 06800 

PTRIA2 

319 

1 

. 97073 

PTRIA2 

320 

1 

. 89447 

PTRIA2 

321 

1 

. 62880 

PTRIA2 

322 

1 

. 43940 

i PTRIA2 

323 

1 

. 43940 

|PTRIA2 

324 

1 

. 62880 

J PTRIA2 

325 

1 

. 89447 

PTRIA2 

326 

1 

. 97073 

PTRIA2 

327 

1 

1 . 06800 

PTRIA2 

328 

1 

1 . 08900 

PTRIA2 

329 

1 

1 . 08900 

PTRIA2 

330 

1 

1 . 06800 

PTRIA2 

331 

1 

. 97073 

PTRIA2 

332 

1 

. 89447 

PTRIA2 

333 

1 

. 62880 

PTRIA2 

334 

1 

. 43940 

PTRIA2 

335 

1 

. 55867 

PTRIA2 

336 

1 

. 70640 

PTRIA2 

337 

1 

1 . 06620 

PTRIA2 

338 

1 

1 . 06487 

PTRIA2 

339 

1 

1. 13913 

PTRIA2 

340 

1 

1 . 06600 

PTRIA2 

341 

1 

1 . 06600 

PTRIA2 

342 

1 

1. 13913 

PTRIA2 

343 

1 

1 . 06487 

PTRIA2 

344 

1 

1.06620 

PTRIA2 

345 

1 

.70640 

PTRIA2 

346 

1 

.55867 

SEQGP 

1 

9 

2 

SEQGP 

5 

<<cl 

6 

SEQGP 

9 

20 

10 

SEQGP 

13 

1 6 

14 

SEQGP 

17 

X 

18 

SEQGP 

21 

28 


SEQGP 

25 


26 

SEQGP 

29 

38 

30 


OF POOR QUALITY 


8 

3 

7 

4 

1 

7 

10 

8 

19 

1 1 

18 

12 

15 

15 

12 

16 

21 

19 

30 

20 

27 

23 

26 

24 

13 

27 

6 

28 

37 

31 

36 

32 
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SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
' SEQGP 
■ SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
SEQGP 
' SEQGP 
SEQGP 
' SEQGP 
SEQGP 
■ SPCl 
SPCl 
SPCl 
SPCl 
SPCl 
' SPCl 
PARAM 
PARAM 
AERO* 

' *AERO 
PARAM 
PARAM 
PARAM 
PARAM 
PARAM 
■ STREAMLl 
’ STREAML2 
+2 10 

STREAMLl 
STREAML2 
+2 :..'0 


37 
41 
45 
49 
53 
57 
61 
65 
69 
73 
77 
81 
85 
89 
93 
97 
101 
105 
109 
113 
117 
121 
1 "^5 
129 
133 
137 
141 
145 
149 
153 
157 
161 
165 
169 
173 
177 
181 
185 
189 
193 
197 
201 
205 
1 
1 
1 
1 
1 
1 

LMODES 
KGG I N 
O 


I REF 
MAX MACH 
MINMACH 
NSEGS 
RPS 

10 

10 

10316.6 

20 

20 

1 0859. 5 


34 

47 

44 

24 

54 

41 

65 

58 

77 

72 

86 

83 

61 

70 

97 

92 

94 

no 

104 

115 

101 

119 

132 
102 

133 
150 
169 
146 
154 
136 
153 
125 
160 
165 
172 
176 
159 
183 
181 
192 
195 
200 
202 

206 

cr 


123456 

5 

1 


60 

0.950 

1.010 

8 

1 33 . 33 
175 
5 

-14.88 

138 

5 

- 12.1 3 


34 

31 

38 

48 

42 

43 

46 

56 

50 

53 

54 

33 

58 

64 

62 

50 

66 

75 

70 

68 

74 

87 

78 

82 

82 

60 

86 

88 

90 

96 

94 

81 

98 

93 

102 

109 

106 

114 

110 

118 

114 

1 1 1 

118 

1 16 

122 

130 

126 

112 

130 

129 

134 

157 

138 

113 

142 

140 

146 

167 

150 

147 

154 

155 

158 

137 

162 

161 

166 

164 

170 

173 

174 

177 

178 

171 

182 

184 

186 

187 

190 

193 

194 

197 

198 

201 

202 

203 

206 

138 

17 

27 

23 

54 

186 

192 

191 

193 

198 


199 

THRU 


0. 91639E+04 


177 163 

7.79 4.032 

1 40 1 29 

17.14 4.675 


35 

23 

39 

46 

43 

40 

47 

57 

51 

52 

55 

66 

59 

63 

63 

42 

67 

74 

71 

59 

75 

85 

79 

78 

83 

79 

87 

99 

91 

95 

95 

80 

99 

89 

103 

108 

107 

106 

111 

117 

115 

121 

119 

127 

123 

139 

127 

122 

131 

131 

135 

144 

139 

123 

143 

141 

147 

179 

151 

151 

155 

166 

159 

148 

163 

162 

167 

168 

171 

174 

175 

126 

179 

180 

183 

185 

187 

190 

191 

194 

195 

198 

199 

189 

203 

204 

36 

14 

37 

1 12 

187 

188 

194 

195 

205 



0.28149E+01 


166 156 

0.322 2.085 

131 121 

0.108 3.508 


36 

14 

40 

45 

44 

32 

48 

55 

52 

49 

56 

67 

60 

62 

64 

76 

68 

73 

72 

51 

76 

84 

80 

69 

84 

71 

88 

98 

92 

91 

96 

1 03 

100 

100 

104 

107 

108 

105 

1 12 

90 

1 16 

120 

120 

128 

124 

145 

128 

134 

132 

142 

136 

156 

140 

135 

144 

143 

148 

124 

152 

152 

156 

178 

160 

158 

164 

163 

168 

170 

172 

175 

176 

149 

180 

182 

184 

186 

188 

191 

192 

188 

196 

199 

200 

196 

204 

205 

10 

45 

146 

185 

189 

190 

196 

197 


O. 91790E-07 


O. 786.9179- 
0. 827. 9179- 


*AERO 


7+2 10 

7+2 20 


CRTOAL Page fs 

OF POOR QUALITY 
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STREAMLl 
STREAML2 
+2 30 
STREAMLl 
STREAML2 
+2 40 
STREAMLl 
STREAML2 
+2 50 
STREAMLl 
STREAML2 
+2 60 
STREAMLl 
STREAML2 
+2 70 


30 

30 

1 1513. 1 
40 
40 

10848. 5 
50 
50 

9745.6 

60 

60 

9163.9 

70 

70 

10952.6 


99 

5 

6.97 

64 

5 

26.04 

37 

5 

40.02 

18 

5 

47. 18 
1 
5 

40.07 


FLFACT 1 1.0 

FLFACT 2 -45.0 


FLFACT 

3 

. 1 

FLUTTER 

1 

KE 

MKAER02 

-45.0 

.05 

MKAER02 

-45. 0 

. 1 

MKAER02 

-45. 0 

. 3 

MKAER02 

-45.0 

.5 

MKAER02 

-45. 0 

.7 

MKAER02 

-45. 0 

.9 

MKAERG2 

-45. 0 

1.2 

MKAER02 

-45.0 

1.5 

MKAER02 

-45.0 

1.55 

PARAM 

MTYPE 

COSINE 

PARAM 

K INDEX 

0 

PARAM 

PRINT 

YESB 

PARAM 

CTYPE 

ROT 


ENDDATA 


101 
18. 27 

75 

18.50 

49 

21. 10 

21 

24.78 

13 

31.02 


103 

4.876 

1 05 
-0. 178 

1 1 1 
4.955 

0.877.9179- 

7+2 

77 

4.529 

92 

-0.312 

108 

6.339 

0.826.9179- 

7+2 

60 
3. 799 

70 

-0. 408 

82 
7. 703 

0.742.9179- 

7+2 

33 

2.815 

44 

-0. 570 

54 

8.894 

0. 698.9179- 

7+2 

14 

1.805 

1 5 

-0.570 

27 

9,716 

0.834.9179- 

7+2 


-■7 .9 1.2 1.5 

3 L 5 


30 


40 


50 


60 


70 


ORIGINAL PAGE IS 
OF POOR QUALITY 


100 



input data deck for 


APPLIED OSCILLATORY AIRLOADS GENERATION 


OSCILLATORY AIRLOADS FOR SR3 ADVANCED TURBOPROP 
CASE 6, TEST READING NO. 273, MACH .798 


NLINES 

7 

FLOTYP 

UNIFORM 10476 

NSEGS 

8 

RPS 

133. 33 

INCANG 

2. 0 

SSOUND 

13128. 

I REF 

60 

NASOUT 

YES 

MXMACH 

.95 

MNMACH 

1.01 


STREAML310 

STREAML320 

STREAML330 

STREAML340 

STREAML350 

STREAML3AO 

STREAML370 

STREAML410 

STREAML420 

STREAML430 

STREAML440 

STREAML450 

STREAML4A0 

STREAML470 

C0RD2R 77 

+C2R 10. 

GRDSET 

GRID 

GRID 

GRID 


175 
138 
99 
64 
37 
18 
1 

9, 179E-8 
9. 179E-8 
9. 179E-8 
9, 179E-8. 

9. 179E-8 
9, 179E-8 
9. 179E-8 
O .0 

-0.618 .0 
77 


177 

1 ^3 

166 

156 

140 

129 

131 

121 

101 

103 

105 

111 

75 

77 

92 

108 

49 

60 

70 

82 

21 

33 

44 

54 

13 

14 

15 

27 


.0 


1.808 

2.376 

2.625 


1 . 839 
2.347 
2.558 


12.250 

12.250 

12.250 


.0 


GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 


GRID DATA IDENTICAL TO THAT FOR FORCED RESPONSE ANALYSIS 


203 

204 

205 

206 
2 

11 
21 
30 
39 
48 
57 
66 
75 
88 
101 
1 14 
127 
140 
150 
160 


0.233 

0. 030 

0. 466 

0. 061 

0.699 

0.091 

-1. 496 

-0. 226 

2. 129 

2. 133 

1.541 

1 . 520 

1.034 

1 . 014 

0. 550 

O. 564 

0. 101 

0. 182 

“0.317 

-0. 138 

-0.677 

-0.383 

-0.998 

-0. 572 

-1.271 

-0. 702 

- 1 . 490 

-0. 779 

-1.621 

-O. 796 

-1 . 642 

-0 - 753 

-1.572 

-0.662 

-1 - 435 

-0.538 

-1. 208 

--0. 374 

-0. 917 

-0. 229 


1.920 

1-920 

1.920 

2. 650 
12.250 
1 1 . 600 
1 1 . 000 
10. 400 

9.800 
9. 187 
8.600 
8 . 000 
7.400 
6. 800 
6 . 200 
5 . 600 
5. 000 
4, 400 
3.715 

3. 180 


+C2R 
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INPUT DATA DECK FOR 
MODAL FORCED RESPONSE ANALYSIS 


NASTRAN BANDIT 


-1, SYSTEM (93) =1, FILES 


(INPT, PLT2) 


EXECUTIVE CONTROL DECK 


ID 

APR 

SOL 

DIAG 

TIME 

% 


NASA, SR3PR0P 
DISP 
8 

14,21,22 
10 % CRAY-1 S 


readfile mfvaaet 

* 

CEND 
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CASE CONTROL DECK 


% 

TITLE = SR3 RESPONSE TO 1 PER REV OSC. AIRLOADS 
SUBTITLE = NASA TEST READING NO. 273 
LABEL = RESPONSE ANALYSIS 
% 

SPC = 1 

MPC = I 

METHOD = 1 

$ 

D I SP ( SORT 1, PHASE) = ALL 
STRESS (SORTl, PHASE) = ALL 

$ 

$ 

% NOTE 

^ THE FOLLOWING DATA IS FOR 

$ HUBTYP <= O < 0 = RIGID, 1 = FLEXIBLE ) 

* FLOTYP =0 < 0 = UNIFORM, 1 = NON-UNIFORM ) 

% 

FREQUENCY = 1 * SID OF FREQ BULK DATA CARD 

% 

SUBCASE 1 

LABEL = K = O MODES, OSCILLATORY AIRLOADS PRESENT 
DLOAD = 1000 

$ 

BEGIN BULK 


OUTPUT FROM 

AIRLOADS 

PROGRAM 
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BULK DATA DECK 


C0RD2R 77 
+C2R 10. 
CTRIA2 1 
CTRIA2 2 
CTRIA2 3 
CTRIA2 4 
CTRIA2 5 
CTRIA2 6 
CTRIA2 7 
CTRIA2 8 
CTRIA2 V 
CTRIA2 10 
CTRIA2 11 
CTRIA2 12 
CTRIA2 13 
CTRIA2 14 
CTRIA2 15 
CTRIA2 16 
CTRIA2 17 
CTRIA2 18 
CTRIA2 19 
CTRIA2 20 
CTRIA2 21 
CTRIA2 22 
CTRIA2 23 
CTRIA2 24 
CTRIA2 25 
CTRIA2 26 
CTRIA2 27 
CTRIA2 28 
CTRIA2 29 
CTRIA2 30 
CTRIA2 31 
CTRIA2 32 
CTRIA2 33 
CTRIA2 34 
CTRIA2 35 
CTRIA2 36 
CTRIA2 37 
CTRIA2 38 
CTRIA2 39 
CTRIA2 40 
CTRIA2 41 
CTRIA2 42 
CTRIA2 43 
CTRIA2 44 
CTRIA2 45 
CTRIA2 46 
CTRIA2 47 
CTRIA2 48 
CTRIA2 49 
CTRIA2 50 
CTRIA2 51 
CTRIA2 52 
CTRIA2 53 
CTRIA2 54 
CTRIA2 55 
CTRIA2 56 
CTRIA2 57 
CTRIA2 58 
CTRIA2 59 
CTRIA2 60 
CTRIA2 61 
CTRIA2 62 


0 

. 0 

-0. 

618 .0 

1 

10 

2 

11 

3 

8 

4 

12 

5 

7 

6 

13 

7 

1 

8 

14 

9 

2 

10 

15 

11 

3 

12 

16 

13 

4 

14 

17 

15 

5 

16 

20 

17 

21 

18 

18 

19 

22 

20 

9 

21 

10 

22 

23 

23 

11 

24 

24 

25 

12 

26 

25 

27 

13 

28 

14 

29 

26 

30 

15 

31 

27 

32 

16 

33 

29 

34 

30 

35 

19 

36 

31 

37 

20 

38 

32 

39 

21 

40 

33 

41 

22 

42 

23 

43 

34 

44 

24 

45 

35 

46 

25 

47 

36 

48 

26 

49 

38 

50 

28 

51 

39 

52 

40 

53 

29 

54 

30 

55 

41 

56 

31 

57 

42 

58 

32 

59 

43 

60 

•r-?- 

61 

44 

62 

34 


. O 


9 

8 

10 

8 

7 

1 1 

11 

7 

1 

12 

12 

1 

2 

13 

13 

2 

3 

14 

14 

3 

4 

15 

15 

4 

5 

16 

16 

5 

6 

17 

19 

18 

20 

18 

9 

21 

21 

9 

10 

22 

11 

22 

22 

11 

12 

23 

23 

12 

13 

24 

24 

13 

14 

25 

15 

25 

25 

15 

16 

26 

26 

16 

17 

27 

28 

19 

29 

19 

20 

30 

30 

20 

21 

31 

31 

21 

22 

32 

32 

22 

23 

33 

24 

33 

33 

24 

25 

34 

34 

25 

26 

35 

35 

26 

27 

36 

37 

28 

29 

39 

38 

28 

39 

29 

30 

40 

31 

40 

40 

31 

32 

41 

41 

32 

3^ 

42 

42 

33 

34 

43 

43 

34 



CTRIA2 63 

63 

CTRIA2 64 

64 

CTRIA2 65 

65 

CTRIA2 66 

66 

CTRIA2 67 

67 

CTRIA2 68 

68 

CTRIA2 69 

69 

CTRIA2 70 

70 

CTRIA2 71 

71 

CTRIA2 72 

72 

CTRIA2 73 

73 

CTRIA2 74 

74 

CTRIA2 75 

75 

CTRIA2 76 

76 

CTRIA2 77 

77 

CTRIA2 78 

78 

CTRIA2 79 

79 

CTRIA2 80 

80 

CTRIA2 81 

81 

CTRIA2 82 

82 

CTRIA2 83 

83 

CTRIA2 84 

84 

CTRIA2 85 

85 

CTRIA2 86 

86 

CTRIA2 87 

87 

CTRIA2 88 

88 

CTRIA2 89 

89 

CTRIA2 90 

90 

CTRIA2 91 

91 

CTRIA2' 92 

92 

CTRIA2 93 

93 

CTRIA2 94 

94 

CTRIA2 95 

95 

CTRIA2 96 

96 

CTRIA2 97 

97 

CTRIA2 98 

98 

CTRIA2 99 

99 

CTRIA2 100 

100 

CTRIA2 101 

101 

CTRIA2 102 

102 

CTRIA2 103 

103 

CTRIA2 104 

104 

CTRIA2 105 

105 

CTRIA2 106 

106 

CTRIA2 107 

107 

CTRIA2 108 

108 

CTRIA2 109 

109 

CTRIA2 110 

110 

CTRIA2 111 

111 

CTRIA2 112 

112 

CTRIA2 113 

113 

CTRIA2 114 

114 

CTRIA2 US 

115 

CTRIA2 116 

116 

CTRIA2 117 

117 

CTRIA2 118 

118 

CTRIA2 119 

119 

CTRIA2 120 

120 

CTRIA2 121 

121 

CTRIA2 122 

122 

CTRIA2 123 

123 

CTRIA2 124 

124 

CTRIA2 125 

125 

CTRIA2 126 

126 

CTRIA2 127 

127 

nTRTA2 128 

1 28 


45 

44 

3S 

35 

36 

45 

47 

46 

37 

37 

38 

47 

48 

47 

38 

38 

39 

48 

49 

48 

39 

39 

40 

49 

50 

49 

40 

40 

41 

50 

51 

50 

41 

41 

42 

51 

52 

51 

42 

42 

43 

52 

53 

52 

43 

43 

44 

53 

54 

53 

44 

44 

45 

54 

56 

55 

46 

46 

47 

56 

57 

56 

47 

47 

48 

57 

58 

57 

48 

48 

49 

58 

59 

58 

49 

49 

50 

59 

60 

59 

50 

50 

51 

60 

61 

60 

51 

51 

52 

61 

62 

61 

52 

52 

53 

62 

63 

62 

53 

53 

54 

63 

65 

64 

55 

55 

56 

65 

66 

65 

56 

56 

57 

66 

67 

66 

57 

57 

58 

67 

68 

67 

58 

58 

59 

68 

69 

68 

59 

59 

60 

69 

70 

69 

60 

60 

61 

70 

71 

70 

61 

61 

62 

71 

72 

71 

62 

62 

63 

72 

74 

73 

64 

64 

65 

74 

75 

74 

65 

65 

66 

75 

76 

75 

66 

66 

67 

76 

77 

76 

67 

67 

68 

77 

78 

77 

68 

68 

69 

78 

79 

78 

69 

69 

70 

79 

80 

79 

70 

70 

71 

80 

81 

80 

71 

71 

72 

81 
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CTRIA2 

129 

129 

CTRIA2 

130 

130 

CTRIA2 

131 

131 

CTRIA2 

132 

132 

CTRIA2 

133 

133 

CTRIA2 

134 

134 

CTRIA2 

135 

135 

CTRIA2 

136 

136 

CTRIA2 

137 

137 

CTRIA2 

138 

138 

CTRIA2 

139 

139 

CTRIA2 

140 

140 

CTRIA2 

141 

141 

CTRIA2 

142 

142 

CTRIA2 

143 

143 

CTRIA2 

144 

144 

CTRIA2 

145 

145 

CTRIA2 

146 

146 

CTRIA2 

147 

147 

CTRIA2 

148 

148 

CTRIA2 

149 

149 

CTRIA2 

150 

150 

CTRIA2 

151 

151 

CTRIA2 

152 

152 

CTRIA2 

153 

153 

CTRIA2 

154 

154 

CTRIA2 

155 

155 

CTRIA2 

156 

156 

CTRIA2 

157 

157 

CTRIA2 

158 

158 

CTRIA2 

159 

159 

CTRIA2 

160 

160 

CTRIA2 

161 

161 

CTRIA2 

162 

162 

CTRIA2 

163 

163 

CTRIA2 

164 

164 

CTRIA2 

165 

165 

CTRIA2 

166 

166 

CTRIA2 

167 

167 

CTRIA2 

168 

168 

CTRIA2 

169 

169 

CTRIA2 

170 

170 

CTRIA2 

171 

171 

CTRIA2 

172 

172 

CTRIA2 

173 

173 

CTRIA2 

174 

174 

CTRIA2 

175 

175 

CTRIA2 

176 

176 

CTRIA2 

177 

177 

CTRIA2 

178 

178 

CTRIA2 

179 

179 

CTRIA2 

180 

180 

CTRIA2 

181 

181 

CTRIA2 

182 

182 

CTRIA2 

183 

183 

CTRIA2 

184 

184 

CTRIA2 

185 

185 

CTRIA2 

186 

186 

CTRIA2 

187 

187 

CTRIA2 

188 

188 

CTRIA2 

189 

189 

CTRIA2 

190 

190 

CTRIA2 

191 

191 

CTRIA2 

192 

192 

nTRT AC? 

1 97; 

1 R3 


82 

81 

72 

87 

86 

73 

73 

74 

87 

88 

87 

74 

74 

75 

88 

89 

88 

75 

75 

76 

89 

90 

89 

76 

76 

77 

90 

91 

90 

77 

77 

78 

91 

92 

91 

78 

78 

79 

92 

79 

80 

83 

79 

83 


93 

92 

83 

80 

81 

84 

84 

83 

80 

83 

84 

94 

94 

93 

83 

81 

82 

85 

85 

84 

81 

84 

85 

95 

95 

94 

84 

100 

99 

86 

86 

87 

100 

101 

100 

87 

87 

88 

101 

102 

101 

88 

88 

89 

102 

103 

102 

89 

89 

70 

103 

104 

103 

90 

90 

91 

104 

105 

104 

91 

91 

92 

105 

92 

93 

96 

92 

96 

105 

106 

105 

96 

93 

94 

97 

97 

96 

93 

96 

97 

107 

107 

106 

96 

94 

95 

98 

98 

97 

94 

97 

98 

108 

108 

107 

97 

113 

112 

99 

99 

100 

113 

114 

113 

100 

100 

101 

114 

115 

114 

101 

101 

102 

115 

116 

115 

102 

102 

103 

116 

1 17 

116 

103 

103 

104 

1 17 

1 18 

117 

104 

104 

105 

1 18 

105 

106 

109 

105 

109 

118 

119 

1 18 

109 

106 

107 

1 10 

1 10 

109 

106 

1 (■)<? 

1 1 O 

1 20 
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UIH1A2 

194 

CTRIA2 

195 

CTRIA2 

196 

CTRIA2 

197 

CTRIA2 

198 

CTRIA2 

199 

CTRIA2 

200 

CTRIA2 

201 

CTRIA2 

202 

CTRIA2 

203 

CTRIA2 

204 

CTRIA2 

205 

CTRIA2 

206 

CTRIA2 

207 

CTRIA2 

208 

CTRIA2 

209 

CTRIA2 

210 

CTRIA2 

211 

CTRIA2 

212 

CTRIA2 

213 

CTRIA2 

214 

CTRIA2 

215 

CTRIA2 

216 

CTRIA2 

217 

CTRIA2 

218 

CTRIA2 

219 

CTRIA2 

220 

CTRIA2 

221 

CTRIA2 

222 

CTRIA2 

223 

CTRIA2 

224 

CTRIA2 

225 

CTRIA2 

226 

CTRIA2 

227 

CTRIA2 

228 

CTRIA2 

229 

CTRIA2 

230 


CTRIA2 231 
CTRIA2 232 
CTRIA2 233 
CTRIA2 234 
CTRIA2 235 


CTRIA2 

236 

CTRIA2 

237 

CTRIA2 

238 

CTRIA2 

239 

CTRIA2 

240 

CTRIA2 

241 

CTRIA2 

242 

CTRIA2 

243 

CTRIA2 

244 

CTRIA2 

245 

CTRIA2 

246 

CTRIA2 

247 

CTRIA2 

248 

CTRIA2 

249 

CTRIA2 

250 

CTRIA2 

251 

CTRIA2 

252 

CTRIA2 

253 

CTRIA2 

254 

CTRIA2 

255 

CTRIA2 

256 

CTRIA2 

257 

CTRIA2 

258 

CTRIA2 

259 


194 

120 

195 

107 

196 

1 1 1 

197 

1 10 

198 

121 

199 

112 

200 

126 

201 

113 

202 

127 

203 

114 

204 

128 

205 

1 15 

206 

129 

207 

116 

208 

130 

209 

117 

210 

131 

■ 211 

118 

212 

118 

213 

132 

214 

119 

215 

123 

216 

122 

217 

133 

218 

120 

219 

124 

220 

123 

221 

134 

222 

125 

223 

139 

224 

126 

225 

140 

226 

127 

227 

141 

228 

128 

229 

142 

230 

129 

231 

143 

232 

130 

233 

144 

234 

131 

235 

131 

236 

145 

237 

132 

238 

136 

239 

135 

240 

146 

241 

133 

242 

137 

243 

136 

244 

147 

245 

138 

246 

149 

247 

139 

248 

150 

249 

140 

250 

151 

251 

141 

252 

152 

253 

142 

254 

153 


143 

256 

154 

257 

144 

258 

155 

259 

145 


1 19 
108 
1 10 
111 
120 

113 

125 

114 

126 

1 15 

127 

1 16 

128 

117 

129 

118 

130 

119 

131 

120 
122 
123 

132 
121 

123 

124 

133 
126 

138 

127 

139 

128 

140 

129 

141 

130 

142 

131 

143 


135 

144 

133 

135 

136 

145 

134 

136 

137 

146 

139 

148 

140 

149 

141 

150 

142 

151 

143 

152 

144 

153 

145 

154 

146 


109 
111 
107 
121 

1 10 

125 

113 

126 

114 

127 
1 15 

128 
116 

129 

117 

130 

1 18 
122 

131 
122 
122 
120 

132 
123 

123 
121 

133 

124 

138 
126 

139 

127 

140 

128 

141 

129 

142 

130 

143 

131 
135 

144 
135 

135 

133 

145 

136 

136 

134 

146 

137 

148 

139 

149 

140 

150 

141 

151 

142 

152 

143 

153 

144 

154 

145 
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CTRIA2 260 

260 

CTRIA2 261 

261 

CTRIA2 262 

262 

CTRIA2 263 

263 

CTRIA2 264 

264 

CTRIA2 265 

265 

CTRIA2 266 

266 

CTRIA2 267 

267 

CTRIA2 268 

268 

CTRIA2 269 

269 

CTRIA2 270 

270 

CTRIA2 271 

271 

CTRIA2 272 

272 

CTRIA2 273 

273 

CTRIA2 274 

274 

CTRIA2 275 

275 

CTRIA2 276 

276 

CTRIA2 277 

277 

CTRIA2 278 

278 

CTRIA2 279 

279 

GTRIA2 280 

280 

CTRIA2 281 

281 

CTRIA2 282 

282 

CTRIA2 283 

283 

CTRIA2 284 

284 

CTRIA2 285 

285 

CTRIA2 286 

286 

CTRIA2 287 

287 

CTRIA2 288 

288 

CTRIA2 289 

289 

CTRIA2 290 

290 

CTRIA2 291 

291 

CTRIA2 292 

292 

CTRIA2 293 

293 

CTRIA2 294 

294 

CTRIA2 295 

295 

CTRIA2 296 

296 

CTRIA2 297 

297 

CTRIA2 298 

298 

CTRIA2 299 

299 

CTRIA2 300 

300 

CTRIA2 301 

301 

CTRIA2 302 

302 

CTRIA2 303 

303 

CTRIA2 304 

304 

CTRIA2 305 

305 

CTRIA2 306 

306 

CTRIA2 307 

307 

CTRIA2 308 

308 

CTRIA2 309 

309 

CTRIA2 310 

310 

CTRIA2 311 

311 

CTRIA2 312 

312 

CTRIA2 313 

313 

CTRIA2 314 

314 

CTRIA2 315 

315 

CTRIA2 316 

316 

CTRIA2 317 

317 

CTRIA2 318 

318 

CTRIA2 319 

319 

CTRIA2 320 

320 

CTRIA2 321 

321 

CTRIA2 322 

322 

CTRIA2 323 

323 

CTRTA2 324 

32^ 


156 

155 

146 

146 

147 

1 56 

148 

149 

157 

158 

157 

149 

149 

150 

158 

159 

158 

150 

160 

159 

150 

150 

151 

160 

161 

160 

151 

151 

152 

162 

162 

161 

151 

163 

162 

152 

164 

163 

152 

152 

153 

164 

165 

164 

153 

166 

165 

153 

153 

154 

166 

167 

166 

154 

154 

155 

167 

155 

156 

167 

157 

158 

175 

206 

175 

158 

158 

159 

206 

176 

206 

159 

159 

160 

176 

176 

177 

160 

160 

161 

168 

168 

177 

160 

168 

178 

177 

169 

168 

161 

161 

162 

169 

170 

169 

162 

162 

163 

170 

171 

170 

163 

172 

171 

163 

163 

164 

172 

173 

172 

164 

164 

165 

173 

174 

173 

165 

168 

169 

179 

179 

178 

168 

169 

170 

180 

180 

179 

169 

170 

171 

181 

181 

180 

170 

171 

172 

182 

182 

181 

171 

172 

173 

183 

183 

182 

172 

173 

174 

184 

184 

183 

173 

178 

179 

185 

186 

185 

179 

179 

180 

186 

187 

186 

180 

180 

181 

187 

188 

187 

181 

181 

182 

188 

189 

188 

182 

182 

183 

189 

190 

189 

183 

183 

184 

190 

191 

190 

184 

185 

186 

192 

193 

1 92 

1 RA 
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CTRIA2 

325 


325 

CTRIA2 

326 


326 

CTRIA2 

327 


327 

CTRIA2 

328 


328 

CTRIA2 

329 


329 

CTRIA2 

330 


330 

CTRIA2 

331 


331 

CTRIA2 

332 


332 

CTRIA2 

333 


333 

CTRIA2 

334 


334 

CTRIA2 

335 


335 

CTRIA2 

336 


336 

CTRIA2 

337 


337 

CTRIA2 

338 


338 

CTRIA2 

339 


339 

CTRIA2 

340 


340 

CTRIA2 

341 


341 

CTRIA2 

342 


342 

CTRIA2 

343 


343 

CTRIA2 

344 


344 

CTRIA2 

345 


345 

CTRIA2 

346 


346 

CYJOIN 

1 



CYJOIN 

2 



EIGR 

1 


FEER 

+EIG1 

MAX 



GRDSET 



77 

GRID 


1 


GRID 


3 


GRID 


4 


GRID 


5 


GRID, 


6 


GRID 


7 


GRID 


8 


GRID 


9 


GRID 


10 


GRID 


12 


GRID 


13 


GRID 


14 


GRID 


15 


GRID 


16 


GRID 


17 


GRID 


18 


GRID 


19 


GRID 


20 


GRID 


22 


GRID 


23 


GRID 


24 


GRID 


25 


GRID 


26 


GRID 


27 


GRID 


28 


GRID 


29 


GRID 


31 


GRID 


32 


GRID 


33 


GRID 


34 


GRID 


35 


GRID 


36 


GRID 


37 


GRID 


38 


GRID 


40 


GRID 


41 


GRID 


42 


GRID 


43 


GRID 


44 



186 

187 

193 

194 

193 

187 

187 

188 

194 

195 

194 

188 

188 

189 

195 

196 

195 

189 

189 

190 

196 

197 

196 

190 

190 

191 

197 

198 

197 

191 

200 

199 

192 

192 

193 

200 

201 

200 

193 

193 

194 

201 

202 

201 

194 

194 

195 

202 

195 

196 

202 

203 

202 

196 

196 

197 

203 

204 

203 

197 

197 

198 

204 

205 

204 

198 

199 



205 




1.808 

1.839 

12.250 

2.376 

2.347 

12.250 

2.625 

2.558 

12.250 

2.877 

2.765 

12.250 

3. 134 

2.966 

12.250 

1.556 

1.589 

12.033 

1.304 

1.339 

11.817 

1.052 

1,088 

11. 600 

1.293 

1.308 

1 1 . 600 

1.791 

1.730 

1 1 . 600 

2.044 

1.937 

11.600 

2.298 

2. 141 

1 1 . 600 

2.555 

2.343 

1 1 . 600 

2.814 

2.542 

1 1 . 600 

3.078 

2.733 

1 1 . 600 

0.737 

0.793 

1 1 . 300 

0.423 

0. 498 

1 1 . 000 

0.725 

0.761 

1 1 . 000 

1.347 

1.264 

1 1 . 000 

1.663 

1.510 

1 1 . 000 

1.981 

1.753 

1 1 . 000 

2.302 

1.993 

11,000 

2.626 

2.228 

1 1 . 000 

2.956 

2.454 

1 1 . 000 

-0. 168 

-0.013 

10. 400 

0. 186 

0.281 

10.400 

0.917 

0.842 

10. 400 

1.288 

1.115 

10.400 

1.661 

1 . 385 

10. 400 

2.038 

1.650 

10.400 

2. 418 

1.910 

10. 400 

2. 806 

2. 158 

10. 400 

-0.702 

-0. 430 

9 . 800 

-0.306 

-0, 118 

9. 800 

0.512 

0. 476 

9 . 800 

0,926 

0.764 

9.800 

1.344 

1 , 048 

9. 800 

1.766 

1 . 326 

9 . 8(’)0 

2.191 

1 . 598 

9. 800 


+EIG1 
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GRID 

45 

2.625 

1.856 

9 . 800 

GRID 

46 

-1. 193 

-0. 766 

9. 187 

GRID 

47 

-0.760 

-0.445 

9, 187 

GRID 

49 

0. 130 

0. 162 

9. 187 

GRID 

50 

0.582 

0.456 

9. 187 

GRID 

51 

1.037 

0.745 

9. 187 

GRID 

52 

1.495 

1.027 

9. 187 

GRID 

53 

1.959 

1 . 303 

9. 187 

GRID 

54 

2.430 

1.563 

9.187 

GRID 

o5 

-1.612 

-1.013 

8 . 600 

GRID 

56 

-1. 149 

-0.691 

8.600 

GRID 

58 

-0. 199 

-0.083 

8.600 

GRID 

59 

0.282 

0.210 

8.600 

GRID 

60 

0.768 

0.498 

8.600 

GRID 

61 

1.257 

0.778 

8.600 

GRID 

62 

1.750 

1.051 

8.600 

GRID 

63 

2.252 

1.308 

8.600 

GRID 

64 

-1.985 

-1. 192 

8.000 

GRID 

65 

-1.496 

-0.874 

8.000 

GRID 

67 

-0.494 

-0.277 

8.000 

GRID 

68 

0.013 

0.010 

8.000 

GRID 

69 

0.525 

0.291 

8. 000 

GRID 

70 

1.040 

0.564 

8.000 

GRID 

71 

1.560 

0.829 

8.000 

GRID 

72 

2.086 

1.078 

8.000 

GRID 

73 

-2.301 

-1.303 

7.400 

GRID 

74 

-1.790 

-0.995 

. 7.400 

GRID 

76 

-0.747 

-0.418 

7. 400 

GRID 

77 

-0.218 

-0. 141 

7.400 

GRID 

78 

0.314 

0. 128 

7. 400 

GRID 

79 

0.850 

0.390 

7.400 

GRID 

80 

1.390 

0.642 

7.400 

GRID 

81 

1.664 

0.760 

7.400 

GRID 

82 

1.937 

0.878 

7. 400 

GRID 

83 

1.324 

0.567 

7. lOO 

GRID 

84 

1.602 

0.681 

7. 100 

GRID 

85 

1.880 

0.795 

7. 100 

GRID 

86 

-2.556 

-1.354 

6.800 

GRID 

87 

-2. 028 

-1.058 

6.800 

GRID 

89 

-0.948 

-0.508 

6.800 

GRID 

90 

-0.402 

-0.245 

6.800 

GRID 

91 

0. 147 

0.010 

6.800 

GRID 

92 

0.701 

0.256 

6 . 800 

GRID 

93 

1.258 

0.493 

6.800 

GRID 

94 

1.541 

0.603 

6.800 

GRID 

95 

1.823 

0.712 

6. 800 

GRID 

96 

1.227 

0.439 

6.500 

GRID 

97 

1.512 

0.544 

6.500 

GRID 

98 

1.798 

0.649 

6.500 

GRID 

99 

-2.716 

-1.339 

6.200 

GRID 

100 

-2. 173 

-1.060 

6.200 

GRID 

102 

-1.066 

-0.542 

6.200 

GRID 

103 

-0.507 

-0.296 

6.200 

GRID 

104 

0.057 

-0. 059 

6.200 

GRID 

105 

0.624 

0. 169 

6.200 

GRID 

106 

1. 195 

0.386 

6.200 

GR I D 

107 

1.484 

0.486 

6. 200 

GRID 

108 

1.772 

0.586 

6, 200 

GRID 

109 

1.195 

0.352 

5.900 

GR I D 

1 10 

' 1.484 

0. 447 

5. 900 

GRID 

111 

1.773 

0.542 

5 . 900 

GRID 

1 12 

-2.748 

-1.254 

5. 600 

GRID 

113 

-2, 198 

-0. 996 

5 . 600 

GRID 

1 15 

-1.083 

-0.519 

5 . 600 

RR T n 

1 1 A 

-O. 519 

-O. ‘794 

5 . 600 
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GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 


1 17 

0. 049 

-0.079 

5. 600 

118 

0. 620 

0. 125 

5. 600 

1 19 

1 . 195 

0.318 

5 . 6o0 

120 

1.485 

0. 408 

5 . 600 

121 

1.775 

0. 498 

5 . 6C)0 

122 

1.216 

0.301 

5.300 

123 

1.504 

0.386 

5.300 

124 

1 . 792 

0.471 

5. 300 

125 

-2.670 

-1.111 

5 . 000 

126 

-2. 124 

-0. 880 

5 . 0C>0 

128 

-1.017 

-0. 453 

5 . 000 

129 

-0.459 

-0.254 

5. 000 

130 

0. 104 

-0.065 

5.000 

131 

0.669 

0.114 

5. 000 

132 

1.238 

0.283 

5.000 

133 

1.523 

0.363 

5.000 

134 

1.809 

0. 444 

5.000 

135 

1.274 

0.278 

4.700 

136 

1.556 

0.355 

4. 700 

137 

1.837 

0.432 

4 . 700 

138 

“2.513 

-0.927 

4. 400 

139 

-1.976 

-0.727 

4 . 400 

141 

-0.891 

-0.359 

4. 400 

142 

-0.344 

-0. 189 

4.400 

143 

0.206 

-0.028 

4.400 

144 

0.758 

0. 125 

4.400 

145 

1.311 

0.274 

4.400 

146 

1.588 

0.347 

4.400 

147 

1.865 . 

0. 420 

4 . 400 

148 

-2.273 

-0.686 

3.715 

149 

-1.742 

-0.524 

3.715 

151 

-0.672 

-0.233 

3.715 

152 

-0. 133 

-0.099 

3.715 

153 

0.406 

0,031 

3.715 

154 

0.946 

0. 159 

3.715 

155 

1.485 

0.291 

3.715 

156 

2.021 

0.433 

3.715 

157 

-2.051 

-0.483 

3. 180 

158 

-1.675 

-0.394 

3. 180 

159 

-1.296 

-0,309 

3. 180 

161 

-0.548 

-0.089 

3. 180 

162 

-0.274 

-0.044 

3. 180 

163 

0.000 

0.000 

3. 180 

164 

0.274 

0.044 

3. 180 

165 

0.548 

0.089 

3. 180 

166 

1,037 

0.214 

3.358 

167 

1.527 

0.328 

3.537 

168 

-0.548 

-0.089 

2.930 

169 

-0.365 

-0.059 

2.930 

170 

-O. 183 

-0.030 

2.930 

171 

0.000 

0.000 

2.930 

172 

0. 183 

0.030 

2.930 

173 

0.365 

0.059 

2.930 

174 

0.548 

0. 089 

2.930 

175 

-1.804 

-0.270 

2.650 

176 

-1. 188 

-0. 182 

2.650 

177 

-0. 750 

-0. 123 

2.740 

178 

-0.550 

-0.072 

2.600 

179 

-0.367 

-0.048 

2 . 600 

180 

-0. 184 

-0 . 024 

2 . 600 

181 

0.000 

0 . 000 

2. 600 

182 

0. 184 

0. 024 

2. 600 

183 

0.367 

0.048 

2 . 600 

184 

0.550 

0.072 

2. 600 

185 

-0.550 

-0. 072 

2. 350 

186 

-0. 367 

-0. 048 

2 . 350 
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GRID 

grid 

GRID 

GRID 

GRID 

GRID 

GRID 

grid 

GRID 

GRID 

GRID 

grid 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

grid 

GRID 

GRID 

GRID 

GRID 

GRID 

grid 

GRID 

GRID 

GRID 

GRID 

MATl 

MPC 

MPC 

MPC 

PTRIA2 

PTRIA2 

PTRIA2 

PTRIA2 

PTRIA2 

PTRIA2 

PTRIA2 

PTRIA2 

PTRIA2 

PTRIA2 

PTRIA2 

PTRIA2 

PTRIA2 

PTRIA2 

PTRIA2 

PTRIA2 

PTRIA2 

PTRIA2 

PTRIA2 

PTRIA2 

PTRIA2 

PTRIA2 

PTRIA2 

PTRIA2 

PTRIA2 


1 

1 

1 

1 

1 

2 

3 

4 

vJ 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 


24 


187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 
2 

11 

21 

30 

39 

48 

57 

66 

75 

88 

101 

114 

127 

140 

150 

160 

1.6 

5 

6 
7 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

t 


E7 


-0. 184 

-0. 024 

2.350 

0. 000 

0. 000 

2.350 

0. 184 

0.024 

2.350 

0.367 

0. 048 

2.350 

0. 550 

0.072 

2.350 

-0. 550 

-0.072 

2.070 

-0.367 

-0.048 

2.070 

-0. 184 

-0.024 

2.070 

0. 000 

0. 000 

2.070 

0. 184 

0.024 

2.070 

0.367 

0.048 

2.070 

0.550 

0.072 

2.070 

-0.699 

-0.091 

1.920 

-0.466 

-0.061 

1.920 

-0.233 

-0.030 

1.920 

0.000 

0.000 

1.920 

0.233 

0.030 

1.920 

0.466 

0. 061 

1.920 

0.699 

0.091 

1.920 

-1.496 

-0.226 

2.650 

2. 129 

2. 133 

12.250 

1.541 

1.520 

1 1 . 600 

1.034 

1.014 

1 1 . 000 

0.550 

0.564 

10.400 

0. 101 

0. 182 

9.800 

-0.317 

-0. 138 

9. 187 

-0.677 

-0. 383 

8.600 

-0.998 

-0.572 

8.000 

-1.271 

-0.702 

7.400 

-1.490 

-0.779 

6.800 

-1.621 

-0. 796 

6.200 

-1.642 

-0.753 

5. 600 

-1.572 

-0.662 

5.000 

-1. 435 

-0.538 

4.400 

-1.208 

-0.374 

3.715 

-0.917 

-0.229 

3. 180 

.35 .0004141 

1. 

0 4 

4 

1. 

0 4 

4 

1. 

0 1 

4 


.01570 
. 02827 
.01897 
. 03380 
. 02043 
- 03623 


.02917 
. 04440 
. 03830 
. 04253 
. 03677 
. 03397 
. 02740 
.01673 
. 00823 
.01970 
. 03550 
. 02390 


• 04250 
. 03487 
. 04743 
.05847 
. 05413 
. 06033 
. 05500 


OF 


IS 


quality 


- 1.0 

- 1.0 

- 1.0 
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PTRIA2 26 

1 

• 05663 

PTRIA2 27 

1 

. 05230 

PTRIA2 28 

1 

. 04877 

PTRIA2 29 

1 

. 04293 

PTRIA2 30 

1 

. 03390 

PTRIA2 31 

1 

, 02023 

PTRIA2 32 

1 

. 00963 

PTRIA2 33 

1 

. 02363 

PTRIA2 34 

1 

. 04273 

PTRIA2 35 

1 

. 03977 

PTRIA2 36 

1 

. 06260 

PTRIA2 37 

1 

. 05863 

PTRIA2 38 

1 

.07160 

PTRIA2 39 

1 

. 06707 

PTRIA2 40 

1 

.07400 

PTRIA2 41 

1 

. 06923 

PTRIA2 42 

1 

.06910 

PTRIA2 43 

1 

. 06850 

PTRIA2 44 

1 

. 06057 

PTRIA2 45 

1 

. 05280 

PTRIA2 46 

1 

. 04220 

PTRIA2 47 

1 

. 02500 

PTRIA2 48 

1 

.01207 

PTRIA2 49 

1 

. 02780 

PTRIA2 50 

1 

. 04730 

PTRIA2 51 

1 

. 05020 

PTRIA2 52 

1 

. 07393 

PTRIA2 53 

1- 

. 07003 

PTRIA2 54 

1 . 

. 07890 

PTRIA2 55 

1 

. 08777 

PTRIA2 56 

1 

. 08443 

PTRIA2 57 

1 

. 08887 

PTRIA2 58 

1 

. 08297 

PTRIA2 59 

1 

.08140 

PTRIA2 60 

1 

. 07273 

PTRIA2 61 

1 

.06277 

PTRIA2 62 

1 

. 05077 

PTRIA2 63 

1 

. 02983 

PTRIA2 64 

1 

.01467 

PTRIA2 65 

1 

. 03187 

PTRIA2 66 

1 

. 04680 

PTRIA2 67 

1 

. 07307 

PTRIA2 68 

1 

. 07707 

PTRIA2 69 

1 

.09197 

PTRIA2 70 

1 

.09173 

PTRIA2 71 

1 

. 10140 

PTRIA2 72 

1 

. 09827 

PTRIA2 73 

1 

. 10277 

PTRIA2 74 

1 

. 09670 

PTRIA2 75 

1 

. 09440 

PTRIA2 76 

1 

. 08500 

PTRIA2 77 

1 

. 07300 

PTRIA2 78 

1 

. 05943 

PTRIA2 79 

1 

. 03483 

PTRIA2 80 

1 

.01723 

PTRIA2 81 

1 

. 03637 

PTRIA2 82 

1 

. 05337 

PTRIA2 83 

1 

.08310 

PTRIA2 84 

1 

. 08787 

PTRIA2 85 

1 

. 10477 

PTRIA2 86 

1 

. 10487 

PTRIA2 87 

1 

. 1 1580 

PTRIA2 88 

1 

. 11257 

PTRIA2 89 

1 

. 11757 

F'TRIA2 90 

1 

. 1 1090 

PTRTA2 91 

1 

. 1 OR 1 7 


PACIE 3 
'OOR QUALITV' 
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PTRIA2 

92 

1 

.09773 

A 

PTRIA2 

93 

1 

. 08390 

PTRIA2 

94 

1 

.06857 


PTRIA2 

95 

1 

.04017 


PTRIA2 

96 

1 

. 02000 


PTRIA2 

97 

1 

.04190 


PTRIA2 

PTRIA2 

98 

99 

1 

1 

. 06107 
. 09537 

OR/GJAJAt Pantr . 
OF Pnr,n „ S 

quality 

PTRIA2 

PTRIA2 

100 

101 

1 

1 

. 10040 
. 12033 

PTRIA2 

102 

1 

.11993 


PTRIA2 

103 

1 

. 13323 


PTRIA2 

104 

1 

. 12900 


PTRIA2 

105 

1 

. 13550 


PTRIA2 

106 

1 

. 12730 


PTRIA2 

107 

1 

. 12473 


PTRIA2 

108 

1 

. 11227 


PTRIA2 

109 

1 

. 09677 


PTRIA2 

110 

1 

. 07883 


PTRIA2 

111 

1 

. 04630 


PTRIA2 

112 

1 

. 02297 


PTRIA2 

113 

1 

. 04790 


PTRIA2 

114 

1 

. 07000 


PTRIA2 

115 

1 

. 10907 


PTRIA2 

116 

1 

. 11527 


PTRIA2 

117 

1 

. 13790 


PTRIA2 

118 

1 

. 13793 


PTRIA2 

119 

1 

. 15300 


PTRIA2 

120 

1 

. 14870 


PTRIA2 

121 

1 

. 15617 


PTRIA2 

122 

1 

. 14720 


PTRIA2 

123 

1 

. 14463 


PTRIA2 

124 

1 

. 13037 


PTRIA2 

125 

1 

.11300 


PTRIA2 

126 

1 

.09187 


PTRIA2 

127 

1 

. 06993 


PTRIA2 

128 

1 

. 04230 


PTRIA2 

129 

1 

.01937 


PTRIA2 

130 

1 

.05410 


PTRIA2 

131 

1 

. 07907 


PTRIA2 

132 

1 

. 12273 


PTRIA2 

133 

1 

. 13063 


PTR1A2 

134 

1 

. 15583 


PTRIA2 

135 

1 

. 15690 


PTRIA2 

136 

1 

. 17377 


PTRIA2 

137 

1 

. 16983 


PTRIA2 

138 

1 

. 17853 


PTRIA2 

139 

1 

. 16917 


PTRIA2 

140 

1 

. 16647 


PTRIA2 

141 

1 

. 15090 


PTRIA2 

142 

1 

. 10750 


PTRIA2 

143 

1 

. 13113 


PTRIA2 

144 

1 

. 12383 


PTRIA2 

145 

1 

. 06393 


PTRIA2 

146 

1 

. 08037 


PTRIA2 

147 

1 

.07120 


PTRIA2 

148 

1 

.09103 


PTRIA2 

149 

1 

.01970 


PTRIA2 

150 

1 

. 03663 


PTRIA2 

151 

1 

.02113 


PTRIA2 

152 

1 

. 03927 


PTRIA2 

153 

1 

.06227 


PTRIA2 

154 

1 

. 09003 


PTRIA2 

155 

1 

. 1 4073 


F-'TRT A7 

1 f=;A 

1 

. 1 4R47 
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• PTRIA2 

157 

1 

. 17940 

- PTRIA2 

158 

1 

. 17887 

PTRIA2 

159 

1 

.20047 

PTRIA2 

160 

1 

. 19430 

PTRIA2 

161 

1 

.20483 

PTRIA2 

162 

1 

. 19300 

PTRIA2 

163 

1 

. 18800 

PTRIA2 

164 

1 

. 17083 

PTRIA2 

165 

1 

. 12680 

PTRIA2 

166 

1 

. 14623 
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77 

77 

92 

92 

92 

92 

92 

92 

108 

108 

108 

108 

108 

108 

37 

37 

37 

37 

37 

37 

49 

49 

49 

49 

49 

49 

60 

60 

60 

60 

60 

60 

70 

70 

70 

70 

70 

70 

82 

82 

82 

82 

82 

82 

18 

18 

18 

in 


3 

3 

1 

1 

2 

2 

3 

3 

1 

1 

2 

2 

3 

3 

1 

1 

2 

2 

3 

3 

1 

1 

2 

2 

3 

3 

1 

1 


3 

3 

1 

1 

2 
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3 

3 

1 

1 

2 

2 

3 

3 

1 

1 

2 

2 

3 

3 

1 

1 


3 

3 

1 

1 


0.95973878E-02 
“102. 15 
0.60110915E+00 
8.85 

0. 14405499E+01 
-171. 15 
0.54113474E-01 
8.85 

0.66993644E+00 
17.25 
0. 16054936E+01 
-162.75 
0.60309493E-01 
17.25 
0.371 52 123E+00 
35. 10 
0.89034558E+00 
-144.90 
0.3344S348E-01 
35. 10 
O. 207353 13E+00 
53.41 
0.49691896E+00 
“126.59 
0. 1 866649 lE-01 
53.41 
0.55394205E-01 
68.02 
0. 13275146E+00 
-111.98 
0.49867364E-02 
68.02 
0.7981 79 15E+00 
25.97 
0. 14759277E+01 
-154.03 
O. 200081 19E+00 
25.97 
0.71490322E+00 
31.01 
O. 13219407E+01 
“148.99 
O. 17920624E+00 
31.01 
O. 27469639E+00 
40.23 
O. 507946 15E+00 
“139.77 
0. 6885870 lE-01 
40.23 
O. 17440697E+00 
49.76 
O. 322499 13E+00 
-130.24 
0.43718948E-01 
49.76 
0.67720352E-01 
59.47 
0. 12522295E+00 
-120.53 
b. 16975598E-01 
59. 47 
0.40443266E+00 
34.14 
O. 63937426E+00 
-145.86 
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DAKEA* 

1 1 




u. 

DPHASE* 

12 


18 

3 

34. 14 

DAREA* 

11 


21 

1 

0, 61634525E+00 

DPHASE* 

12 


21 

1 

36.35 

DAREA* 

11 


21 

2 

0.97439037E+00 

DPHASE* 

12 


21 

2 

-143.65 

DAREA* 

11 


21 

3 

0. 200502 13E+00 

DPHASE* 

12 


21 

3 

36.35 

DAREA* 

11 


33 

1 

0.38572686E+00 

DPHASE* 

12 


33 

1 

41.14 

DAREA* 

11 


33 


0.609B0195E+00 

DPHASE* 

12 


33 

2 

-138.86 

DAREA* 

11 


33 

3 

0. 12548009E+00 

DPHASE* 

12 


33 

3 

41. 14 

DAREA* 

1 1 


44 

1 

0. 1 62657 19E+00 

DPHASE* 

12 


44 

1 

47.73 

DAREA* 

11 


44 

2 

0.25714743E+00 

DPHASE* 

12 


44 

2 

-132.27 

DAREA* 

11 


44 

3 

0.52913708E-01 

DPHASE* 

12 


44 

3 

47.73 

DAREA* 

11 


54 

1 

0. 366740 13E-01 

DPHASE* 

12 


54 

1 

.56.96 

DAREA* 

11 


54 

2 

0.57978552E-01 

DPHASE* 

12 


54 

2 

-123.04 

DAREA* 

11 


54 

3 

0. 11930355E-01 

DPHASE* 

12 


54 

3 

56.96 

DAREA* 

11 


1 

1 

0.30489490E+00 

DPHASE* 

12 


1 

1 

30.96 

DAREA* 

11 

ORiQiNAL PAGE ES 

1 

2 

0. 43567597E-^00 

DPHASE* 

12 

1 

2 

-149.04 

DAREA* 

11 

OF POOR QUALfTY 

1 

3 

0.99471353E-01 

DPHASE* 

12 

1 

3 

30.96 

DAREA* 

11 


13 

1 

0.23555727E+00 

DPHASE* 

12 


13 

1 

33.94 

DAREA* 

11 


13 

2 

0.33659678E+00 

DPHASE* 

12 


13 

2 

-146.06 

DAREA* 

11 


13 

3 

0. 76850090E-01 

DPHASE* 

12 


13 

3 

33.94 

DAREA* 

11 


14 

1 

0.7861 237 lE-01 

DPHASE* 

12 


14 

1 

38,54 

DAREA* 

11 


14 

2 

0. 11233222E+00 

DPHASE* 

12 


14 

2 

-141.46 

DAREA* 

11 


14 

3 

0. 25647 129E-01 

DPHASE* 

12 


14 

3 

38.54 

DAREA* 

11 


15 

1 

0. 10624694E+00 

DPHASE* 

12 


15 

1 

43.11 

DAREA* 

11 


15 

2 

0. 15182032E+00 

DPHASE* 

12 


15 

2 

-136.89 

DAREA* 

11 


15 

3 

0.34662853E-01 

DPHASE* 

12 


IS 

3 

43. 11 

DAREA* 

11 


27 

1 

0.48804804E-01 

DPHASE* 

12 


27 

1 

49.89 

DAREA* 

11 


27 

2 

0.69739050E-01 

DPHASE* 

12 


27 

2 

-130.11 

DAREA* 

11 


27 

3 

0. 15922470E-01 

DPHASE* 

ENDDATA 

12 


27 

“T 

49.89 


NOTE; ALL DATA FROM AERO* CARD UPTO ENDDATA CARD ARE 
OUTPUT FROM AIRLOADS PROGRAM 
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